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MEETING OF 1958 OCTOBER 10 
Dr W. H. Steavenson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 

John Evan Baldwin, Cavendish Laboratory, Cambridge (proposed by 
F. Graham Smith); and 

Norman Arthur Routledge, Department of Geodesy and Geophysics, 
Cambridge (proposed by M. N. Hill). 


The election by the Council of the following Junior Member was duly 
confirmed :— 
James Martin Hall, 70 Bishop’s Road, London, S.W.6 (proposed by J. M. 
Bruckshaw). 


One hundred and forty-five presents were announced as having been received 
since the last meeting, including :— 
G. de Vaucouleurs, L’exploration des galaxies voisines (presented by the 
author); and 
J. W. Dungey, Cosmic electrodynamics (presented by the author). 


MEETING OF 1958 NOVEMBER 14 
Dr W. H. Steavenson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


Frederick Donald Abbott, Dept. of Astronomy, The University, St. Andrews, 
Scotland (proposed by E. Finlay Freundlich); 

Reginald Edds, 1o Orme Court, London, W.2 (proposed by H. C. King); 

Harold Hicks, 34 Kingsley Road, Plymouth, Devon (proposed by E. A. Beet); 

Joachim Christoffel Cornelius van Loggerenberg, 3b Abraham Grayling 
Street, Wilgehof, Bloemfontein, S. Africa (proposed by A. D. Thackeray); 

Cecil Newman-Sanders, The Little House, Hatch End, Middlesex (proposed 
by H. C. King); 

Derek William Sida, Dept. of Mathematics, The University, Leeds (pro- 
posed by T. G. Cowling); 

Julius Dirk Willem Staal, Quarterdeck, ‘T'wo-mile-Ash, Horsham, Sussex 
(proposed by H. C. King); and 

Bengt Strémgren, Institute for Advanced Study, Princeton, New Jersey, 
U.S.A. (proposed by H. Bondi). 
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The election by the Council of the following Junior Members was duly 
confirmed :— 


Dennis Keith Richards, 9 High Holborn, Sedgley, Worcs. (proposed by 
H. E. Howle); and 

Robert Reginald Shobbrook, 75 Romway Road, Leicester (proposed by 
E. Finlay Freundlich). 


Eighty-five presents were announced as having been received since the last 
meeting, including :— 
R. S. Ball, The story of the heavens (presented by Col. J. A. Watson); 
J. C. Begg, Essays on thoughts and worlds (presented by the author); and 
International Astronomical Union and International Union of ‘Theoretical 
and Applied Mechanics, Proceedings of third symposium on cosmical gas 
dynamics (presented by the I.A.U. and the I.U.T.A.M.). 


MEETING OF 1958 DECEMBER 12 
Dr W. H. Steavenson, President, in the Chair 


The President announced the death of Dr John Jackson, a former President 
of the Society, and paid a tribute to his memory, the Fellows standing. 
The election by the Council of the following Fellows was duly confirmed :— 


Geoffrey Basil Anderton, 10 Taybank Road, Westering, Port Elizabeth, 
S. Africa (proposed by H. Welsh); 

Arnold Noel Argue, The Observatories, Madingley Road, Cambridge 
(proposed by D. E. Blackwell); 

Leslie Barraclough, 52 School Street, Low Moor, Bradford (proposed by 
G. Fielder); 

Cecil Maxwell Cade, 42 Latton Green, Commonside Road, Harlow, Essex 
(proposed by H. Spencer Jones); 

Barry Arthur James Clark, 11 Ryan Street, Coburg N. 13, Victoria, Australia 
(proposed by W. E. L. Clapham); 

Brian Clegg, 29 Primley Park Crescent, Leeds 17 (proposed by G. Fielder) ; 

William Davidson, Battersea College of Technology, Battersea Park Road, 
London, S.W.11 (proposed by W. H. McCrea); 

Peter Donnelly, King’s Arms Hotel, Wigton, Cumberland (proposed by 
W. L. Rae); ' 

George Alfred Eveleigh, 5 Hale Road, Walton, Liverpool 4 (proposed by 
J. R. Platt); 

Pierre Gouin, University College Geophysical Observatory, Addis Ababa, 
Ethiopia (proposed by D. O’Connell); 

Ronald Green, 12 Masson Street, Turner, Canberra, Australia (proposed 
by A. W. Rodgers); 

Ernest David Jones, Alba S.A., Caixa Postal 738, Curitiba, Parana, Brazil 
(proposed by C. A. Ronan); 

Peter James Jooste, 58 Strathleven Road, London, S.W.2 (proposed by 
H. Wildey); 

Raymond Leopold Kerdel, St. Peter’s Rectory, Broken Hill, N.S.W.., 
Australia (proposed by C. S. Middleton); 
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Otto Ernst Walter Rudolph Kurth, Dept. of Astronomy, The University, 
Manchester 13 (proposed by Z. Kopal); 

John Wishart Macvey, 15 Adair Avenue, Saltcoats, Scotland (proposed by 
A. E. Roy); 

James Gault Mitchell, c/o Marischal College, Aberdeen, Scotland (proposed 
by R. V. Jones); 

Thomas William Olle, P.O. Box 174, ‘The Hague, Netherlands (proposed by 
Z. Kopal); 

Miroslav Plavec, Astronomical Institute, Czech Academy of Sciences, 
Ondrejov, Czechoslovakia (proposed by Z. Kopal); 

James Dickie Raeside, D.S.I.R., P.O. Box 733, Dunedin, New Zealand 
(proposed by A.F.A.L. Jones); 

Richard Edwin Ridley, 23 Woodlands Avenue, Wolstanton, Newcastle, 
Staffs (proposed by ‘T. E. Sadler); 

Richard Langley Sears, Goethe Link Observatory, Indiana University, 
Bloomington, Indiana, U.S.A. (proposed by J. B. irwin); 

Uli W. Steinlin, Astronomisch-Meteorologische Anstalt, Universitat Basel, 
Binningen, Switzerland (proposed by W. P. Bidelman); 

Peter Antony Eabry Stewart, 11 Fairhaven Road, Bristol 6 (proposed by 
A. F. Collins); 

John H. Waddell III, Upper Air Research Observatory, Sacramento Peak, 
Sunspot, New Mexico, U.S.A. (proposed by J. T. Jefferies); and 

Arne A. Wyller, Sproul Observatory, Swarthmore College, Pa., U.S.A. 
(proposed by M. K. Vainu Bappu). 


‘The re-election by the Council of the following Fellows was duly confirmed :— 

Colin Kerr Grant, c/o Bank of Adelaide, 11 Leadenhall Street, London, 
E.C.3 (proposed by-T. F. Gaskell); 

David George Kendall, Magdalen College, Oxford (proposed by H. H. 
Plaskett); and 

Cecil Francis Retter, 40 Parkhurst Road, Torquay (proposed by T. A. 


Brown). 
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THE SCATTERING OF RADIO WAVES IN 'THE SOLAR CORONA 


A. Hewish 
(Received 1958 July 29) 


Summary 


Observations of the radio emission from the Crab nebula during June, 
when it passes within a small angular distance of the Sun, have shown that the 
radiation suffers considerable scattering during its passage through the 
solar corona. In this paper an account is given of measurements carried 
out each year over the period 1952-1958. ‘The following new results have 
been obtained: 

(i) measurements at a wave-length of 7-9 m have suggested the presence 
of refraction effects, in addition to scattering, which may provide information 
about the mean electron density in the corona; 

(ii) extension of the measurements to the shorter wave-length of 1-9 m 
has enabled the large scattering at the distance of closest approach to be 
determined more accurately; 

(iii) a pronounced sunspot-cycle variation has been observed in regions 
of the corona at a distance of 8R@; and 

(iv) the scattering has been shown to exhibit anisotropy. This result 
strengthens the earlier conclusion that the scattering is caused by filamentary 
irregularities aligned in a magnetic field and suggests that the direction 
of the field is approximately radial at distances of 15-20R (>. 





1. Introduction.—The possibility of investigating the solar corona by studying 
the refraction effects which arise when a radio star is viewed through it 
was first tested experimentally in 1951 and 1952 (Machin and Smith 1952, 
Vitkevitch 1958). In these experiments the Crab nebula was observed during 
June when it passed within 1°-5 of the Sun’s southern limb. The measurements 
were made with phase-switching interferometers aligned on an east-west axis 
and the recorded amplitude of the radio source decreased gradually as it 
approached the Sun, the reduction first becoming apparent at an angular 
distance of about 5°. 

The results could not be explained by simple refraction or absorption, 
but were consistent with a scattering process in which the apparent diameter 
of the source was increased to such an extent that it became partially resolved 
by the interferometer (Hewish 1955). ‘The spread of the scattered waves was 
found to be considerably greater than the angle through which the radiation 
would be deviated by the radial gradient of electron density and the scattering 
therefore masks any simple refraction effects. 

The scattering was explained by the presence of irregularities in the corona 
having a scale which could lie in the range 1-105 km. It was suggested that the 
irregular structure was maintained by a general magnetic field in a similar 
manner to the visible polar plumes (van de Hulst 1950). 

In the present paper an account is given of more recent observations which 
have been made each year since 1953. In addition to confirming the earlier 
work a number of new results have been obtained (see Summary, (i) to (iv), 
above). 

2. The observations.—All observations have been of the Crab nebula during 
the period June 6-23 when it passes within 4-6Ro of the Sun’s southern limb. 
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Phase-switching interferometers were used at a variety of spacings and 

A wave-lengths in order to study the scattering over as wide a range of angular 
distances as possible. Observations at long wave-lengths and with large 
interferometer spacings are necessary to study the small scattering in the 
outermost corona, while observations at shorter wave-lengths must be used 
to measure the largest scattering at the distance of closest approach. Interfero- 
meters employing a spacing of only a few wave-lengths, for which the scattering 
is always much smaller than the angular separation of the interference lobes, 
have also been used to study changes of integrated intensity of the source. 
While most of the interferometers have been aligned on an east-west axis, 
measurements were made in 1956 using NS and EW axes and in 1958 using 
axes in NE, NW and EW directions. Also in 1958 observations were made 
using a 3 000 ft array which gave a fan beam of width 1° to half power at 7-9 m. 
Use of the fan beam enabled the scattering to be determined directly from the 
broadening of the response. ‘The 1958 measurements, which were carried out 
in conjunction with J. A. Hégbom, will be discussed at greater length elsewhere. 
Some preliminary results are included in the present account as they remove 
an ambiguity concerning the direction of the scattering anisotropy. 

Rejection of radio emission from the Sun itself, which is considerably more 
intense than that from the Crab nebula except at the longest wave-length used, 
represents the main experimental difficulty. ‘The use of interferometers greatly 
reduces the response due to the quiet Sun and the most precise observations 
were carried out with the large four-aerial interferometer, first at 3-7 m and 

: TABLE I 
25 Wave-length Interferometer Aerial system 
ng spacing 
its 1952 79m ; 49 A EW 150 ft x 16 ft corner reflectors 
nie 3°7m 105A ,, 40 element linear array 
it 1953 79m 49A ,, As in 1952 
lar ms Sa. 5 150 ft x 16 ft corner reflectors 
3°97 m TOs A, As in 1952 
yn, _- is7A + Four-element radio telescope 
ter 1954 79m CH x As in 1953 
ed ” 49A ,, 9” 
vas 3°7m 157A ,, moo” 
ion 1955 79m SA. x As in 1953 
ing s SOR . 5. ‘Twelve-element Yagi arrays 
1956 79m 8A ,, As in 1953 . 
yna - 52A NS 100 ft x 30 ft corner reflectors 
the 5 56 A EW As in 1955 
lar 377m TOS A, As in 1953 
19m 300A ,, Four-element radio telescope 
ich 1957 I'9m 300A ,, As in 1956 
lier 1958 79m 146.4 x Six-element Yagi arrays 
iv), _ 43 A NW Yagi array plus 3000 ft array 
- 52 A NE - o» ” ” ” 
ing 9% +}° fan beam 3000 ft array plus 200 ft x 40 ft 
nb. corner 
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later at 1-9m. ‘This system allowed the solar response, in the absence of 
sunspots, to be reduced to about 1 per cent of that of the radio source. 
Details of the aerial systems used for the different experiments are 
summarized in Table I, and some typical records are shown in Fig. 1. 
CRAB NEBULA 
11:56:56 
1335 GMT 
19m 300) 
SUN 
t t A NI P 13:6 56 
1305 1310 1315 1320 1325 GMT 
SUN 
4 
casieatts neliaiaaimedneltt: tail dk jh tates Militant 13: 6:56 
7.9m 56% 0500 0530 | 0600 0630 ST SUN 
po Vt 4 ssa 
_ (PPP 
SUN 
‘ 
tK, 13 6 54 
ieee eee 
nee 0500 0530 0600 0630 5.7 
e™ shibeg 
ee h a 
SUN 
7.9m {on becom é ‘ - HA 
; > a 1: 6:58 
“ial ee ef 
| | | Ie 3 | 
0510 0520 0530 : ‘ 0950 $1 
7 Sy : ‘ 20 
ae as 5, eee 17.6558 
Sun 
Fic. 1.— Some typical records obtained with the different systems during the 
course of the 1956 occultation, and two fan beam records obtained in 1958. 
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Fic. 2.—The normalized response of the Crab nebula as observed with interferometers 
of different spacing and wave-length during the occultations of 1952—1958. 
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Fic. 3.—The position of the Crab nebula relative to the Sun at noon. 
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‘The basic results are illustrated in Fig. 2 in which the apparent intensity 
of the source is plotted as a function of the date of observation. ‘The position 
of the Crab nebula on the corresponding dates at noon is shown in Fig. 3. 
The general similarity of the occultation curves from year to year is immediately 
apparent, although greater reductions were observed towards sunspot maximum. 
It is also apparent that the reductions increase both with the wave-length, 
and with the interferometer spacing at a given wave-length. The detailed 
interpretation of these curves will be considered in the following section. 


3. The interpretation of the observations 

(a) The scattering theory.__It was shown previously that the results obtained 
in 1953 were explained by a scattering mechanism in which the integrated 
intensity of the source was unchanged, but the radiation emerging from the 
corona could be represented by an angular spectrum of the form exp (—¢4?/¢%)). 
The half-width (¢)) of the angular spectrum depends upon the square of the 
wave-length, the size of the irregularities, and on the deviation of the electron 
density from its mean value. ‘The present results have entirely confirmed this 
theory, which is also in agreement with other observations (Vitkevitch 1955, 
Blum and Boischot 1957). A contrary hypothesis has been suggested by Slee 
(1956) corresponding to a mechanism previously considered (Hewish 1955, 
p. 244, Case I) and shown to disagree with the observations. ‘This point has 
been discussed elsewhere (Hewish 1957). 

The amplitude of the record obtained when radiation having an angular 
spectrum given by exp —(¢?/4%)) is incident upon a phase-switching interfero- 
meter of spacing d and wave-length A is proportional to exp —(z?d?%,/A?). 
A comparison of the occulted and unocculted responses of the radio source 
enables ¢, to be determined. 

The values of ¢) derived from the different observations are shown ir 
Fig. 4, plotted against angular distance from the Sun. It is interesting to note 
that in 1956 and 1958 the scattering along an EW line was less than that 
along NS, NE, and NW lines*. In 1958 the largest scattering at closest 
approach was determined from the fan beam records and it was interesting to 
find that the results at greater angular distance agreed well with the 
interferometric data. 

In Fig. 5, values of log dy, obtained simultaneously at the different 
wave-lengths during 1956, have been plotted against log A. It is seen that the 
points lie reasonably well on straight lines of slope 2. Such a result, which 
is not accounted for by Slee’s hypothesis, confirms the scattering theory over a 
4:1 range of wave-length corresponding to a 16:1 variation of scattering 
angle. 

(6) Modifications to the simple scattering theory at long wave-lengths.—The 
simple scattering theory assumes that the scattered distribution is isotropic. 
When the scattering is so great that the scattered distribution extends over an 
appreciable range of radial distance then the shape of the distribution depends 
upon the radial gradient of scattering. ‘This effect must certainly occur at 
7-9 m at the distance of closest approach when the scattered distribution extends 
over at least two solar diameters. Some computations for this wave-length 
have been made on the assumption that the scattering is isotropic and follows 


* Only the NE scattering has been plotted in Fig. 4. The NW scattering was similar and has 
been omitted to avoid confusion. 
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a radial law with spherical symmetry. ‘The degree of scattering was obtained 
from the 1-9 m observations of 1956, appropriately scaled, since at this short 


wave-length such effects are believed to be negligible. 


Brightness contours 


were computed for the scattered distribution at angular distances of 6Ro and 


12Ro and are shown in Fig. 6. 


It is seen that the scattering is only slightly 


anisotropic at 12Ro but significantly so at the shorter distance. 
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Fic. 4.—The scattering half-width $, measured at different wave-lengths plotted as a function 


of angular distance from the Sun. 


The arrows denote lower limit. 
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Fic. 5.—The scattering half-width ¢, measured at three wave-lengths simultaneously during 
1956 plotted as a function of the wave-length on a logarithmic scale. The broken line indicates 
a slope of 2 as predicted by the scattering theory. 
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Fic. 6.—The theoretical scattered distribution at 7-9 m computed for isotropic scattering 
in a spherically symmetrical corona. 


(c) Evidence for refraction.—Observations at 7-9 m using an interferometer 
spacing of 8 A exhibit marked decreases of integrated intensity when the source 
lies within an angular distance of 8Rg. ‘The relevant observations have been 
re-plotted in Fig. 7 to demonstrate the effect more clearly. ‘These occultation 
curves can only be accounted for by a real decrease of integrated intensity 
since the scattered distribution would otherwise need to subtend an angle 
greater than 6°; it is known from measurements at’ the same wave-length at 
greater angular differences (but still within + 3°) that the scattering irregularities 
are not distributed widely enough to produce such an effect. 
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Fic. 7.—The reduction of integrated intensity observed at 7:9 m compared with theory. 


While the simple theory of scattering gives no change of integrated intensity, 
some reduction will occur when the angle of scattering away from the direction 
of incidence is greater than the angle subtended at the Earth by the scattering 
region. An estimate of this effect, based on the observed scattering at closest 
approach at 1-9 m in 1956, is shown in Fig. 7. It appears that this mechanism 
alone cannot entirely explain the larger decreases actually observed. 

It is possible, however, that some reduction of intensity might be caused 
by refraction phenomena arising from the radial gradient of the mean electron 
density. In the absence of scattering, refraction would cause a divergence 
of the incident rays away from the Sun’s centre, which would give rise to a 
caustic surface. If the Earth passed through this caustic a complete occultation 
would occur, preceded by a sharp increase of intensity (Link 1951). Calculations 
based on the electron densities quoted by Blackwell (1956) show that the 
Earth would pass through the caustic when the source was at an angular distance 
of about 7Ro. In fact, as shown in Fig. 8, the observed scattering is appreciably 
greater than the computed angle of refraction at distances exceeding 5Ro and so 
a straightforward occuitation of this type is not possible. At smaller angular 
distances the refraction phenomena may well predominate, however, owing 
to the steeper radial variation of the refraction angle. Occultation of a portion 
of the scattered distribution nearest the Sun is then a possibility. ‘The precise 
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reduction of intensity which might arise cannot be estimated without more 
knowledge of the form of the scattered distribution close to the Sun, but it 
seems difficult to account for the reductions shown in Fig. 7 without increasing 
the refraction angle. A two-fold increase of the mean electron density would 
make the refraction angle and the scattering angle comparable at about 6Ro 
and a sufficiently great reduction might then occur, but this question must 
remain open in the absence of further data. 











= | L 
5 10 15 Ro 
Fic. 8.—The computed angle of refraction (8) due to the radial gradient of the mean electron density 
compared with the measured scattering angle at 7-9 m (1956). 
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Fic. 9.—The solar cycle variation of scattering at 7-9 m. 
The point for 1957 was derived by scaling a measurement at 1-9 m. 


(d) Solar cycle variations.—In Fig. g the scattering at an angular distance 
of 8Ro has been plotted as a function of the solar epoch. It is seen that the 
scattering falls to its least value in 1953-54, the period of sunspot minimum, 
and increases rapidly to a maximum in 1957-58. At an angular distance of 16Ro 
no such variation is apparent and there seems to be no significant change during 
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the solar cycle. Since at 16Ro we are concerned with scattering in the equatorial 
region, and at 8Ro with regions nearer the pole, these results suggest that the 
distribution of the scattering irregularities follows the changing form of the 
visible corona, which tends to become more spherical towards sunspot 
maximum. 
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Fic. 10.—The scattering computed for a dipolar model of the general magnetic field compared 
with observations carried out in 1956 and 1958 at 7-9 m. The arrows denote lower limits of the 
observed scattering. 

(e) The scattering anisotropy.—The observations carried out in 1956 and 
1958 (outlined in Section 2), using interferometers aligned in different directions, 
showed that the scattering was more pronounced in a direction parallel to the 
solar axis. In 1956 only two axes were used so that the observations could 
either be interpreted as a scattering of the radiation over an elliptical area, 
or simply. along a line in a direction intermediate between the interferometer 
axes. This ambiguity was removed in 1958 by the use of three axes and some 
ellipses defined by the three scattering angles are shown in Fig. 10*. The 


* It is likely that the ellipses for other days will be known after a more complete analysis of 
the 1958 data has been carried out. 
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results for the two years are seen to be quite compatible since, for 1956, any 
ellipse which has the observed limits of scattering in the two directions as 
tangents represents a possible interpretation. 

It was shown in Section 3 (4) that anisotropy arising from the radial gradient 
of scattering was small at angular distances exceeding 12Ro and it follows that 
the irregularities themselves must scatter anisotropically. 

The origin of coronal irregularities whose size is smaller than the mean 
free path raises difficulties unless it is assumed that they are maintained by 
a magnetic field (Hewish 1955). Ifa general magnetic field pervades the outer 
corona to distances of at least 20Ro, then material will diffuse rapidly along the 
direction of the field but not across it, and it is reasonable to suppose that this 
process will produce a filamentary structure similar to the visible polar plumes. 
Any individual filament will scatter perpendicularly to its axis and it is likely 
that the integrated scattering along a line of sight will show anisotropy 
dependent upon the mean direction of the magnetic field. It is shown in the 
next section that such a model accounts for the observations, provided that 
the lines of force are more radial than those of a dipole field. 

4. Scattering by irregularities aligned in a magnetic field—The present 
observations are not sufficient to deduce, uniquely, the precise shape and 
distribution of the scattering irregularities throughout the corona. It is, however, 
possible to suggest a type of coronal model in reasonable agreement with the 
data. Such a model must account both for the variation of scattering with 
angular distance from the Sun and for the shape of the scattered distribution. 
In computing the scattering the following quantities are relevant: (i) the size 
and shape of the irregularities; (ii) AN/N, the variation of the electron density 
from its mean value; and (iii) the radial variation and latitude dependence 
of both (i) and (ii). 

In the model adopted it has been assumed (a) that the irregularities are 
aligned by dipole magnetic field in the range 5-20Ro, (6) that the mean equatorial 
electron density is that given by Blackwell (1956) and that the contours of 
constant density are ellipses of axial ratio 3: 1, (c) that the value of AN/N 
is constant and (d) that the average size of the irregularities is independent 
of latitude near the Sun’s surface. The regions of higher and lower density 
thus take the form of slightly conical filaments along the field lines whose area 
of cross-section is approximately the same at their base. 

Since, as we suppose, the irregular structure arises from diffusion along the 
field of density variations generated at the base of the corona, conditions (a) 
and (b) may not be quite independent. The radial variation of density along 
visible streamers is, however, known to be more rapid than can be accounted 
for by an outward flow alone (van de Hulst 1950) and it seems not unreasonable 
to assume that the shape of the filaments is determined by the magnetic field, 
while their density is controlled by other factors. 

To compute the scattering for various lines of sight the model corona was 
divided into ten concentric shells, the bounding surface between adjacent 
shells being swept out by the rotation of a particular loop of force about the 
dipole axis. ‘The scattering was assumed to be constant along each element 
of a line of sight within a given shell, and its value was estimated from the 
mean density, cross-sectional area and orientation of the filaments. Finally, 
the contributions due to each shell were integrated to yield the approximate 
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elliptical distributions shown in Fig. 10, where the lines of sight correspond 
to 1956 June 10-14*. Since neither AN/N, nor the cross-section of the 
irregularities, is known absolutely only relative values of the scattering may 
be derived. Modifications due to the radial gradient of scattering have also 
been ignored, which implies that the results are not applicable to a wave-length 
of 7-9m for angular distances less than about 10Ro (cf. Section 3 (4)). The 
latter effect will tend to distort the ellipses into a crescent form. 

Comparing these theoretical results with the observations of 1956 and 
1958, it is seen that the eccentricity of the ellipses agrees reasonably well; their 
inclination, however, is not in agreement with the expectations of a dipole model. 
The observations suggest that the field lines extend more or less radially with 
no tendency to turn into the equatorial plane. 

In order to check the radial variation of scattering the relative scattering 
in an EW direction only has been plotted as a function of radial distance. 
This curve has been fitted to the experimental results for 1956 as shown in 
Fig. 11. The experimental curve is a composite one in which observations 
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Fic. 11.—The radial variation of EW scattering observed in 1956 
compared with that predicted for a dipolar model. 


at all wave-lengths have been scaled to a wave-length of 3-7m. It is seen 
that there is fair agreement although there is no evidence for the flattening 
of the theoretical curve at about 10oRq. This flattening arises from the rotation 
of the axis of the ellipse and its absence again suggests that the magnetic field 
does not bend towards the equatorial plane as must a dipole field. 

It is concluded that the type of model discussed, with some modification of 
the magnetic field, gives a fair representation of the observed scattering. It 
seems likely that more detailed measurements will provide valuable information 
on the direction of the solar magnetic field at great distances from the Sun. 

It should also be noted that the model discussed is one in which material 
flows outwards from the base of the corona. The alternative possibility, in 
which the outer corona consists of matter flowing inwards, makes it difficult 
to account for the present observations. 

* For a line of sight crossing the polar axis the projected direction of the filaments is entirely 
north-south and purely east-west scattering might be expected. This does not occur, owing to 


the strong isotropic scattering from the particular filaments tangential to the line of sight in front 
of, and behind, the Sun. 
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MEASUREMENT OF THE PROFILE 
OF A LUNAR WRINKLE RIDGE 


Gilbert Fielder 


(Received 1958 June 17) 


Summary 


A lunar wrinkle ridge has been measured uniquely by means of a method 
which involves measuring the changing length of a shadow cast across it by 
a neighbouring mountain. Probable errors in the positions of points are 
+330 m in the plane of the surface, and only +14 m along normals to the 
surface. The profile of the wrinkle ridge shows that it cannot be regarded 
as a wave front of solidified lava. 





1. Principle of the method.—\n essentials the method is that which McMath, 
Petrie and Sawyer (1937) proposed and used to obtain cross-sections of parts 
of Theophilus. Cine photographs are taken of the surrounds of a lunar mountain 
which is casting a shadow either at local sunrise or at local sunset. ‘The length 
of the shadow is measured on every good negative and, the instant of each exposure 
being known, it is possible to reduce the measurements to differences in vertical 
height of the country across which the shadow-tip travels. 

2. Equipment employed.—The telescope used for these observations was the 


f/30, 60 cm refractor of the Observatory of the Pic-du-Midi. At the prime 


focus, the image of full Moon is about 16 cm in diameter, and the depth of focus 
is about 1 mm at A=560 mp. A passband of only about 150A, which centred 
on this wave-length, was necessary partly because of the characteristics of the 
film which was selected and partly because the secondary spectrum of the 
objective is such that the longitudinal aberration is zero at 585 my, and the 
curve rises rather steeply to either side of this point, so that to have utilized 
light of shorter wave-length than 550 my, or of greater wave-length than 630 mp, 
would have necessitated the use of even smaller passbands and exposures increased 
beyond reasonable limits. 

Since the experiments were concerned essentially with photography of the 
terminator regions, where contrasts are very low, it was necessary to select an 
emulsion having as high a contrast as possible, compatible with 4 grain size 
smaller than the best seeing disk. ‘Tests showed that resolution was limited by 
scintillation, and not by grain-size. Of eight emulsions tested, each in conjunction 
with five different filters, four different developing agents, and differing times of 
development at constant temperature, the best proved to be Ilford 5Ggr in spite 
of the rather long exposure time of 3 sec which had to be employed. ‘This was 
used in conjunction with the Chance glass filter OY2 (yellow) and was developed 
for 40 minutes in Microphen or Promicrol at 20 °C. 

A 35mm film holder for 30m of film was adapted to automatic lunar 
cine-photography. ‘The obscurator used was a self re-cocking Zeis Ikon tri-metal 
leaf shutter, which yielded a clear aperture 2-5 cm in diameter when fully opened. 


38 








548 Gilbert Fielder Vol. 118 


Elementary considerations showed that a sufficiently large rate of taking 
pictures was one per minute. A Ir.p.m. synchronous motor was therefore 
arranged to operate the shutter by means of a cam, shown in Plate 8. This motor 
also wound on the film. The cam closed and opened a microswitch which was 
wired in series with a solenoid and 24 volt d.c. supply. In the closed position 
of the shutter the solenoid core, which was linked mechanically with the shutter 
lever, was held out of the solenoid by the shutter spring. When the microswitch 
was closed, the core plunged into the solenoid, opening the shutter. The duration 
of the exposure depended only upon the pre-set position of a timing lever which 
was attached to the shutter; the speed was continuously variable from 1 sec up 
to 1/100 sec, and there was a setting for a ‘‘time’’ exposure. 

‘The smaller of the visible cog-wheels in Plate 8, fixed to the spindle of the 
motor and carrying the cam, had a quarter of its total number of teeth removed. 
Thus, in operation, the motor was running continuously but the driven (larger) 
cog-wheel was at rest for 25 per cent of the time; it was an easy matter to arrange 
the system so that the shutter was opened and closed only within the interval 
during which the film was stationary. In addition, the geometry of the cog-wheels 
and sprocket was chosen so that a suitable frame-spacing was obtained. 

The film-holder was mounted on a specially constructed optical bench of 
steel which was fixed rigidly to the eyepiece end of the telescope. Focusing 
was achieved with the aid of a reflex system by sliding the film-holder along this 
optical bench. The observer viewed the image of the Moon in the plane of the 
scratches on a clear strip of film placed in position for normal photography, 
behind a pressure plate in the film-holder. This was made possible by having a 
hole o-g cm in diameter in the pressure plate. The images were viewed with 
a compound microscope mounted on the lid of the film-holder. When the 
images of the scratches on the film and of the Moon were sharply defined and 
coincident, the film-holder was clamped rigidly to the optical bench. The 
film-holder was then loaded with a roll of film, and the assumption was made 
that the depth of focus was of sufficient magnitude to accommodate any change 
of focus during a night’s work. This assumption was not disproved. 

The time of the first exposure was noted and successive frames were spaced 
by nearly equal intervals of one minute of time. (Slight corrections to the 
synchronous motor had to be effected because of mains’ fluctuations.) A finder 
was used for guiding. 

3. Reductions of the film.—The instrument which was selected for measuring 
shadow lengths on tHe negatives of the film was a Mark II Automatic Recording 
Microdensitometer, manufactured by Joyce, Loebl and Company, of Newcastle- 
upon-Tyne. Orientation of each transparency was effected rapidly, since the 
translucent table which carried it could be rotated about a vertical axis as well 
as displaced in two horizontal mutually perpendicular directions. It was an 
easy matter to align a pointed shadow with the line of scanning. 

The measured shadow lengths were converted into seconds of arc, s._ It may 
be shown that, at a given instant, 








op Se, a -/ae A _a- . (1) 
cosec p—c 


where A (degrees) is the altitude of the Sun at a point P, on the Moon’s surface, 
defined by the selenographical coordinates (A, 8,) of the mountain peak which 














of 
1g 
‘is 
he 
y; 


th 
ne 
id 
he 
de 
ve 








XUM 


Montuty Notices oF R.A.S. 





Vor. 118, Prate 8 


Gilbert Fielder, Measurement of the profile of a lunar wrinkle ridge 














Ww 


by 
B, 


fr 


ar 





No. 6, 1958 Measurement of the profile of a lunar wrinkle ridge 549 


is casting the shadow, p is the equatorial horizontal parallax of the Moon, the 

quantity (cosecp—c) is the distance between the observer and the point P, 

¢ is the phase angle, and h is the difference in height of the top of the mountain 

peak and the tip of the shadow, measured in units of the Moon’s mean radius. 

The principal terms in c have been discussed in a paper by MacDonald (1931). 
Solving equation (1) for h, 


h= 5 (sina 55), (2) 
where 
R= 56207 sing 
~ cosecp—c 


4. Results and estimated errors.—The particular peak of concern in this paper 
is that near to the southern termination of the lunar Alps. It is designated B 
by Wesley and Blagg (1935) and has selenocentric coordinates Ap= + 6°-03, 
Bp = +31°:00, taken from the same 1.A.U. map. The wrinkle ridge in question 
is that which, at a mean distance of some 30 km east of peak f, runs roughly 
from north to south. The profile of the wrinkle ridge has been plotted in Fig. 1. 

The measured section, which cuts roughly east to west, indicates that there 
are minor ridges running along the top of the broad-based wrinkle ridge. The 
width of the base of the ridge is of the order of 20 km. Also, the country on both 
sides of the wrinkle ridge is lower than the crests of the ridge by at least too m. 
These facts seem to indicate that the wrinkle ridge is not a frozen wave of lava, 
as Baldwin (1949) has suggested. 

Profiles of the eastern slopes (broken lines) of this wrinkle ridge were not 
obtained because the hills themselves shielded the sunlight from these regions 
during sunrise. However, the western slopes of the ridge are seen to be vertical, 
in places, and undercuts in these parts are not excluded. 

An estimate of the most-probable diameter of the seeing disk gave ~ 0-33 arc. 
As a first approximation, therefore, it is thought reasonable to take the error in s 
to be 6s= +0":16 (i.e. +330 m). 

The estimated error in h is then given by 


ah= SH as+ % 544% oe, (3) 
From equation (2), 
0h .sinA_ s 
a EE (4) 
In equation (4) s/E? is essentially positive. Hence, 
oh sin A 
(5) ax “— (5) 


Since A and E are quantities which have been computed to the desired degree 

of accuracy, 5d =dE~o0, and using these conditions, together with condition 

(5), in equation (3), 

sin A 

E 

Using 5s= +0”:16, it follows that the estimated error in A does not exceed 

+14m. It should be noted that the assumption has been made that any 

systematic errors having their origin in the photographic emulsion have been 
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(3h) max = 8s. (6) 
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neglected, and hence the quoted error should be considered to be a provisional 
estimate, only. 

Peak 8 towered to a height of 1390 m above the arbitrary zero level shown 
in Fig. 1; but the exact radius of the arbitrary reference surface is not known 
and it is necessary to refer the height of the peak to some point on the ridge. 
‘The peak is 1290+14m above the highest point of the wrinkle ridge, and 
1440+ 14m above the lowest point measured. On subtraction, these figures 
give 150+28m for the greatest difference in the altitudes of points on and 
adjacent to the wrinkle ridge. 
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Fic. 1.—Profile of a wrinkle ridge. 
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A SHORT METHOD FOR THE DISCOVERY OF NEPTUNE 
R. A. Lyttleton 
(Received 1958 September 25) 


Summary 


‘The mathematical basis is described of a method for the discovery of 
Neptune that involves far less numerical calculation. The time of helio- 
centric conjunction can be found solely from considerations of the 
discrepancy in the longitude of Uranus. ‘This information alone makes 
prediction possible within less than 15° on the basis of Bode’s law. By 
finding the distance in a circular orbit appropriate to the best fit of the 
observations a prediction can be made comparable in accuracy with that 
achieved by Le Verrier. It is shown how by suitable combination of the 
equations of condition the number of unknowns can be reduced to three 
(as compared with eight in the original methods) for any assumed mean 
distance, and the same process removes certain awkward features that 
would otherwise enter for orbits near the 2 : 1 resonance. 





1. Introduction.—It is an outstanding question of interest in relation to the 
discovery of Neptune whether in fact there could have been devised any much 
less laborious method that would have achieved the prediction with an accuracy 
comparable with that attained by Adams (23°) and Le Verrier (1°). Such a 
method is described in the present paper as to its theoretical basis; the details 
of its numerical application have been given in full elsewhere*. The whole 
of the present investigation was inspired by discussions between the author and 
J. E. Littlewood following the 1946 centenary celebrations of the original 
discovery, and the procedure here described in Section 2 for. estimating the 
time of heliocentric conjunction, which is the first step of the present method, 
is due principally to himf. 

It is helpful first to explain briefly why the original methods were so extensive, 
and for this the one used by Le Verrier will be sufficient{. All perturbations by 
known planets having been dealt with, the unknowns that have to be solved for 
consist of the four corrections necessary to the existing elements assigned to 
the orbit of Uranus, da (or 5m), de, de, and 5a; the four elements of the (coplanar) 
orbit of Neptune, a’, e’, e’, and w’; and also the mass m’ of the unknown planet ; 
making apparently nine in all. ‘The quantities e’ and w’ can be replaced by 
h'=e' cosa’, k'=e' sina’, as is usual, and h’ and k’ then occur only linearly 
in the perturbations in longitude of Uranus§. But a’ and ¢’ do not occur linearly, 
while m’ multiplies everything throughout. The mean distance a’ was assigned 
(with considerable inaccuracy, it will be recalled) by means of Bode’s law, 
but e’ occurs trigonometrically as a e’, pg 2e’, etc. To determine a first 


* Vistas in Astronomy, Vol. III, 1959, Pergamon Press, London. 

+ For Prof. Littlewood’s own account of the problem see his A Mathematician’s Miscellany, 
pp. 117-134, Methuen, 1953. 

} For a more detailed account of Le Verrier’s method see: F. Tisserand, Mécanique Céleste, 
I, Ch. 23, Paris 1889; W. M. Smart, Celestial Mechanics, p. 259, Longmans, 1953. 

§ Tisserand, Mécanique Céleste, Vol. I, p. 378. 
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approximation to «’, with a view eventually to having linear equations, Le Verrier 
took in turn 40 values for e’ spaced at g° intervals covering the whole possible 
range from 0° to 360°, and for each value solved by least squares 18 consolidated 
equations of condition relating to 1690-1845, now in seven variables. It was 
then possible to select that value of e’ giving the best fit to the whole series of 
observations. Using this, he then proceeded to improve the solution to arrive 
finally at a set of ‘‘best’’ elements for Neptune. 

The quantity e’ is closely associated with the longitude of the unknown planet, 
and the first step of the method to be described in this paper leads to an equivalent 
but even more specific piece of information concerning the unknown planet, 
namely the time of heliocentric conjunction with Uranus. 


2. Determination of the time of conjunction.—The fact that the rate of change 
of angular momentum, dh/dt, must vanish at conjunction would enable the 
instant to be found were it not that h=r?dv/dt is not itself sufficiently well 
determined because of its strong dependence on the radius vector r. But an 
alternative procedure is possible that depends only on considerations of the 
longitude, and is one that involves little calculation in application. 

In undisturbed elliptic motion we have for the heliocentric longitude, in 
standard notation, 

v=nt+e+4+2esin(nt+e—m), (1) 


to the first order in the eccentricity. We note that for Uranus e=o0-047. If the 
elements are in error as a result of small perturbations depending on the mass m’ 
of the unknown planet (Sun=1), the resulting error in vw will be given by 


Av =tAn + Ae + 2Ae sin (nt + € — a) + 2e(Ae —Aw) cos (nt +e«—a) 
+ 2etAn cos (nt+e—wm), (2) 


and on the right-hand side, from consideration of the equations of motion, 
the first four terms are of order m’, while the last term is of order em’ and so 
is much smaller than the others. Hence if we agree to neglect this last term 
we can write, with an obvious notation, 


Av=m' (a) + bt +ccosnt + dsin nt) (3) 


for any discrepancy in longitude due to perturbations, correct to the present order. 

Consider now the elliptic orbit EF best fitting the observations of Uranus 
(after allowing for all known perturbations). The observations will not in fact 
be properly representable by E because they contain effects of the unknown 
body m’. Writing v,(t)= calculated heliocentric longitude of Uranus in ellipse E, 
v(t)=actual heliocentric longitude determined from observation, then the 
well-known discrepancies in longitude, which constitute the observational 
material to be explained, are given by 


6(t)=0(t)—2,(¢). (4) 

The values of 5(t) from 1690-1840 were the quantities utilized by Adams* and 

have been used in applying the present method. Le Verrier had observational 
data from 1690-1845. 

Next, denoting by ¢) the instant of heliocentric conjunction that we wish to 

find, let E, denote the instantaneous, or osculating, elliptic orbit of Uranus 


* Coll. Works, I, p. 11, Cambridge, 1896. 

















rh 


~a fn @. = 


a. 





XUM 


(8 
er 


le 


aS 
of 


ve 


1€ 





XUM 


No. 6, 1958 A short method for the discovery of Neptune 553 


corresponding to tj, and v(t) the longitude of Uranus in this instantaneous 
orbit. Also let the perturbations produced by m’ since time ty be denoted by w(t), 


so that 
w(t) = v(t) — v(t). 


3(t)=(vo— vy) +, (5) 
and in this both terms are of order m’. Since we are proposing to neglect all 
terms of order em’, then in calculating w(t) both orbits E and E, can be regarded 
as circular, and therefore w(t) will take equal and opposite values at times equally 
separated on the two sides of conjunction. So if 7 measures time from conjunction 
w(ty—7)= —w(ty+ 7), or writing 

W(r) = w(t) =w(ty + 7) (6) 


then W(7) is an odd function of 7. ‘The analytical form of W is given below 
by (11) but is not required for the present purpose of finding fy. 

Now to the first term on the right in (5), the form (3) will clearly be applicable, 
and so 


Then we have from (4) 


5(t) = m'{ay + bt) + br + ccos (nity +n7) +d sin (nty+n7); + W(r), 
which may be written, introducing new constants, in the form 
8(ty +7) = A+ B(1 —cosnr) +{Cr+ Dsinnr + W(7)}. (7) 
Since W(r) is an odd function, we also have 
8(t) — 7) =A + B(1 —cosnr) — {Cr+ Dsinnr + W(7)}. 


Hence, since 4(t))=A, it follows that 
neglecting terms of order em’ the function 
given by 
5(ty + 7) — 28(ty) + 8(ty — 7 
a(x) — Sltort7)=28(to) +3(l0~7) 
I —cosnt 


=B 
(8) 





will be constant. 

This result provides the following 
simple rule for determining the approxi- 
mate time of conjunction: Select an 
instant ty, calculate the second difference 
5(to + 7) — 28(ty) + (ty 7) for a number of 
values of +t, divide each by (1—cosnr), 
and the resulting quantity will be (approxi- 
mately) constant if ty has been selected at 
conjunction. 

The result of applying this rule to the 

é 1 = PTs 3 = Fic. 1.—Graph showing degree of 
observationa data 6(t) (Adams, loc. cit.) is constancy of p(7). The ordinate represents 
shown in Fig. 1. The quantity p(7) turns — the root mean value of (p—>)* plotted at ty 
out to be most nearly constant for a value > 1813, 1816, 1819, 1822, 1825, and 1828. 

f f fetes ics The different curves correspond to different 
0 to 0 about 1822-3. ‘This time for the weightings. For each the minimum occurs 
heliocentric conjunction in longitude differs at about 1822-3. 
by only just over six months from the + equal weightings of p(7); 


j P < weightings increasing with 7; 
actual value 1621°74. © weights increasing and then decreasing. 


— (rP--P/ 
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3. Prediction purely from a knowledge of conjunction.—On the basis of the 
unknown planet being in a circular orbit, it is a simple matter once the time 
of conjunction is available to calculate its longitude for any assumed size of 
orbit. Even on the crude assumption of Bode’s law, viz. a’/a=2, the longitude 
so predicted for the date of discovery is only about 13° behind the actual position. 
The actual planet would have been well inside a zodiacal belt 30° long by 10° 
wide centred on this place, which was the kind of region Airy suggested should 
be searched by Challis. 

The foregoing procedure can be carried out with a minimum of arithmetical 
labour, taking at most a few hours, but the question arises, in view of the serious 
failure of Bode’s law, which had it held would have made discovery almost 
triflingly easy by the present method, whether any better estimate of a’/a can 
be reached without unduly extensive calculation. For this purpose, Bode’s 
law is dropped and instead we assume always a circular orbit for the unknown 
body, which is at least an equally valid assumptivn to make for any planetary 
orbit. 

4. Equations for finding an improved value of a’'.—On this basis then, the 
perturbations in longitude of Uranus are given by* 


P(t)=—m' y F,sini{(n’ —n)t+«'—e} 
I 


+ a 
+m'e > G,sin [i{(n'’—n)t+e—€}+nt+e—a] (9) 
—o 
wherein 
(x2 Zz () 
es ee al 
2;°(1 —2; z?(1-32;7) da 


z,=1(1—n'/n), and A" is a certain infinite series in a/a’; while a similar form 
gives G;. Unless a definite numerical value is adopted for a/a’ there is clearly 
therefore no possibility of arriving at linearized equations, and this plainly must 
have been the great attraction of Bode’s law. 

The assumption of a circular orbit and knowledge of the time of conjunction 
immediately reduce the number of unknowns associated with Neptune to two, 
namely m’ and a’, and also remove altogether the awkward feature otherwise 
produced by the appearance of «’ trigonometrically, for once conjunction is 
known ¢ can be measured from it. Also since e is small, the terms in m’e in (9) 
are small compared with those solely in m’, and the perturbations can be 
calculated as if Uranus itself moved in a circular orbit. The principal terms 
in m’e can however readily be included for higher accuracy. Accordingly, if 
from now on the time is measured from conjunction, the whole expression for 
the discrepancies in heliocentric longitude will be of the form 


dv = de + thn + da sin nt + 88 cos nt + P(t) (10) 
wherein now 


P(t)= —m' > Fsini(n—n')t = —m' 5 F,siniD, (11) 
Il | 


say. If, for any assumed value of a/a’, direct comparison of (10) with the 
observations is made, the variables are now only five in number, de, 5, 5a, 58, 
and m’, already a considerable reduction on the number used by Adams and 


* Tisserand, Vol. I, p. 365. 
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Le Verrier, though practical application would require solution for a number 
of values of a/a’. But it is possible to eliminate the variables da and 5f, as we 
shall next explain. 


5. Further reduction of the number of unknowns.—In view of the unavoidable 
feature that a series of values of a/a’ must be taken, it is convenient to take advantage 
of a simple regrouping of the equations of condition that has the effect of removing 
da and 58, as well as a certain other important consequence. For this, we note 
the simple identities 

SIN (4 4+4)4+ SIM (x—9)— SIM y= (2 cosd—1) SM x (12) 
cos cos cos cos 
in which, if @=60°, the right-hand sides vanish. 

Now for Uranus 27/n=84-015 years, and accordingly 60° of longitude is 
described by the planet in almost exactly 14 years. Hence if instead of a particular 
observation (discrepancy in longitude), f(t) say, we adopt 


f(t+ 14) +f(t— 14)-f(t), 
where ¢ is now measured in years, then the terms in da and 6f will disappear. 

This device is of course no more than a regrouping of the equations of 
condition, incidentally reducing their number since, to some extent, 14 years 
are lost at each end of the range of available observations. But as will be shown 
later it simultaneously has the important effect of making the method safely 
applicable near the resonance n/n’=2 (for which a/a’=0-63 approximately) 
where in fact the coefficient F,, can become very large. 

Under this regrouping, the terms de + t6n come out the same, but the terms 
in P(t) have to be adjusted as to their coefficients. For the term in sin D, for 
instance, 

D+6=(n—n’')(t+14) where n.14=60 
=D+(1—n'/n)60°. 
So, writing v=n'/n, we have 6=(1—v)60', and hence to transform P(t) the 
coefficient F, has to be multiplied by 2 cos {(1 — v)60°}— 1, the coefficient F, has 


to be multiplied by 2cos{(1—v)120°}—1, and so on. With these conversion 
factors applied, we can write, denoting the new coefficients by /;, 


O(t)= P(t+14)+P(t—14)—P(t)=m' > f,siniD, (13) 
and the revised equations of condition take the form 
de + tn + O(t) = dv(t + 14) + v(t — 14) —dv(t), (14) 


and these, for each selected a/a’, involve now only three unknowns, namely de, 
dn, and m’. 

The procedure adopted was therefore to set up these modified equations of 
condition for each of a short series of values of a/a’, and solve by least squares 
to find what value of the ratio of axes gives the best fit to the observations. 


6. Effect of the transformation on the coefficient F,.—The following table 
shows the values of F,, F,, F;, F, for an appropriate set of values of a/a’. 
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TABLE | 


Coefficients of terms in P(t)=m'XF; siniD 


ala’ sin D sin 2D sin 3D sin 4D 
"50 — 5°738 1°429 o'l4I 0'031 
0°55 — 7°432 3°380 0'299 0068 
o'60 — 9911 15°712 0672 O'154 
0°65 —13°718 — 38-422 1692 0°366 
0°70 — 19°944 — 18-655 5°461 0°967 


It will be noticed that F, becomes very large in magnitude as the resonance 
n'|n=2, or a/a'=0-63, is approached from either side. ‘The expression (9), 
or (11), for the perturbations in longitude is not of course applicable for exact 
resonance, at which the coefficient F, would become infinite. Now to avoid 
excessive calculation of the various F; for different a/a’, it would obviously be 
convenient if we could interpolate within the above table, but in the case of F, 
interpolation from only a few values is not possible with much accuracy because 
of the resonance. (All these coefficients are nowadays available in tabular form*, 
but as they were not so to Adams and Le Verrier any comparison of methods 
must suppose the F; to be calculated as part of the work.) 

But the conversion factor 2cos@—1 vanishes at the resonance when the 
coefficient F, has its infinity through the factor 1 —;=2v—1 in its denominator. 
But it is easily seen that the expression 

2cos {(1—v)120°}—1. 
2v—I 





has limit 7/\/3 as v3}, and the coefficient f, passes smoothly through the 
resonance value without singularity, as the following table of the transformed 
coefficients shows: 

Tasce II 


Coefficients of terms in Q(t) =m' Xf; sin iD 


ala’ sin D sin 2D sin 3D sin 4D 
"50 — 3°208 —o'814 —0'266 —0°087 
0°55 — 4°664 — 1290 — 0469 —0'176 
0°60 — 6°877 — 2°052 —o'821 — 0°345 
0°65 — 10°382 — 3°299 — 1°436 — 0666 
0"70 — 16°165 —5°471 —2°558 — 1284 


The smooth trend of these coefficients makes sufficiently accurate interpolation 
for other intermediate values of a/a’ easily possible, for it is mainly the general 
form of the curve represented by Q(t) as a/a’ changes that matters, and our 
aim is to find that value of a/a’ giving the closest fit to the (adjusted) observations. 
Another feature of this table is that the coefficients corresponding to sin 5D 
and sin6D can also be roughly estimated for any given a/a’ from inspection of 
the general run of the earlier coefficients. 


7. The solution giving the best fit.—The equations of condition (14) have 
been set up and solved for a series of seven different values of a/a’ ranging from 
0:50 (Bode’s law) to 0-65. ‘The mean square value of the resulting sets of residuals 
for each solution is shown plotted against a’/a in Fig. 2, from which it is seen 
that there is a pronounced minimum at just about a’/a= 1-6, and there is obviously 


* Yale Trans., V1 and VII, Cambridge, 1932. 
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little point in endeavouring to improve on the round figure of 1-6 for the ratio of 
axes. ‘This corresponds to a mean distance for the unknown planet of 30°71 a.u. 
compared with the true value of 30-07. 

The longitude thereby predicted for “ 
the date of discovery is readily found to 
be 329°-4 compared with the actual value 
of 328°-4. ‘The error of the present method 
is therefore about 1° ahead of the true 











m.s.d. (secs. of arc) 





























position, whereas that of Le Verrier was \ See 
° ° . = _ ae 
just under 1° behind. : a 
Consistently with the closeness of a’ to \. 
the true value, the present method also bts — hy — 1 


gives a closer value of the mass of the Ia 

unknown body than was obtained by either Fic. 2.—Curve of mean square residuals 
Se Venter as Ades ‘Tho wits Goad computed Sor a number of values of a’la. 
: & or: are UNG The best fit is seen to occur at about a’la=16. 
is m’/©=1/25 000 in round figures. 


8. Prediction at a date other than 1846.—The orbital elements (and mass) 
obtained by the various methods, adopting the better of the two solutions given 
by Adams, are as follows. 


Yasce III 


Adams Le Verrier Present method Neptune 
a 37°25 36°15 30°71 30°07 
e 0°1206 0'1076 ° 00086 
wo 299° 285° —_ 44° 
O/n’ 6670 9350 25500 19300 
Longitude at 330°°9 327°°4 329°°4 328°°4 


discovery (equinox 1950) 


It is a simple matter to find the longitude that these solutions would predict 
at any other time than 1846-73, the time of discovery. ‘The three curves of 
Fig. 3 show the difference between the longitude calculated in turn from the 
above values and the actual longitude of Neptune. It is evident from the curves 
that some element of good fortune entered into the predictions to make them 
as good as they were. Le Verrier’s solution appears superior to that of Adams, 
and his prediction would have been better still at a somewhat earlier date 
(supposing the same material to have been available, which could not in fact 
have been so) since the error changes sign at about 1842. Both solutions place 
the unknown planet ahead of the true position in the earlier part of the range, 
so that with a’/a assumed nearly equal to 2, the hypothetical planet will move 
more slowly than the actual one and therefore must eventually coincide with it 
in longitude. For Adams’ solution this does not occur till as late as 1856. But 
their orbits were intended to fit the observations prior to 1840 in Adams’ case, 
and 1845 in Le Verrier’s, so comparison can only properly be made for times 
earlier than discovery. ‘This shows that the errors implied by Adams’ solution 
are on the whole roughly twice as great as those of Le Verrier’s. (A possible 
reason for this is discussed below.) On the other hand the circular solution 
arrived at by the present method shows an accuracy comparable with that of 
Le Verrier’s orbit, but persisting over a much longer range in time. ‘The three 
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curves show that +1° accuracy would be obtained by Adams’ solution only 
for about 5 years, by Le Verrier’s solution for about 12 years, while for the present 
solution it would obtain for about 24 years. 

The orbits arrived at by Adams and Le Verrier both exhibit the curious 
feature that their difference in longitude from the true position is at its least 
fairly near the time of prediction—six years earlier for Le Verrier and nine years 
later for Adams. Why this place of “‘ best fit’’ should have come near the crucial 
time seems something of a puzzle when the observations utilized extended back 
some 150 years. It may well be that the latest observations were of particularly 
strong influence, but in any event it would seem that some element of good 
fortune attended the degree of closeness of the actual prediction. 
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Fic. 3.— Showing the amounts by which the orbits of Adams, Le Verrier, and the circular 
orbit solution differ from the true position of Neptune. 
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g. Other points.—The weaker agreement of Adams’ orbit with the actual one 
seems very probably to arise from what appears to be an invalid step made by 
him at the very outset of his work. Adams says (loc. cit., I, p. 11) ‘‘It is easily 
seen that the series expressing the corrections of mean longitude in terms of 
the corrections applied to the elements of the orbit, is more convergent than 
that which gives the correction of the true longitude, and the same thing is true 
for the perturbations of the mean longitude, as compared with those of the true. 
The corrections found above were accordingly converted into corrections of 
mean longitude by multiplying each of them by the factor r?/ab.’” Now whereas 
such a relation as that here claimed subsists between the true longitude and 
the mean longitude in an undisturbed orbit in virtue of the equation 


dv 
2 = = h 
“-" nab, 
which may be written 
(r2/ab) dv =n dt =dl, (15) 


there would seem no validity in supposing that in perturbed motion the 
variation in true longitude Av and that in mean longitude A/ are similarly related. 
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No such relation seems known in celestial mechanics. We have in any case that 
v= | arvay, l=nt +«, 
and so for a variation A, in which all the elements are changed, 
Av= | (r-*Ah — 2hr-3Ar) dt, Al=tAn+ Ae, 


and that the ratio of these cannot reduce to a form independent of the variations 
of the elements, such as (15), seems evident. But the actual case is more complex 
still because the discrepancies involve also the corrections to the orbit of Uranus. 
The validity of this step has also been queried by E. W. Brown*. 

In his paper E. W. Brown also was concerned with the possibility of finding 
some simple criterion for the prediction of an unknown planet, and the method 
he proposed therein rests on consideration of the series —m’' > F;siniD for 
the principal perturbations in longitude, naturally enough, as here. Brown’s 
discussion shows how the general trend of this function (with changing a/a’), 
rather than its precise value at all stages, is sufficient to estimate the position 
of the unknown planet, though he applies the method on the basis that the first 
three terms of this function are adequate. Independent preliminary determin- 
ation of the instant of conjunction is not included in his method, and one gets the 
impression that this could be a considerable defect when it comes to application. 
There is difficulty in being certain, however, from the actual work given in 
his paper, since the present author has not found it possible to recover the 
numerical values he gives for several of the various coefficients. At all events, 
Brown himself reaches the conclusion that his method predicts conjunction as 
occurring close to 1840, which in fact is almost go” from the true position, and 
would lead to a longitude at the time of discovery some 27° ahead of the actual 
place of Neptune, even assuming a circular orbit at exactly the right distance. 
Nevertheless, the general theory behind Brown’s method would appear to be 
sound, though, handicapped as it would be in detailed application by lack of 
knowledge of conjunction, to reach a prediction of accuracy comparable with 
those of Adams and Le Verrier might very well require considerable numerical 
labour. 

St. John’s College, 

Cambridge : 

1958 September 19. 

* M.N., 92, 93, 1931. See Section 14. 





STELLAR GROUPS. III. THE SPACE MOTIONS 
OF SOME MAIN SEQUENCE DWARFS 


Olin #. Eggen 


(Communicated by the Astronomer Royal) 


(Received 1958 September 8) 


Summary 


Using photoelectrically determined magnitudes and colours, photo- 
metric parallaxes have been determined for 359 stars with known radial 
velocity and proper motion that have been spectroscopically classified as 
main sequence dwarfs on the Mt Wilson or Yerkes Observatory systems. 
These parallaxes have been used to compute the U, V, W-vector space 
motions of the stars. ‘The stars were subdivided according to the size of the 
W-motion (perpendicular to the galactic plane) and the motions plotted in the 
U, V-plane in Figs. 2 to 5. These figures suggest that the distribution of 
space motions, at least for small values of W, is not random. 





Introduction.—In previous papers of this series (Eggen 1958a, b; Papers I 
and II) several widely scattered stars were assigned to five groups—the Hyades, 
Sirius, ¢ Herculis, « Indi and 61 Cygni groups—on the basis of the apices and 
total amounts of their apparent motions. ‘This method—which we will call the 
‘““convergent-point method ’’—of choosing group members has the advantage 
that the parallaxes, which are accurately known for very few stars, are not used. 
After the members of the above-mentioned groups were chosen by the convergent- 
point method, the distances computed for the individual stars, from the group 
motion, were shown to agree satisfactorily with those derived from the observed 
parallaxes, within the very wide limits of uncertainty of most of the latter. Also, 
several tests, discussed in Papers I and II, indicated that far too many stars 
satisfied the requirements for group membership to be the result of chance 
coincidences of apparent motion. Nevertheless, the reality of the groups will 
remain in doubt until the space motions of the individual stars can be examined 
and the design of the present paper is to explore the possibility of using photometric 
parallaxes for this purpose. 

Photometric parallaxes.—In a previous paper (Eggen 1956) photometric 
parallaxes were computed for 285 main sequence stars with annual proper motion 
exceeding 0”:5. In addition to the correctness of the fundamental assumption 
that main sequence stars later than, say, Go have a unique colour—luminosity 
relation, the accuracy of such parallaxes depends upon (1) the accuracy of the 
colours and magnitudes, (2) the freedom of the observed colour from contamination 
by interstellar reddening or the presence of a companion star, and (3) the validity 
of the assumption that a star is a main sequence dwarf. There is some evidence, 
mainly from open clusters (Eggen 1955), that (a) spectral classification on either 
the Mt Wilson or Yerkes system—at least in the range from Go to K5—can 
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distinguish the main sequence dwarfs from the subdwarfs and subgiants; and 
(6) the colour—luminosity relation for these main sequence dwarfs, in the same 
spectral range, is given by M, = +2™-56 + 5-84(P— V),,, for all values of (P—V), 
between +0™-5 and +1™-24. ‘Therefore, using all of the magnitudes and 
colours determined on the (P, V),,-system, or reducible to that system, photo- 
metric parallaxes have been determined for the 359 stars in Table I. ‘The 
parallaxes for stars with annual proper motion greater than 0”:5 have been taken 
from the previous publication (Eggen 1956) where a comparison with the 
trigonometric values is given; a few values have been altered on the basis of 
more recent magnitudes and colours. 

Space motions.—T able I contains ail stars for which (1) spectral classification of 
dG1 to dKs inclusive, or G1 V to Ks5 V inclusive, has been assigned by Joy (cf. 
Wilson 1953) or by the Yerkes observers, respectively; (2) radial velocities are 
known; and (3) photoelectric magnitudes and colours, + 0™-5 <(P—V),,>1™-24, 
are available. Also included are a few objects which, although just outside the 
spectral range given above, have colours in the accepted range. ‘The radial 
velocities were taken from Wilson’s catalogue (1953) except for a few published 
since that compilation was made. Many of the Yerkes MK classifications are 
by Roman (1955) but some additional classifications by Morgan and others have 
also been taken from the literature. In the main the Yerkes and Mt Wilson 
classifications agree for stars classified at both observatories, but two or three 
subdwarfs, on the Mt Wilson system, have been classified as main sequence 
objects by Miss Roman, who has also classified as subgiants or giants a half-dozen 
stars later than G5 which were called dwarfs at Mt Wilson; in all cases of doubt 
about the luminosity class, the star has been omitted from Table I. Visual 
binaries have also been excluded except when the difference in visual magnitude 
is known to be less than one or two tenths, in which case ‘‘ AB’’ follows the 
HD number in Table I and the parallax was computed for the mean component ; 
double-lined spectroscopic binaries have been excluded. The proper motions 
were, in decreasing order of preference, taken from the GC, Yale zone catalogues, 
Cincinnatti catalogues, and some unpublished Greenwich results. Of the stars 
which otherwise meet the requirements stated above, colours and magnitudes 
are available for about go per cent of those with annual proper motions in excess 
of 0”+5, 50 per cent of those with smaller proper motions but with radial velocity 
exceeding 50 km/sec, and only about 25 per cent of the remaining stars. 

Because of the manner in which the stars in Table I have been selected, the 
parallaxes undoubtedly represent, within the framework of the assumptions, 
the most consistent set available at present. ‘The uncertainties in the proper 
motions and in the radial velocities, which are all of quality a, b, or c (Wilson 
1953), will cause some blurring of the space motions as will uncertainties in the 
colour, where the probable error of o™-o1 leads to an error of 2 per cent in the 
parallax, but gross inaccuracies in the space motions will only arise from an error 
in the assumption that the star is a main sequence dwarf. The space motions, 
which are all referred to the Sun, are given in the vectors: U (directed away 
from the galactic centre; /=148°, b=0"), V (in the direction of galactic 
rotation; /=58°, b=o°), and W (perpendicular to the galactic plane; b=go’). 

All of the stars in Table I are plotted in the U, V-plane in Fig. 1; the filled 
circles represent stars with positive values of W and the cross at U= +11 and 
V = —17 km/sec represents the velocity of the local standard of rest with respect 
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TABLE I 
Space motions of 359 main sequence dwarfs 
HD/BD .Z*  u _’ wi wpe 7, vu _v oo wp 
(0” oor) km/sec (0”-001) km/sec 
166 75 + 16 — 20 — 10 — 7°699 73 + 64 + 23 — 16 
245 21 + 29 — 96 — 45 24916 71 - 5 - 2 — 16 
834 ~—s5I ‘a = 3 - 2532934 + 76 -297 + 35 
3443 64 + 89 — 49 — 18 25680 60 +29 —-12-—- 7 
3628 33 +108 — 39 + 43 +32°719 22 +163 —-193 — 54 
3651 I1l — 36 -— 20 + 11 26965 193 — 96 — 12 — 46 
3765 65 — 21 — 78 — 22 27274 82 +32 - 4 + 4% + 
3795 63 + 26 — 50 + 51 =+24°659 26 +79 — 28 + 37 
4628 = 124 o — §3 — 12 28343 62 — 34 + 12 + 15 
4747 61 ‘+ - - 5 + 38 +H -~-@ +3 + 
5133 69 + 37 -— 20 — 6 29400 28 — 33 — 62 + 43 
5351 36 ‘+e ~&% - @ 29587 29 “Hy - gs +s 
+ 63°137 67 + 82 — 63 + 29 30003AB 38 +22 — 14 — 6 
8262 38 ~s ~ + 38 30455 37 ‘— = -s 
9407 51 —- 46 — 4 ° 30495 87 + 20 — 6 ° 
9540 68 + 9 —- 22 + 3 31501 29 +590 —- 71 + 52 
9770AB 50 +20 —- 2 — jo 32147 = 125 + § = 563 = ¥3 
10126 Ss ae: A a a ei 32923 74 +a - 2b +2 
10145 a4 +95 — 52 — 11 +55°960 ra — 49 —133 —II0 
10307 94°5 + 33 — 25 + 1 +29°847 18 — 65 — 15 — 6 
10436 = 73 ~“@a-R = —% 34673 68 +h - 8 +S 
10476 142 — 35 — 24 + 4 36003 82 —- 74+ 2 — 18 
10780 106 +22 — 16 — 5 36283 20 +29 -— 77 — 62 
13043 938 -_ - + gy 36435 54 “ee - ¢ > & 
13435 92 +13 — 12 + 36 36443 22 « =— 55 — 94 — 12 
13445 85 +106 — 79 — 20 37008 44 — 46 + 25 -— 52 
13783 22°§ — 53 + 14 —I05 37394 go +12 —-21 - 13 
14214 60 +48 + 9 + 4 37706 46 —- 43 —- 4 -— 23 
14412 64 + 13 + 34 -— 10 38230 53 —- 25 — 65 + 10 
+70°169 27 + 62 — 8 + 8 38392 »=—s«115 —- 18 + § -II 
+ 4415 18 —165 -— 78 — 66 39881 50 —- 15 — 55 + 10 
16160 147 +71 + 1 + 33 +26°1067 22° -— 90 — 43 — 75 
16591 69°95 + 19 — 29 — 37 42618 41 — 644 — 14 + 8 
16619 20°95 + 19 — 36 — 40 42807 53 - 6 — 26 —- 5 
+ 48°739 2 ~~ -ag - 43147 “4 ~ | -  - SS +. 
17155 39 64 ~ oe 12 43834 112 ~~“ ~~ eB +. 
17925 105 + 16 -— 19 — 10 45289 48 + 86 -— 22 — 4 cr 
18702 38 +184 — 64 17 47752 58 — 51 — 19 —- 2 
18757 45 + 64. -— 79 -— 25 +47°1355 57 + 50 — 44 — 26 
19305 68 —- 31 -— 60 — 1 50281 118 —- 1 + 14 - 21 
20165 47 —- 1 — 58 + g +40°1758 46 +49 — 44 +11 
20630 ~=Cs III +22 —- 4 - § 50806 71 + 9 —- 79 — I0 
20766 68 + 86 — 54 + 22 51219 37°55 — 43 — 48 — 34 
20794 170 + 75 -— 93 — 26 51608 31 — 20 — 35 — 2 
21197 80 * *+@ = 55458 33 ~~ +7 -B 
22049 321 + 4+ 7-19 56274 27 +92 + 3 — 61 
22496 78 + 36 — 24 — 13 +33°15S05 44 -7 -# + 3 
23439AB 31:5 +112 —155 — 93 +46°1264 22 + 47 — 50 — 66 
23484 70 + 33 — 16 — 16 +50°1435 14 — 74 — 18 — 20 
24002 30 —- 43 -—-100 — 11 60298 36 —130 — 20 — 47 +; 
24331 35 +90 + 8 — 10 — 1°1792 21 - 4-77 +11 
24451 66 +47 -—- 14 — 2 61606 66-5 - 28 —- 2 —- 9g 
24616 65 + 50 — 69 — 62 +54°1175 45 + 28 -— 50 — 16 
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/ 7, 
HD/BD se” eae) 
62613 59 
63685 40 
64606 39 
65430 = 42 
6558344 
66171 21 
68017 43 
68788 33 
+31°1781 51°5 
69830 83°5 
70352 _-28°5 
+ 33°1694 26 
71881 31 
72614 33 
7267378 
73393 _26°5 
73667 46 
74377 38 
74576 91 
74842 45 
75530 ~— «18 
75632AB 83 
75732 105 
78558 32 
79096 59°5 
79969AB 56 
84035 70 
84737 86 
85512 -94°5 
86728 87 
87884 34 
88371 20 
233719 33°5 
89668 36 
89777 16 
go711 38 
+56°1458 81 
+ 46°1635 3a 
+42°2163 17 
95128 go 
97584 8971 
98281 41 
100623 85 
10134925 
IOI50I 102 
101581 73 
102158 21 
102345 105 
104304 98 
104556 77 
104988 30 
+22°2442 13 
106116 39 
106156 36 
106210 33 
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TABLE I (cont.) 


V 


km/sec wd 
°F 7g 
= 3I ° 
—s7 = F 
= So = 3S 
—Iai = 41 
=- 3 =- 4 
— 66 — 44 
=— 52 = 23 
= 97o — 30 
—- 58 -— 10 
= 5 + 22 
—106 — 66 
= 4 wr 4 
=@G + @ 
- § = 2 
—- 50 — 6 
— 4 — 34 
— 78 — 48 
= a ° 
— 10 — 20 
—120 — 0% 
+ 8 — 19 
—m = 1 
= = 
= £ + © 
= 7 = FS 
~ ie = § 
= £ + 2 
a = a 
= 34 + 23 
= & = F9 
— 2 + 21 
—103 + 16 
= 5 — 3° 
=— 4 + 10 
— =~ 1% 
—- 2s + F 
-2a + 9 
—143, — 90 
—- 1 +4 
oe = 
+ tt + 3% 
+ 2 + a7 
— i = 10 
—- 22 + 8 
- ia + 7 
— 92 + 28 
— 38 + 6 
—- 12-13 
=o = ag 
~*~: = & 
— 50 — 26 
—- 4 + 26 
— 26 — 20 
— 60 — 51 


t 
HD/BD (0”-001) 
+29°2279 58 
107582 21 
107596 58 
108754 17 
+ 9 2639 58 
109200 55 
110463 44 
110833 77 
III515 26 
112758 44 
112943 36°5 
113101 19 
114060A 26 
114703 30 
115404 91 
115617 = 132 
116012 37 
116495AB 51°5 
117635 43 
119070 20 
+ 72692 28 
120690 54 
+13°2721 50 
122742 64 
124752 29°5 
125072 109 
125161B = 33°5 
125455 49 
+ 30 2512 74 
126053 61 
127506 45 
127871 42 
+34°2541 46 
129747 20 
+17°2785 53 
130669 25 
130871 32 
13099259 
131511 96 
131582 = 35 
132142 39°5 
134319 20 
134331 36°5 
135101B = 37°5 
— 3°3746 32 
135204AB 44°5 
136202B) ss 41°5 
136274 31 
13635273 
136834 49 
137826 20 
138549 29 
139323 60 
139777 5° 
140901 72 


7 


bt+tetttt 1 t+ettgei tei tteiges 


t++etei +¢el ti ¢ttesti 
COW OAL Dsl 
nuon vO O 


+++ 


$2 


106 


N 


> WwW 
Nb WwW 


wm NN 
- hu 


sIDN 
Cw 


- 
we 
Www 


563 
V 
km/sec 
- 39 - 3 
—113. — 40 
— 27 + 27 
—164 — 93 
—- 6 -— 3 
= 20 — 36 
a ae. 
= i? + 3g 
— 96 — 56 
—35 + % 
Teme Ser > Se 
- - & 
= © + 3 
=— Fe + 9 
~ fh = & 
7 ee =e 
a 
= 9 = 38 
— 17 — 16 
~@ - = 
— 66 — 15 
= = 9 
—- 68 - 4 
—- 13 — 22 
+ 9 -—- 1% 
—- 12 — 27 
—- 13 - 14 
—- 50 + 8 
— $5 = 9 
— 18 — 42 
—- 25 + 8 
— 4 = 18 
—- 51 — 18 
- 4- 6 
— 62 + 27 
-e -F 
— 66 + 13 
- 4 —- %9 
- 15 - 17 
— 81 + 28 
—- 59 + 26 
—- 16 + 4 
— 2, = ee 
— 36 + 16 
— 48 ° 
—-128 - | 
- 10 — 19 
—- 8: + 23 
— 43 + 36 
- - 
— 68 + 16 
— 28 + 20 
-G& —- 3% 
— 22 ° 
— 31 + 6 
39 





564 

HD/BD te” 201) 
143291 32 
143761 75 
144287 47°5 
144515 30 
144579 56 
144628 58 
144872 39°5 

+35°2774 30 
145417 43 
145675 85 
145958AB 39 
146233 79 
148530 22 
148653AB 46 
149661 105 
149957 41 
150248 37 
150437 30 
151188AB 36 
151288 112 
151337 61 
151877 30°5 
152391 60°5 
152792 43°5 
154345 52°5 
154363 86 
154577 56 
155456 36 

+42°2810 30 
155885AB 157 
156274 117 
156802 26°5 
157214 67 
157750 26 
157881 138 
158633 65°5 

+ 31°3025 15 
159222 48 
159868 4! 
160346 80 

+71°851 38 
160691 107 
161198 40°5 
161848 25 
162756 28 
164922 56°5 
165401 42 
166348 82 
166620 87 
169822 30 
170493 67 
170357 20 
170657 72 
171314 = 45 
171627 93 


P+t+lL tei ttl teetei tei +i 


l+++1 


Stellar groups: III 


V 
km/sec 
— 7 
— 26 
— 12 
— 72 
— 65 
— 12 
+ 44 
— 18 
— 153 
=e 
- 13 
— 85 
— 78 
— 12 
- 1 
—- 8 
— 28 
— 10 
— 16 
— 13 
- 39 


i 
ww <I 
wan - 


—129 
— 80 


— 18 


— 97 
— 104 
— 26 


TABLE I (cont.) 


W 


b+l+e+4+4+ 


Soe Th ih eg 


| 


b+++14 


L+t++i 


i ia 


b+++ 1441 


oe 


L+++1 


28 
19 
18 

9 
II 
24 

3 
28 
44 

a 
19 
21 
41 
14 


= 


43 
14 
13 
I 
15 
4 
4 
6 
7 
34 
19 
20 


HD/BD 


(0” oor) 
172310 23 
+31 3330 39 
174080 67 
229599 50°5 
175541 46 
230409 Ue 
178126 37 
178428 = 59 
179957AB 46°5 
+ 41°3306 25°5 
+33 3433 33 
182488 78 
183255 46 
231683 22 
184467 81 
184700 18 
185144 180 
186427 61 
186858AB 48 
225732 48 
189087 36°5 
189567 58 
190007 92 
190067 47 
190404 52 
190406 67 
191069 25 
191408 140 
191785 55 
192031 22 
192310 }=6128 
193202 66 
194640 8958 
195987 55 
196761 61 
196794 38 
197076 = 49°5 
197214 43 
197484 16 
+ 19-4499 40 
199476 = 33 
+12°44909 = 42 
200779 73 
201091AB 296 
202123 26 
202457 54°5 
204587 = 53°5 
204814 34 
205067 30 
205390 = 65 
205855 =. 29°5 
+ 24°4460 21 
+ 4°4762 13°5 
209100 8280 
209742 33°5 


T pt 


b+++ter tert 
N 


+l ttt 
a 





Vol. 118 


km/sec 


_ 19 


+172 


oe 1 eal Vo aaa I 


4 


bth i tei ti tee reer ei si 


28 


_ 








ei i. i ae oe oe Oe ee Ue: he ee ee, ee ee a, ee ee, i, 


wwe ww Tw ww PP VS FP WT VFN OOO NS VF WV we wee wry | y* 


_ 





TABLE I (cont.) 


























Tens V 7 V 
pt pt [ 
—ae (0”-001) - km/sec ad ri aed (0”-o01) - km/sec . 
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214749 79 +1 - 7 = 2 221914 32 + 97 -— 40 + 13 
214759 47 ‘a - 8 + 4 222237 89'5 —- 32 — 75 — 23 
216054 33 — 61 -— 16 + 10 222794 39 + 28 — 84 + 50 
216777 —s +e - |S - BS 222935 34 +193 - 8 = $5 
216803 =-130 + 2-7-1 223238 29 +45 -— 36 -— 2 
217014 85 + 11 ~— 28 + 16 223498 26 + 88 — 28 + 19 
+ 68°1345 24 +109 — 75 + 5 223778 + =112 +11 —- 6 — 1 
217357. +4110 —- 2 +15 — 4 224383 25°55 + 61 — 74 + 3 
219134 86175 +45 - 37 -I11 224618 21 +229 —136 — 6 
219175B 25 +o - 3 - 8 
' +229 ' ' ' ' ~ +169 t T t J T ' ' ' a 
1) 
[ ot ? | 
0 o 4 4 
uf 1 
L 
00 
F 1°) 
+ -279 
50} 
| O 
O}-4+_+_ + 
5 
- 50+ 
-100+ 
r 
r 
Vv F 7150 ; 





Fic. 1.—Values of U and V for the 359 stars in Table I. Stars with positive values of W are 
plotted as filled circles. Circles with 6 km|sec radius surround the values of (U, V) for the Hyades, 
¢ Herculis, « Indi, and 61 Cygni groups. The velocity vectors are all referred to the Sun; the 
cross at U= +11 and V=—17 km/sec represents the velocity of the standard of rest for the nearby 
stars. The arrows indicate stars falling outside the boundaries of the diagram. 
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to the Sun (Woolley 1958a). The values of (U, V, W) for the Hyades, ¢ Herculis, 
« Indi and 61 Cygni, taken from Papers I and II are the following: 


U 
+92 
+79 


+54 
+40 


V 


— 53 
—39 
—45 
—18 
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W 
— 6km/sec 
+ 6 
— 26 
-— 2 


These values of U and V are shown in Fig. 1 as the centres of circles with a radius 
of 6 km/sec. 

It is of interest to examine the distribution of points in Fig. 1 for the presence 
of regularities and therefore, in an effort to reduce any blurring effect that might 
arise from the inclusion of all values for W, the material in Fig. 1 has been 
subdivided as follows: 




















Fig.2: W= oto —20km/sec, 
Fig.3: W =+1 to +20, 
Fig. 4: W < —20, 
Fig.5: W > +20. 
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Fic. 2.—Values of U and V for the stars with W from 0 to —20 km/sec. The symbols are as in 
Fig. 1. A few stars fall outside the limits of the figure. 


A few stars falling outside the limits of each of these diagrams are not plotted. 
The values of |W|>50 km/sec in Figs. 4 and 5 are indicated by filled circles. 
Since the Sun’s motion in W, with respect to the local standard of rest, is 
—7 km/sec (Woolley 1958 a), the limits for W in Figs. 2 and 3 become +7 to 
—13 and +8to +28 km/sec, respectively, referred to that standard. ‘The values 
of (U,V) for the four groups of high velocity have been marked in Figs. 2 to 5 
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Fic. 3.—Values of U and V for the stars with W between +1 and +20 km|sec. 




















as in Fig. 1; a few stars fall outside the limits of the figure. 
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FiG. 4.—Values of U and V for stars with values of W< —20 km/sec. 
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indicate stars with W < —sokm|sec; the other symbols are asin Fig.1. A few stars fall outside the 


limits of the figure. 
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as in Fig. 1. For comparison with Fig. 2 the values of U and V derived from the 
trigonometric parallaxes of the stars in Gliese’s (1957) list of objects within 
20 parsecs of the Sun and with W=o to —20 km/sec have been plotted in 
Fig. 6; a few stars fall outside the limits of this diagram. 
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Fic. 5.—Values of U and V for stars with W > +20km|sec. The filled circles indicate stars 
with W > +50 km/sec; the other symbols are as in Fig. 1. A few stars fall outside the limits of the 
Giagram. 




















Fic. 6.—Values of U and V, for stars with W between o and —20 km|sec, computed from the 
trigonometrical parallaxes given by Gliese in his list of stars within 20 parsecs of the Sun. 
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In comparing Figs. 2 and 6 it must be remembered that, although the material 
used in both diagrams is nearly complete for stars within 20 parsecs, which are, 
in the main, low velocity stars, the larger number in Fig. 6 results from the use 
of all spectral types. On the other hand, the stars in Fig. 2 that are more distant 
than 20 parsecs are largely high velocity objects. The apparent regularity in 
the distribution of the low velocities suggested by Fig. 2 may indicate the blurring 
effect of inaccurate parallaxes used in constructing Fig. 6, but the reality of the 
two chains of stars, with V near — 5 and — 18 km/sec, separated by a lane, needs 
confirmation by a larger sample of stars; observations of magnitude and colour 
for this purpose are now in progress. In any case, the apparent connection of 
the Hyades with the chain at V near — 18 km/sec may be partly accidental since 
finer subdivisions in the values of M can break up the chain and remove the 
Hyades from the figure. This effect is better illustrated, perhaps, in the case of 
61 Cygni. The value of (U,V) for 61 Cygni, for which W= —6 km/sec, falls 
in a rather densely populated region of Fig. 2 while the values for ¢ Herculis 
and « Indi, both of which have W’s outside the range shown in the figure, fall 
in unpopulated regions. The ten stars, including 61 Cygni, which fall in the 
range, U= +80 to +100 km/sec and V = — 40 to — 60 km/sec are the following, 
arranged in order of increasing photometric parallax (in units of 0”-oo1): 


HD Tot U V W 
216777 24°5 +97 —55 ~wg 
219175B 25 +94 —55 —II 
10145 2 +95 —5§2 sy 
199476 33 +86 — 60 —- 7 
112943 36°5 + 83 —50 “3 
64606 39 +80 —57 — 7 
183255 46 +95 — 46 —12 
130992 59 +89 —49 —19 
3443 * (64 +89 —49 ~ 
Mean ae +90 —53 —12 
61 Cyg 296 +92 —53 — 6 


The mean values of U and V for the nine stars are nearly identical with those 
for 61 Cygni but the values of W differ by 6 km/sec. Since the apex of the apparent 
motion of 61 Cygni was used in selecting—by the convergent-point method— 
the group members in Paper II, and the value of (A, D) used in that method is 
sensitive to W, more group members perhaps would have been found if the 
mean value of W given above had been used. However, tolerances far wider 
than those used in Paper II would be necessary to include, by the convergent- 
point method, all of the above stars in the 61 Cygni group. The question arises 
as to the spread in U and V, as well as in W, that can be tolerated in considering 
stars as members of groups. ‘The answer, of course, depends upon both the 
definition of ‘‘group’’ and the relation of the groups to the general dynamics of 
the galaxy (Woolley 19584) and will not be discussed here. 

A suggestive clumping of stars in Fig. 2 occurs near (U, V)=(26, — 50). 
Some of these stars had previously been associated, by the convergent-point 
method (the # Centauri group), with a group moving toward A=6" 12™, 
D= —27°:2 with a velocity of 60 km/sec. The six stars that form this group 
in Fig. 2 are the following: 
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HD/BD Tot U V W 
90711 38 +27 —54 —19 
+54°1175 45 +28 —50 — 16 
159222 48 +30 —50 — 6 
136834 49 +23 —53 —18 
120690 54 +31 — 48 — 9 

+ 30°2512 74 +25 —5I —19 


The following additional stars, among those within 20 parsecs of the Sun, may 
also belong to the same group: 


Name Te U V W Ter 
6 Cen 60 +26 —48 -— 5 76Y(8), 39C(7) 
+ 37°748 64 +28 —47 —19  83M(r1o), 51V(16) 
GAZ 82°1111 110 +30 —48 —19 93M(12), 104 Yk(10), 106G(12) 
+ 18°638 112 +27 —54 —I!I 99M(t0), 104V(12) 
Wolf 359 = 400 +26 — 46 —19  404W(6), 397V(4) 


The values of 7, given above are those used in computing the space motions; 
the trigonometric parallaxes, with the observatory designations and weights, are 
taken from the Yale catalogue (Jenkins 1952). 

In comparing Figs. 2 and 6 it is interesting to note the near absence of stars 
in the former that might belong to the Sirius group (Paper I) for which 
(U, V, W)=(— 14,0, —12) while the latter shows a heavy concentration of 
group members. The scarcity of dwarf members of this group, with colours 
in excess of +0™+5, was noted in Paper I and seems to be confirmed by the 
present comparison since bluer stars are contained in the material used to 
construct Fig. 6 but not in that used for Fig. 2. 

Since the division point of W=o, used to separate the stars in Figs. 2 and 3 
is an arbitrary one, we might expect some of the features of Fig. 2 to appear also 
in Fig. 3. This seems to be the case, but there are some curious differences—the 
chain of stars in Fig. 2 with V near —5 km/sec and with small values of U does 
not appear in Fig. 3, indicating the relatively small range of W for these stars. 
Members of the Hyades group, for which W= —2, appear in both figures but 
they are accompanied in Fig. 3 by a chain with V near —25 km/sec. Members 
of the « Indi group, absent from Fig. 2, appear in Fig. 3 as expected since the 
value of W for this group is +6 km/sec. The distribution of the points in Figs. 4 
and 5 is more featureless than that of Figs. 2 and 3 but there are too few stars 
to permit any detailed discussion. Some of the regularities in the distribution 
of points in Figs. 2 and 3 have previously been suspected by Dyer (1956), in 
the distribution of the velocities of M-type dwarfs, and by Vyssotsky (1951), 
in the distribution of the velocities of the A-type stars and the K-type giants. 

Summary.—Photometric parallaxes have been derived for a selected list of 
359 main sequence dwarfs with colours, (P—V),, between +0™-5 and +1™-24. 
These parallaxes have been used to compute the space motions (U,V, W) of 
the stars. Although the stellar groups, whose members were chosen by the 
convergent-point method in Papers I and II, appear in the distribution of space 
motions of these stars, that distribution shows (a) the limitations of the 
convergent-point method of defining group membership and (b) the need for 
a dynamical definition of stellar groups. 
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The space motions in Table I were computed at least once by the author 
and once by either Misses L. Mather or A. Heather, using desk calculators, 
or Messrs A. E. Carter and D. A. Harragan, of the Nautical Almanac Office, 
using Hollerith punch card machines. I am greatly indebted to the Astronomer 
Royal for many stimulating discussions. 


Royal Greenwich Observatory, 
Herstmonceux Castle, 
Hailsham, Sussex: 
1958 September 1. 
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ENERGY LEVELS AND TRANSITION PROBABILITIES: OF Fe IV 


R. H. Garstang 
(Communicated by the Director of the University of London Observatory) 
(Received 1958 September 17) 


Summary 


Estimates have been made of the quantum mechanical parameters of the 
3d® configuration of Fe tv, and the energy levels of this configuration have 
been calculated. ‘Transition probabilities have been computed for magnetic 
dipole and electric quadrupole radiation for transitions between the levels of 
the 3d5 configuration. A search is made for [Fe Iv] lines in the spectrum of 
RR Telescopii, and a number of identifications are suggested. 





1. Introduction.—The absence of any term analysis of the Fetv spectrum 
prevents precise predictions of the wave-lengths of forbidden lines which might 
be of astrophysical interest. The positions of a few terms were estimated by 
Edlén (1) and, more recently, by Glad (6). Edlén pointed out that the *G-4*F 
multiplet should occur in the visible region of the spectrum, but that as it 
required a rather high excitation energy it would probably not be very strong. 
Swings and Struve (12) suggested that a line A5040-1 in RY Scuti should be 
attributed to [Fe1v]; they pointed out that this line also occurs in AX Persei 
(12), NGC6572 and NGC7027 (18) and in Nova RR Pictoris (1925) (13) 
and suggested that it might belong to the *G—‘F multiplet. Swings (14) 
suggested that some unidentified lines in RY Scuti might also be due to 
[Fe1v]. Thackeray (15) noted the presence of a number of unidentified lines 
in RR Telescopii, between A 4835 and 45288 and in the infra-red, some of 
which he suggested were due to [Fe1v]; excitation conditions in RR Telescopii 
were favourable to the production of this ion. In a later paper (16) Thackeray 
quoted Edlén’s identification of lines at AA 4868, 4900, 4906, 4918, 7189 and 7221 
as due to [Fe iv], but the transitions involved in these lines were not specified. 
The object of the present paper is to examine whether quantum mechanical 
calculations can contribute to the identification and analysis of the Feiv 
spectrum. We shall compute all the energy levels of the lowest configuration 
(3d°) of this ion, using as a starting point the predictions of a few terms by Glad 
(6), and then obtain the transition probabilities of forbidden lines between 
these levels. Having found the lines which should be strong we shall attempt 
to identify them in the spectrum of RR Telescopii. ‘This paper has the further 
object of providing transition probability data for a type of configuration (d°) 
which has not been studied previously in this connection; it also fills an 
important gap in the data for ionized iron, and forms the fifth and penultimate 
paper in the writer’s programme on transition probabilities of forbidden lines 
of elements in the first long period (2, 3, 4, 5). 


2. Method of calculation The theoretical formulae for the electrostatic 
energies of levels in the 3d° configuration have been given by Racah (g) and 
by Ishidzu and Obi (8). Corrections of the form «L(L +1) and BO were added 
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to the theoretical formulae along the lines of the work of ‘Trees (e.g. (17)) and 
Racah (10). The matrices of spin-orbit interaction for the d° configuration 
have been computed by Greyber (7); a copy of his results was sent to the 
writer by R. E. Trees, together with four corrections to Greyber’s work found 
by Trees. 

The centres of gravity of five terms ®S, 4G, *P, *D and 4F, were predicted 
by Glad (6). Our theoretical formulae contain five parameters, and these five 
parameters were chosen to fit these five terms. This is a somewhat dangerous 
procedure, in that it forces a fit which may be unjustified, and no check on this 
is possible if no other terms are known. No better procedure suggests itself 
at the present time. The spin-orbit parameter (, was estimated by linearly 
interpolating ¢? between Feri and Fev (3). The adopted parameters are listed 
in Table I. Using these parameters the complete energy matrices were set up. 
Their eigenvalues and eigenvectors were calculated. The smaller matrices 
were treated by hand computation; the three largest matrices were handled 
on an electronic computer. The line strength matrices were obtained in the 
usual way from Shortley’s formulae for magnetic dipole radiation (11) and the 
present writer’s results for electric quadrupole radiation (2). ‘The eigenvectors 
of the energy matrices form the transformation matrices from LS-coupling to 
intermediate coupling; this transformation was carried out for the two types 
of radiation. The final stage in the calculations was the conversion of the line 
strengths, obtained in atomic units, into transition probabilities. For this work 
the electric quadrupole radial integral 


tai | ” 2 P2(3d) dr 
0 
was calculated by D. F. Mayers from his unpublished self-consistent field, 
with exchange, wave function. He obtained s,= 1-14, and this value was used 
in the present paper. ; 


3. Results—The final energies are given in Table II. It is difficult to 
estimate the accuracy of these energy levels. Previous experience suggests 
that many of the intervals between the levels of one term should be accurate 
to about 50cm~', but that some will be substantially worse; most of the 
intervals between the terms should be accurate to about 500 cm~', but a few 
will almost certainly be less accurate. ‘This accuracy is rather low for the purpose 
of identifying predicted wave-lengths of spectral lines with observed transitions. 
In particular, the lines of a multiplet which extends over a very small range of 
wave-length may not be predicted in the correct relative positions in the spectrum. 
This is largely due to the vanishing of diagonal spin-orbit interactions in the 
half-closed configuration 3d°._ The observed term-intervals are second-order 
perturbations. ‘The spin-orbit interaction still causes first-order perturbations 
iff the wave functions, so that we should get reasonable estimates of the transition 
probabilities of forbidden lines between levels of different terms. Fortunately, 
three of the terms which are of most interest (2D, ?F and ?H) are predicted to 
have large term-intervals, and it is highly probable that, for example, lines from 
the level *H;, to the 4G term will be at shorter wave-lengths than *H,,-*G 
transitions. This is of importance when searching for identifications in celestial 
spectra. 
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TABLE I 


Adopted parameters (in cm~) 


A 37656 
B 1075 
Cc 3884 
a 104 
B 182 
Ca 468 
TABLE II 
Calculated energy levels 
Term J Energy Term J Energy 
6S 24 ° *H 43 55959 
2 2 6 
53 56435 *7 
4G 23 32213 - 
33 32243 van 2G 33 56872 wis 
43 32263 _ “4 43 57283 
53 32254 
- 2} 60316 28 
“p 3 35259 _ ga 34 60344 
ui $5873 _ 4 
™ atte 25 } 65501 
‘4p 4 38389 
13 38453 * 3D 1} 73182 53 
2) 38468  _'5 2$ 73235 
7 — 163 
32 38305 
2G 33 82193 _. 
2] 5} 47218 8 4} 82171 
6} 47296? 
i 4 98588 
i 1 os c — 49 
2 14 49589 8 I 99539 
24 4g10g 
2D 1} 107018 
2 1 24 106 — 39 
2F 2} 51864 _ 646 2 979 
33 51238 
4p 1} 52595 _ 3 
2} §2552 coi 
33 $2449 _ 
43 52414 F 


The transition probabilities of selected multiplets are listed in Table III. 
We have included all transitions from and between levels up to *D; most 
multiplets from the next five terms are given, only a few very weak multiplets 
being omitted. From the next five terms only the strongest multiplets are 
given, and from the highest two terms all multiplets have been omitted. In all, 
75 multiplets have been omitted; this omission seems justified for reasons of 
economy, having regard to the unlikely need for data concerning these weak 
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transitions from high-energy levels. It is, as always, difficult to assess the accuracy 
of these transition probabilities. As mentioned above they are mostly due to 
first-order perturbations, and should be moderately reliable. ‘There is, however, 
one cause of uncertainty which does not usually apply but which is of importance 
here. The unknown errors in the energies may introduce errors in the transition 
probabilities in two ways. First, the assumed positions of the terms enter into 
the calculation of the transformation matrices, and thence into the line strengths. 
Secondly, the conversion from line strengths to transition probabilities involves 
the third and fifth powers of the frequencies of the transitions. ‘The first source 
of error will chiefly affect transitions between different terms. ‘There is little 
that can be done to eliminate this source of error until an observational analysis 
of the Feiv spectrum becomes available. ‘The second source of error affects 
transitions within a term and between close terms, but has a relatively small 
effect on transitions in the visible spectrum between well-separated terms. 
There are a number of transitions for which the transition probabilities are small 
because of the frequency factors (v* or v°) but the line strengths from which they 
were computed are large. In order that these line strengths may be obtained 
in a few cases which may be of interest, the corresponding transition probabilities 
have been retained in Table III. 

Among the results it is worth noticing the high degree of metastability 
possessed by the 4G levels. The line strengths of the ®S—‘G transitions are so 
small that they could not be computed without a special investigation; the 
strengths came out as zero to the last figure retained in the calculations. An 
attempt was made to carry one section of the calculations to higher numerical 
accuracy and this led to a value of about 10~° for the probability of the magnetic 
dipole transition °S,,-*G,,. A similar attempt for the °S,,-‘G,, transition led 
to an upper limit of 10~’. Although these values are rather uncertain—one 
cannot say without a great deal of extra work whether they are sensitive to changes 
of the parameters in the calculation—it does seem that the 4G levels of the 3d° 
configuration are highly metastable. 

The weakness of the *P—*F multiplet is not entirely unexpected: it has zero 
electric quadrupole line strength in LS-coupling in the 3d° configuration. But 
one might have expected rather larger magnetic dipole strengths. There is 
probably little prospect of finding it in celestial spectra. 

There are no strong sextet-doublet transitions and, as Edlén pointed out, 
no transitions at all to the sextet ground state which are in the accessible region 
of the spectrum. 


le 
2 
‘ 


4. Tentative identification of [Fe1v] lines.—One of the striking results which 
emerges on examination of Table III is the intensity of some of the magnetic 
dipole transitions. Indeed, almost all the strong transitions requiring less than 
8 volts energy are magnetic dipole transitions, with the exception of the *G—'F 
transition, and even in this multiplet magnetic dipole transitions make substantial 
contributions to the total transition probabilities. Such searches as have been 
made for [Fe 1v] lines in celestial objects have of necessity been confined to those 
multiplets whose positions were roughly predicted by Edlén (1), together with 
an examination of the behaviour of unidentified lines with varying excitation 
conditions (for example in novae). Now that we can pick out from Table III 
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Transition probabilities for selected multiplets of [Fe IV] 


Transition 
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All entries in this Table are probabilities of spontaneous emission. 
The letters ‘ v.s.’ denote a line with a very small strength. See text for details 
of multiplets omitted from this Table. 


A 


7°6 x 1075 
7°6 x 1074 
6:3 x 1074 


0058 
0°034 
1°4 x 107% 
0°039 
0018 
3°5 X 107% 
07022 


0038 


A, 
<107° 
<1ze~* 
<16°* 


75x 107" 
9°9 x 10719 
3°8 x 10-7 
16x 10-** 
6°5:x 10°™* 


3°9 x 1075 
I°5 X 107-5 
v.s. 


iox20°° 
2°3 x 107° 
6:4 x 1071 
8-5 x 107'° 
13% 10°° 
3°4 x 10710 


4°9 x 10718 
1'2 x 190738 
7°4x 10715 


1°8 x 107 
6:2 x 107* 
1'o X 10° 
I°r X 107% 


I'°oX 1077 
6:6 x 1077 
6:0 x 1077 
7°2X 1077 
ir x 107° 
6°7 x 1077 
$2 x 107° 
6°8 x 1077 
I°5 X 107° 
1°7 xX 1077 


6:9 x 1078 
9°9 X 1077 
1°8 x 1078 
3°3 x 10° 
3°4x 10% 
3°1 X 107° 
1°5 X 107? 
55 x10°° 
2°9x 10 °° 
4°9 X10 





Units: reciprocal 
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| TABLE III (continued) 
Transition A, A, 
‘DD = 3 $-1 os 3°r X 107" 
cal 33-24 1°3 x 1074 8-7 x 10-4 
ails. ty} I'o xX 107° 12x 70°** 
4-2} = 1*9 x 107!9 
13-2} 1°3 x 1077 1°99 x 107!® 
IG] 33-54 = 1°0 X 1078 
42-53 370 x 10-3 3°3 x 1077 
: 5h-st 51x 1073 3°2 x 1077 
43-6} - 1-6 x 10% 
53-63 1°6 x 10-4 7" x 107° 
*J-*] 53-6} 5°9 X 1078 rg X 107" 
‘GD ~—s 2} -14 O°010 3°7 x 1073 
33-13 ar 1°3 Xx 107% 
23-2! O13 3°3 x 1074 
34-23 o'12 5°9 x 1074 
- os aer* 
J 43-2} 2°1 X10 
4sp2D }-1} 0°065 7:0 x 10-* 
13-1} 0°25 2:0 x 1075 
2$-13 0:069 7°3 x 10-* 
}-2} ee rx to"? 
13-2} O°10 5°9x 107° 
24-2} 0°54 2°8 x 10-* 
3D-D = = 23-14 1°5 x 107° 6-4 x 1071! 
‘GF —s 24-23 or21 4°5 x 107* 
34-24 0°47 8-1 x 1073 
44-2} we 4°1 x 107% 
23-3 O°o1s 2-7 x 10-§ 
34-33 0080 2°3 x 10-4 
43-3} 0°50 4°6 x 10-4 
5-3} 14 x 10~* 
‘p—"F 4-2} - 24x10 
14-2} 0'013 rs xX 10-* 
2$-2} 0°087 7-1 X 107 
1}-3} _ 2:0 x 10°5 
24-3} 5 x 10-4 25x 1077 
‘D-3F 4-2} wee 54x 10-4 
14-2} 0'021 3°0 X 10-4 P 
2}-2} 0°022 771x104 
33-2} 72x 10-4 I'9 X 10-4 
13-3} aes 2°3 x 1075 
2}-33 0034 2°8 x 10% 
33-32 o'10 3°3 x 10-§ 
;D-F 14-2} 0°035 1°8 x 1077 
24-2} 0°090 92x 1078 
13-33 — 5°9 x 107° 
23-3} 0°023 7x 107° 
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578 
‘Transition 
FF 34-2} 
‘G-‘F = 24-14 
33-13 
23-2} 
33-23 
4}-2} 
24-34 
34-3} 
44-3} 
54-34 
24-4} 
33-43 
43-4} 
54-4} 
‘p—*F oe 
13-1} 
23-1} 
4-2) 
13-2} 
24-2! 
14-3! 
23-3 
24-4} 
‘D—‘F j-1} 
ij-1} 
24-1} 
3h-1h 
}-24 
14-2} 
24-2) 
34-24 
14-3) 
24-3} 
33-34 
2}-4} 
34-44 
2D-"F sy 4-1} 
: ahr} 
1-2} 
2i-2) 
13-3} 
23-33 
23-4} 
‘F-"F = 43-3} 
4)-2} 
33-24 
33-1) 
24-1} 





TABLE III (continued) 


A m 


3°0 x 1073 
O14 


or19 
9°4x 1074 


O°014 
0°13 
0°039 


2°3.x 107° 
0°073 
oll 


1°5 x 1074 
I°5 X 1074 
v.s. 


2°4 x 1073 
5°3 x 1073 


I'l X 1079 


o'10 
016 
0039 


0°023 
018 
0'028 
9°4 x 107% 
o'10 

O14 
O'015 
1°8x 1074 
6:9 x 1073 
0'022 


o’o2!I 


1°4 x 107% 


5°8 x 1075 


4°6 x 10-* 


Ay 
10 X 107!9 


018 
0-071 
0°064 
Ol 
0-061 

4°11 x 107% 
0067 
O'ls 
0°040 

8-7 x 1075 
3°8 x 10° 
0045 
o'21 


41x 10-4 
1°6 x 1073 
5°6x 1074 
1°3 x 10-4 
4°6 x 10-# 
1°4 Xx 107% 
6-5 x 10-4 
9°9 x 10-4 
8-8 x 1074 


8-3 x 107° 
0016 
4°5 x 10-3 
1°9x 10-4 
97% 10°° 
41x10°% 
oO-012 
2°1 x 107% 
8-5 x 107% 
g°2 x 10-3 
ool! 
5°7 x 107% 
0°023 


9°5 x 10719 
3°5 x 1078 
3:1 x 1078 
4°0xX10-° 
1°6 x 10719 
1°8 x 1078 
7:0 x 107! 


9°8 x 10718 
3°5 x 10715 
1'2 x 10718 
4°8 x 10715 
8-8 x 10718 
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Transition 
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TABLE III (continued) 
Am 


0°13 
0042 
0°56 
0°61 
0°47 


1°5 X 107° 


0060 
o'18 
0°094 
2°4 x 1073 
6-1 xX 1073 
v.S. 
oO 4 
9°4 X 10 


V.S. 
V.S. 
0°077 


0°064 
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V.S. 
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GHG ah-4 
2p 1 
2F-2F —2}-2 

33-2 
23-3 
33- 
4F-2F 13-2! 
24-2 
33-2 
43-2 
13-3 
23-3 
33-3 
43-3 
3F-§F 23-3 
‘pS }- 
14- 
2}— 
Ip—s I 
1}- 
2} 
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14-1 
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1}-2 
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3)-2 
1 
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24-1 
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2}-2 
2FD = 25-1 
32-1 
23-2 
33-2 
MG ‘51-3 
53-4 
63-4 
sD-jG 13-3 
24-3 
24-4 
sG-3G 34-3 
4-3 
33-4 
43-4} 
2F2G =e 2}- 
3}-3 
2}-4 
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An 
6:0 x 1074 
7°9 X 10~* 
0°047 
O'051 
1°3x 107* 
0098 
0-016 
0'027 


0°026 
0'013 
Ol! 

28x 1077 
OrTE 
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0:078 
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0'027 
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a number of other multiplets which might be expected to be strong and for which 
the predicted wave-lengths are in the accessible region of the spectrum, a further 
search for [Fe 1v] lines is justified. 

Of the multiplets listed in Table III, 18 lie in the accessible spectrum. 
Six of these are rather weak, and are unlikely to be identifiable without precise 
wave-lengths. (These are *G—*I, *P—‘F, 4D-*H, *G-2G, *D-*S and 2?G-3G.) 
The strongest lines of the remaining 12 multiplets are listed in ‘Table IV, together 
with the wave-lengths predicted for their lines and the quantities (2/,,,,..+1)A 
which are proportional to the strengths of the lines. Lines with (2J +1)A <0-25 
have been omitted from Table IV; data for these lines can easily be obtained 
from Tables II and III. 


TABLE [IV 


Stronger [Fe IV] lines in the accessible spectrum 











Transition Predicted Possible identification in 
A Intensity RR Telescopii (1954) 
f (2J +1)A>0°25 A Int. Notes 
4‘G-2D 24-2} 5920 08 5911 3 I 
33-25 5931 0-7 5911 3 I 
‘P-2D t-y4 6978 03 
14-14 6936 Sze) 6997 4 2 
24-14 6899 0°3 
14-24 7178 06 7221 3 3 
2}-24 7138 3°2 7189 3 3 
‘G-3F 24-24 5089 1°3 5033 10 4 
33-24 + 5097 2°9 5041 15 4 
33-33 5265 o°6 
4-34 5270 4°0 5233 16 4 
‘P-3F 2}-2$ 5963 05 5911 3 I 
‘D-3F = 23-34 7831 03 
33-33 7732 08 
IG—'F 24-13 4905 13 4868 8 5 
33-13 4912 0°3 
24-24 4916 I'5 4868 8 5 
33-23 4923 o"7 
43-23 4928 o"4 
33-32 4948 1°6 4900 8 5 
43-3 4953 1'5 4900 8 5 
53-33 4951 0°3 
\ 44-4} 4961 12 4918 4 5,7 
53-43 4959 3°2 4906 12 5,7 
DF 3-1} 7039 "4 
14-14 7071 o'7 
23-24 7100 1'2 7221 3 3 
32-33 7070 o'9 
34-45 7088 1°6 7189 3 3 
40* 








XUM 


582 R. H. Garstang, Energy levels Vol. 118 
TABLE IV (continued) 
Transition Predicted Possible identification in 
A Intensity RR Telescopii (1954) 
(2J +1)A>0-25 A Int. Notes 
‘G“H 33-4} 4217 1°3 4209 7 6 
43-4} 4220 0-4 
54-44 4218 5°6 4206 10 6 
43-53 4137 7°3 4152 8 6,7 
53-53 4135 5°6 4152 8 6,7 
75 1}- 3 3297 1°5 
sD-3;D 13-2} 4239 °"4 
24-24 4144 0°6 
3F-§D 23-1} 7742 "4 
23-2} 7741 I'l 
33-23 7757 0-6 
sF-3G = 23-34 4571 1-2 
33-33 4577 1°8 
33-43 4575 1°6 


Notes to Table IV 

1. The three identifications of A 5911 are all plausible. 

2. Plausible. 'Thackeray’s suggestion of [Ti 111] is doubtful. 

3. A7221 and A 7189 might belong to either ‘P-2D or to ‘D-‘F. If the former is 
correct, the computed transition probability of 47189 may be too high. According to 
Thackeray, Edlén considers both lines as definitely [Fe tv]. 

4. Highly probable identification. 

5. The four lines AA 4868, 4900, 4906 and 4918, which Edlén ascribes to [Fe tv], almost 
certainly belong to the *G—*F multiplet. The identification of individual transitions is 
difficult. It seems likely that A 4906 is the 5}-4} transition, but the other suggested 
identifications are very uncertain. The small spread of wave-length in this multiplet 
hinders the identification of the lines. 

6. 44152 may represent either or both of the transitions 44-5} and 5}-5}. The 
identifications of A 4209 and A 4206 are both plausible. 


7. If 44152 is a blend of *G,,, 5,—-"H5, then the two 4G levels must be relatively close 
in energy, and in that case A 4906 and A 4918 cannot both be correctly identified. If A 4152 
is a blend, the identification of A 4918 is a problem. If A 4152 is not a blend, it should be 
accompanied by a companion of comparable strength. There is no obvious way out of this 
difficulty. 


The star RR Telescopii in 1954 was undoubtedly (16) the best object in 
the sky for a search for [Fetv] lines. We have examined Thackeray’s list of 
wave-lengths and intensities, with special attention to unidentified lines and to 
lines which increased substantially in intensity from 1953 to 1954. ‘The results 
of this search are contained in Table IV, and some detailed comments are 
appended to that table. There are a number of lines in RR Tel whose wave- 
lengths agree with those of certain predicted lines, within the probable limits 
of error of the calculations, and for which the observed and predicted intensities 
are in tolerable agreement. There is a good chance that some or all of these 
lines are due to [Fetv]. There is of course no certainty about the suggested 
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identifications: only a laboratory analysis can decide as to the correctness of the 
suggestions. An attempt to use the data of Table IV to provide a partial term 
analysis of Fev was not very successful; the number of identifications made 
is inadequate. 


5. Conclusion.—This paper has presented data on the energy levels and 
transition probabilities of forbidden lines of Fe 1v, and a number of identifications 
with lines in RR Tel have been suggested. ‘The numerical data in this paper 
are of rather lower accuracy than data in earlier papers by the writer because 
the calculations are not based directly on observed energy levels. Nevertheless, 
even if the suggested identifications are not correct, the results of this paper 
should give a reliable indication as to the transitions of importance in ions with 
3d° configurations. It is hoped that the work described above will stimulate 
some experimental spectroscopist to make a renewed attempt to produce and 
analyse the Fe iv spectrum. 
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Note added in proof.—Since the above paper was written Professor Edlén 
has kindly sent to the author further details of the identifications referred to by 
Thackeray (16). Analysis of measurements in the vacuum-ultraviolet which 
had been made in connection with work on Fe enabled six quartet terms to 
be established. Among these are the 4G, ‘D and ‘F terms of the 3d° configuration, 
and these allowed the identification of individual forbidden transitions. 
Professor Edlén has given permission for his results to be published here. The 
energy levels found by him are listed in Table V; to facilitate comparison they 
are in the same order as in Table II. In all cases the calculated term intervals 





TABLE V 


Relative energy levels of Fe IV derived by Edlén from a partial analysis of the 
vacuum-ultraviolet spectrum 


Term J Energy Term J Energy 
3d5 4G 33 a 3d*4s 4D 2} 106100, 
“4 9 _ 6 34 106601 
53 ° 
3d"4p ‘F 2} 164093, 1g 
33 164311, 
3d° 4D 2} 6698 44 164604793 
i — 163 - 
32 6535 
3d*4p *D 1} 169843, 6 
3d 4F 24 20599 _ 2} 170089, 
34 2045514 3h 170366777 
43 20380 


agree with the observed ones within the accuracy of 50 cm~! mentioned in the 
paper, and the intervals between terms agree within 500 cm~!. The wave-lengths 
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of the forbidden lines in the 4G—*F and *D-*F multiplets of the 3d° configuration 
have been calculated from the observed energy levels and are listed in Table VI. 





TaB_e VI 
Wavelengths of forbidden lines calculated from observed energy levels of the 3d° 
configuration 
Transition A Identification in RR Tel 
Int (1954) 
‘GF 34-2} 4867°4 4868-2: 8: 
43-23 4866°5 
33-33 4901°8 4900°2 8 
43-32 4900°8 4900°2 8 
53-33 4887°4 
33-44 4919°9 
42-43 4918-9 4918"1 4 
53-41 4905°4 4906°2 12 
‘DF 23-2} 7191°7 7189 3 
33-23 7108°4 
23-33 7267°0 
33-33 7181-9 
23-43 7306°9 
33-43 7220°8 7221 3 


In Table VI we have also included those identifications in RR Telescopii which 
are now firmly established. It is gratifying to note that four of the suggested 
identifications in the *G-*F multiplet have been confirmed. The problem 
mentioned in note 3 to Table IV is resolved: 47189 and A7221 belong to *D-*F. 
There remains the problem of identifying *P,,-?D,,, which if the computed 
transition probabilities are correct ought to be of considerable strength; 
possibly it contributes to 47171. 
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THE PHOTO-ELECTRIC LIGHT CURVE OF RT SCULPTORIS 
G. G. Cillié and E. M. Lindsay 


(Received 1958 August 13) 


Summary 


With the 60-inch Rockefeller reflector of the Boyden Observatory photo- 
electric light curves have been obtained in blue and yellow light for the 
eclipsing variable RT Sculptoris. ‘The two maxima are unequal and are not 
flat-topped in any way. ‘There has been a slight shortening of the period of 
this star during the past 50 years and assuming the decrease to be linear we 
find the following light elements which best fit the observations: 


Min. = JD 242373651145 + (0%-51156702 — 1:0568 x 10 ME)E. 





The variability of the star, RT Sculptoris, CoD —26°179, CPD — 26°43, 
position oo! 31™308, —26°13'-5 (1900) was discovered on Harvard plates by 
Ida Whiteside (1) who gave the following light elements: 

Min. = 24117364114 +.04511574E. 
The star varied between photographic magnitudes 9-64 at maximum and 10°47 
and 9:84 at the two minima. The first maximum was flatter and of longer 
duration than the second and it was considered that the star was of the 8 Lyrae 
type. Miss Whiteside’s material consisted of 276 photographic observations 
over the period 1891-1907, and since her first observation was on J D 2411736°436 
her epoch must have been found by using the determined period and extrapolating 
backwards. From 1891.to 1899 only thirty estimates were available. The 
residuals from these and a few others were larger on the whole than from the 
rest. It was considered, however, that they did not show any systematic 
difference which would indicate a changing period and that they were due to 
the plates having been taken with a different instrument. 

Shapley (2) discussed Miss Whiteside’s observations, determined orbital 
elements, and gave ranges of primary and secondary minima of o™-71 and o™-20 
respectively. 

In 1925 Schilt (3) published some further observations of this star from 
measurements with the thermopile-microphotometer on 326 plates taken by 
Hertzsprung on sixteen nights at the end of 1924. The quality of many of the 
plates was rather poor due to spots and irregularities in the film and the 
determination of the magnitudes of the variable was uncertain in some cases 
as a consequence of the large irregularities in the galvanometer readings of the 
comparison star. Schilt found that in comparison with the earlier observations 
the secondary minimum had shifted 04-026 towards the following primary 
minimum. He also found that the first maximum is higher and broader than 
the second one and that the descent towards secondary minimum is steeper 
than the rise from secondary minimum, a result that is confirmed by our 
photo-electric observations. Schilt determined a ‘‘normal minimum’’ at 
JD 2424147298. This minimum, together with that of Miss Whiteside’s, 
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gave an “‘improved period”’ of 0%-51156935. Thus his light elements for 
RT Sculptoris are: 
Min. = 2411736°114 + 0"51156935£. 
In 1928 Dugan (4) gave the following light elements: 
Min. = 2423736751145 +.09-511569633E. 

These were determined by using Miss Whiteside’s original epoch and two 
well-determined minima by himself, using the polarizing photometer on 1923 
December 7 (JD 2423761) and 1924 November 26 (JD 2424116). The 
dispersion in Dugan’s observations were rather larger than usual and their 
distribution not normal. Systematic errors of some unknown type were present. 
Dugan concluded that any irregular variation in the comparison star or in one 
of the components of the variable would be pretty well ironed out by observations. 
In attempting to compute the physical dimensions of the system he found it 
impossible to obtain elements which predict the light curve satisfactorily. 


TaBLe I 
Primary minima of RT Sculptoris 


Min. = 242373651145 Min.=2411736°114 


Observed +04-51156471E +(0°51157183 
Date of run time of prim. Wt. — 10245 x 10-!F")E’ 
min. Calc. min. O-C Calc. min. O-C 
2AZ ove BAZ ons d CT ae a 
1952 July 28/29 4222°5659 4222°564 +O°O0OII 4222°567 —o'0018 


1954 September 5/6 4991°44630 
September 28/29 5014°46611 
October 14/15 5030°32572 
October 15/16 5031°35012 
October 16/17 5032°37130 


4991°4466 —0:0003 4991°4465 +0-0001 
5014°4670  —o-0009 5014°4666 +0:0004 
5030°3256 +0°0001 5030°3251 +0°0006 
5031°3487 +0°0014 5031°3483 +0-0018 
5032°3718 —0°0005 5032°3713 0°0000 


Un nN AW 


In communications to Dugan, Adams and Joy gave the spectral classification 
of RT Sculptoris as A5n and they found that the two spectra were of nearly 
equal intensity, that of the star receding between primary and secondary minima 
being perhaps a little fainter. 

As Dugan pointed out, and as will be seen from the above, all observers have 
remarked on some peculiarities of the observations. It is also to be noted that 
any ‘‘improved’’ periods have been based on Miss Whiteside’s original epoch 
and obtained on the’assumption that the period is constant. In view of the 
uncertainties with regard to this system, the star was put on the programme of 
the Linnell-King photo-electric photometer attached to the 60-inch Rockefeller 
reflector of the Boyden Observatory in 1951. All told, 1986 comparisons—1064 
in blue light and g22 in yellow light—were made in the course of 1951, 1952 
and 1954 with the nearby star CoD — 26°193 (9™-1) = CPD — 26°49 (g™-8). 

Well-observed primary minima of RT Sculptoris were obtained on six nights, 
and using Dugan’s fundamental epoch we obtain the best mean-square fit with 
the following light elements: 


Min. = 2423736°51145 +0%-51156471E. 
In Table II we summarize the various determinations of the period of 
RT Sculptoris up to date, the first two being determined with the fundamental 
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epoch of Whiteside and the last two with that given by Dugan. From Table II, 
as indeed we suspected from rough plots of our observations, it is evident that 
the period of RT Sculptoris has been slowly decreasing during the past 50 years. 
So far as one can judge, the decrease has been steady and roughly linear with 
time. Hence we try to fit a formula of the following type to the observed primary 
minima, 
Min. = 7) + (P)—a@E’)E’, 
where 7) is some fundamental epoch, a a constant, E’ the number of cycles 
calculated from 7, and Py may be termed the “‘original period’’. Using 
Whiteside’s fundamental epoch, Schilt’s ‘‘normal minimum’”’ and our own 
minimum of 1954 October 16/17, we obtain the light elements: 
Min. = 2411736-114 + (0°51157183 — 1:0245 x 10-1E’ )E’. 
The minima predicted with this formula at the time of our observations and the 
residuals when these are compared with the observed times are given in the 
last two columns of Table I. The actual period of RT Sculptoris at any 
particular time (or epoch E’) is given by 
Py —2aE’ =04-51157183 — 2:049 x 10-1 E’. 
At the time of our observations in 1954 October, the period of RT Sculptoris 
was 0751156249 and we use this to compute a ‘“‘normal minimum”” at that 
time. It comes out at JD 2435019°58244 at the epoch E’=45514 calculated 
from Whiteside’s original epoch. 


TaBLe II 
Determinations of the period of RT Sculptoris 


Observer Period 
d 
Whiteside O°511574 
Schilt - "51156935 
Dugan *511569633 
Cillié and Lindsay *51156471 


We now recalculate the constants in the formula for the minima of 

RT Sculptoris, namely, 
Min. = 7," + (P,”—a@’E")E’, 

taking as 7)"=2423736°51145, which is Dugan’s fundamental epoch, and 
making use of Whiteside’s fundamental epoch and our own “‘ normal minimum ’’. 
This gives for RT Sculptoris, 

Min. = JD 242373651145 + (0451156702 — 1:0568 x 10-"E" JE”. 
Schilt’s normal minimum gives an O-C residual of —o*-o017 with this 
formula, which is not surprising and about the same residual found by Dugan 
when he tried to fit Schilt’s normal minimum to his own observations. 

Our observations of the primary minima of RT Sculptoris prove that the 
period of the star has in all probability been steadily decreasing during the past 
50 years. No doubt this change of period has something to do with Dugan’s 
difficulty in interpreting the light curve satisfactorily. Secondary minima of 
RT Sculptoris fall half-way between primary minima according to our 
observations. The weighted mean phase of secondary minimum is 0-4984. 
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In Tables IJI and IV we give the normal points in blue and in yellow light 
for the light curve of RT Sculptoris. These have been computed in the usual 
way by using our final light elements. The period changes very slowly, so that 
even for the duration of our observations the period of the star changed very 
little and we had no difficulty in combining our observations into normal points. 
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Fic. 1.—Blue light curve of RT Sculptoris. 
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Fic. 2.— Yellow light curve of RT Sculptoris. 
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TaBLe III 
Blue light curve of RT Sculptoris 
Phase M.e. of Phase M.e. of 
(period) Amp normal pt. Wt. (period) Amp normal pt. Wt. 
070008 +o™-470 +0™-o902 43 0°5621 —0-9073 +0™-902 18 
‘OIQI “426 "003 39 *5792 "099 "002 18 
0385 296 "002 24 6051 "134 003 26 
"0585 155 "003 28 6335 ‘I77 "003 28 
"0829 +0°024 "005 21 6698 "205 "003 26 
*1006 — 0060 005 25 "7193 234 "004 15 
"1214 "104 "005 16 "7395 ‘241 "003 27 
*1396 "149 "003 23 *7582 239 "003 24 
“E717 208 003 23 -7796 228 "002 26 
"2093 258 "002 20 “8000 "224 "003 33 
*2480 ‘289 "003 20 “8195 *208 "003 20 
"2912 285 "003 15 8389 185 "002 2 
°3245 267 ‘002 17 8613 ‘149 "002 e36 
°3534 "250 "002 27 *8798 "105 "002 37 
+3856 “2579 "003 28 -9g018 —0'043 “002 40 
"4152 "145 "004 21 ‘Q177 +o-018 003 25 
"4399 "084 004-24 ‘9397 "145 004 34 
-4606 OSI "004 18 9589 278 003 36 
+4818 ‘020 ‘002 28 0'9783 +0°422 +0°003 37 
“5011 ‘006 "005 19 
“5213 "024 003 19 
0°5402 —0°'040 +0'003 19 
TABLe IV 
Yellow light curve of RT Sculptoris 
Phase M.e. of Phase M.e. of 
(period) Amy normal pt. Wt. (period) Amy normal pt. Wt. 
0°0001 +oM-581 +0™-003 31 0°5402 +0-096 +o™M-004 18 
‘0188 "539 003 27 +5629 066 003 16 
"0394 408 +903 19 *5799 +0033 “002 18 
"0595 292 "004 22 6066 —0°005 "003 22 
“0804 "158 "005 15 6343 "049 "002 19 
-0998 ‘070 "004 23 6628 ‘070 "002 25 
“r2I3 +0°008 "005 13 *6976 ‘099 "003 16 
"1399 — 0033 "003 18 "7179 “E72 "003 15 
*1703 076 ‘003 22 "7403 ge "003 24 
*2099 “$22 "003 20 7583 ‘110 003 19 
-2466 146 ‘003 19 “7800 093 "004 22 
"2900 "145 "003 18 “8000 *O0gO "003 28 
3259 ge "002 17 8196 “084 "004 16 
“9592 “157 “002 24 8401 "055 "002 30 
3851 ‘081 "002 24 *8608 —0'024 ‘002 29 
"4152 —o-'o016 "005 23 8783 +0006 003 28 
"4399 +0°053 “006 21 *QOI5 ‘074 003 32 
"4589 "094 *005—s«d12 ‘9178 ‘139 7004-26 
“4810 ‘129 ‘002 28 -9408 ‘270 "003 ay 
"5012 "149 "005 18 ‘9604 *405 003 28 
0°5201 +0°131 +0°004 16 09789 +0°524 +0°003 30 


Our light curve confirms previous observations in connection with this very 
interesting star. We find primary minimum deeper by o™-48 in blue light 
and o™-44 in yellow light than secondary minimum. The descent into secondary 
minimum is steeper than the subsequent ascent, as first pointed out by Schilt. 
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The maxima do not appear to be flat-topped in any way, but they are slightly 
unequal. ‘‘ Primary maximum ’’, i.e. the maximum following primary minimum, 
exceeds ‘‘secondary maximum’’ by o™-o05 in blue and 0™-035 in yellow light. 
In this, as indeed in other respects, there seems to be a considerable similarity 
between the eclipsing systems RT Sculptoris and V 525 Sagittarii (5). The 
only significant difference we find is the gradually changing period of 
RT Sculptoris. Over the last 50 years the period of this star has decreased 
by something less than 1 second*. 

In conclusion it is a pleasure to record our thanks to the staff of the Boyden 
Observatory for kind assistance while the programme was being carried out. 


* The individual observations will be printed separately and distributed with Armagh Contri- 
bution Number 25. 


Department of Mathematics, Armagh Observatory, 
University of Stellenbosch, N. Ireland. 
Union of South Africa: 
1958 August. 
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RADIO EMISSION FROM ORION 
H. Rishbeth 


(Communicated by F. Graham Smith) 
(Received 1958 July 30) 


Summary 


The Orion region has been surveyed at two wave-lengths by the Sydney 
‘* Mills Cross’ radio telescopes. ‘The radio emission from bright nebulae 
in the region appears to be of thermal origin; but a faint extended source in 
the region may emit non-thermally, and is perhaps associated with an extensive 
dark nebula. Radio and optical observations of the Orion Nebula, M42, 
are compared with the aid of a simple model. 





1. Introduction—The Orion region has long been recognized to be of great 
astrophysical interest, and a large literature dealing with it is now in existence. 
The outstanding object in the constellation is the Great Nebula, M42; many 
other emission nebulae excited by early-type stars are found in the vicinity. 

The region offers a good opportunity for the comparison of optical and radio 
observations, both of which enable the ‘‘emission measure’’ to be determined 
provided that the radio-frequency radiation is of thermal origin. Differences 
between the optical and radio emission measures may thus be attributed to optical 
obscuration or to non-thermal radio emission. Evidence already available (e.g. 
(1)) suggests that emission nebulae emit thermally at radio wave-lengths, in 
contrast to other types of radio source which can only be explained in terms 
of a non-thermal emission mechanism. 

Many radio observations of the Orion region have been made in recent years, 
and isophotes at decimetre wave-lengths have been published (2, 3, 4). This 
paper describes the results of high-resolution surveys at metre wave-lengths, 
obtained with the Mills Cross radio telescopes of the C.S.I.R.O. Radiophysics 
Laboratory near Sydney, N.S.W. 

2. Observations.—The principles of the Mills Cross are now well known, 
and as the equipments have been described in detail (5, 6) only the salient details 
will be given here. ‘The 1500 ft cross operates at a wave-length of 3:5 m 
(85:5 Mc/s); the second instrument, whose overall dimensions are about 
3500 ft, operates at 15:2 m (19°7 Mc/s). ‘The crosses are transit instruments 
and possess pencil-beam responses of approximately Gaussian form. ‘The 
Orion region transits about 30° north of the zenith at Sydney, where the aerial 
responses are slightly elongated in declination; the half-power responses for 
the two systems are about 50’ x 60’ and 85’ x 100’ respectively. 

The survey at 3°5m covered an area of over four hundred square degrees ; 
the records were of the scanning type previously described (7). Most of the area 
was covered twice, and a number of simple records of greater sensitivity were 
obtained to provide additional information about particular features. ‘The 
reduction procedure outlined in another paper (8) was employed in the 
construction of isophotes, which are presented in Fig. 1. ‘The discrete sources 
are substantially those of Mills and Slee (7). 








x 
~ 
YS 
SS 
= 
J 
& 
= 





*,6— pup .9+ 
suoyouysap puv SSSo pnv SoSo suotsuarsp qysia Aq pausof apSuvjra4 ayy w!yji2 pautwjuor st ((q) § uotjrzag) Sura 11 mann 
(6) SOU Osalg 2oupoys 243 wosf uayny aav (pappys) avjngau yavp puv (pI}jop) suoisas 1] FY “wuvaqg jvlaav fo anoquo? sanod-fjvH = Vy 

*a24nOS papuayxy = © *dI4NOS parwjosr4ut_ = + 
“uaz0ys aav saranos pogigutag ayy fo saanjpsaqiuay potsav yvad ayy 3ngq ‘pajqtwu0 a4V $2I4NOS 3]24ISIP 
SUIPUNOAANS SANOJUOD PIso]) AuD]Yy "My, COOL yu = *sagvuIps007 OSO61 = *Y15ua]-aav2 wWS.€ yw uo0lsa4 UOtG_ ay} fo sajoydos—T—'I “D1. 


OrrO 00S0 ozs0 a, Orso 0090 0290 


OC 
yor 



































— =e 


0505 aud 0555 and declinations 


m& WECMON 5S \O)) tS Contained within the rectangle formed by right ascensions 


acer reerte SD SL Li Fiti 


+6° and —9°. 











Radio emission from Orion 593 





No. 6, 1958 


At 15m wave-length, broadcast interference and ionospheric effects were 
troublesome and fewer records were obtained. Corrections for ionospheric 
refraction were estimated from daily observations of intense sources, and have 


been applied to the isophotes of Fig. 2. 
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Fic. 2.—Isophotes of the Orion region at 15 m wave-length. 1950 coordinates. Unit 1000 °K. 
A is the half-power contour of the aerial beam. 
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3. Theory.—The theory of thermal radio emission by free-free transitions 
in ionized hydrogen has been discussed by Piddington (10) and Mills, Little 
and Sheridan (11). A few results are given here. 

The optical depth 7 of a region of ionized gas is proportional to its emission 
measure EF, defined as E=f N*.dl, where N is the electron density (in cm~*) 
and / is the path length through the gas (in parsecs). ‘The equation is: 

7+=A29(A). ET? (1) 
where g(A) is a slowly-varying function of wave-length A and 7’, is the electron 
temperature, which is of the order of 10*°K in emission nebulae. " Typical 
values of the quantity E/r for various radio wave-lengths are given in ‘Table I. 


TABLE I 
Values of E/r 
E/r 
Wave-length T,=10000°K T,=20000°K 
15'2m 7°9 X 107 21 xX 10° 
9°65 1°6 x 104 4°4x 104 
1°9 5°8 x 10! 1°6 x 105 
0°5 gt x 105 2°5 X 10° 
orl 2°6 x 107 7°0X 107 


For optically thin objects, sufficiently extensive to fill the aerial beam, a 
simple relation connects the average emission measure and the increase of aerial 
beam temperature (A7\,) due to the source. It may easily be shown that under 
these conditions A7T.,=7.'T,, and substitution in (1) of the values 7',= 10000 °K, 
A=3'5m and g(A)=5:0 gives 


E=1-6AT,,. (2) 
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At 3°5m wave-length, the emission nebulae in Orion are seen against a 
background of cosmic radio emission, whose brightness temperature is about 
1000 °K. ‘The above equation must be modified to take account of this, and 
becomes 


E=1-8AT,, (3) 


The radiation from an emission nebula, small compared to the aerial 
beamwidth, is more conveniently specified by its flux density S, given by: 
° 2 
s | E.dQ= —— (4) 
where dQ, dV are elements of solid angle and volume respectively. D is the 
distance from the solar system. If the ionization is contained in a homogeneous 
sphere of radius a, as in the model of an Hi region proposed by Strémgren 
(12), then 
2 3 
[Se- (5) 
3D 
where /*= N*a*, J then depends upon the temperature of the exciting star: 
it is the quantity tabulated by Strémgren as a function of the spectral type of 
the star. 
4. The Orion Nebula 
(a) A simple model of the Orion Nebula, M42.—A simple thermal model of 
M42 can be constructed from available data. Spherical symmetry is assumed 
and the electron temperature is supposed uniform throughout the nebula. ‘The 
distribution of electron density N is taken as a function of radial distance r such 
that 
N?= N,* e-"!4 + N,* ee", (6) 


‘Two terms are necessary in order that the model shall resemble the optical 
appearance of M42. The first term represents the central nucleus of the nebula, 
whose brightness falls off very rapidly with increasing distance from the 
Trapezium cluster of exciting stars. ‘The second term represents the extensive 
outer portion of the nebula. Its radial decrease of brightness is much slower, 
and is specified by the parameter 6 whose value is derived from the optical 
isophotes of Michelson (13). Nz, is chosen to give the observed flux density 
at 3-5 m wave-length, and then the parameters Nj, a of the central portion are 
determined from the measurements at centimetre and decimetre wave-lengths. 

Values of these quantities are given in Table II, and the radial variation of 
emission measure is shown in Fig. 3: the two curves refer to models with electron 
temperatures of 10000 °K and 20000°K. ‘The radio spectra of the models are 
plotted in Fig. 4, together with the observed flux densities which are listed in 
Table III. 

At centimetre wave-lengths the nebula is optically thin, but S decreases 
slowly with decreasing wave-length because of the logarithmic factor g(A) of 
equation (1). The probable error of the point at 3:2 cm is large, as the flux 
density quoted by Kaidanovski and others (17) is equal to the equipment 
sensitivity. A measurement by Haddock (1) suggests that the 3-2cm flux 
density is higher, and closer to the 9-4. cm value. Recent work by Westerhout 
(18) gives the flux density at 22 cm as 520 x 10°-*° w. m.~* (c/s)! + 20 per cent, 
in substantial agreement with the points shown on Fig. 4. 
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TaBLe II 
Model of M42 (equation (6), Section 4 (a)) 


Quantity T, = 10 000 °K T,=20 009 °K 
Assumed distance from the solar system (D) 500 pc 500 pe 
Emission measure at centre 7°6 x 108 1°32 X 107 
N, 6200 cm-* 8500 cm-* 
No 83 cm-* 78 cm-* 
a 0°10 pe 0°092 pc 
(angular equivalent) (41”) (38) 
b 1'o pc IO pc 
(angular equivalent) (6'-9) (6’-9) 
Mass: 
Central portion 30 Mo 30 MO 
Outer portion 400 MO 400 MO 


If Dis changed, then N,, N2 vary as D-'/*; a, bas D; and the mass as D5”, 











pc 
Fic. 3.—Radial variation of emission measure E (cm-® parsec) of the model of M42, for two 
values of electron temperature. For the lower scale ( parsecs) a distance of 500 parsecs from the 
solar system is assumed. 
T,=10 000 °K. 
moececee T,,=20 000 °K. 
said arelecharena Radial brightness variation at 1:-9m 
wave-length for T,=10 000 °K (1) 
and T ,=20 000 °K (2). 


At wave-lengths greater than about 3 m the brightness temperature of the 
galactic background is appreciable in comparison to the electron temperature, 
so the apparent flux density is reduced. At 15 m the nebula should absorb the 
background radiation: the effect would be near the limit of detection of the 
equipment, and is not in fact observed. 

There does not seem to be any non-thermal component in the radiation 
from M42. If such a component existed, it would be most prominent at long 
wave-lengths. An upper limit of 100 x 10-** w. m.~* (c/s)! can be set at 7-°9m 
(G. R. Whitfield, private communication), whereas the upper limit at 15m is 
about 200 x 107° w. m.~?(c/s)~!. With present techniques it is not significant 
that the spectrum of the model does not include both the 1-9 m and 3°5 m points. 
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Taste III 
Radio observations of the Orion Nebula, M42 (IAU source o5S0A) 
Observer Wave-_ ‘Freq. R.A. Dec. Flux 
length 1950 1950 1o"™ 
m Mc/s ‘i a the 0 + +  Ww.m.~?(c/s)-* 
Rishbeth (present obs.) 3°5 85°5 05 32°5 +O2 —5 2743 85 
Baldwin, Archer (private 
communication) 1'9 160 05 32° to2 —5 2444 65 
Seeger, Westerhout, van de 
Hulst (2) 0°75 400 05 32 —5 230 
Piddington and Trent (14) 0°50 600 os 32 —5} 300 


Denisse, Lequeux, Le Roux 














4) 0°33. 900.05 32 —5 30 420 
Hagen, Lilley, McClain (15) 0-21 1420 05 32 —5 20 420 
Haddock, Mayer, Sloanaker 

(16) 0'094 3200 05 33°0 +02 —5 27+5 45° 
Kaidanovski, Kardashev, 

Shklovski (17) 0°032 9400 05 33 —5 27 270 
Visual position of Trapezium 

stars — = 05 32°8 AG 25 = 

_— F500 

gree... ™ 
O>~O 

400} ‘\ F400 

300; s 
5 ° 300 
200; +200 
100; %* F1OO 

0.03 0.1 0.3 | 3 10 


A metres 


Fic. 4.—Observed radio flux densities S of the Orion Nebula, M42, from Table III. Unit 
107 ?%y,m.~*(c/s)-!.. The observation at 7-9 m is an upper limit. The curves show the spectra of 
the model of equation (6) for two values of electron temperature. 
—————  T,=10 600 “K. 
----- eee T ,=20 000 °K. 
jevemataotre These portions of the curves are dependent on the 
background temperature. 


Information concerning the angular size of the M42 source is given by the 
interferometric measurements of Mr S. Archer, of the Cavendish Laboratory, 
at 1'gm. At this wave-length (see Fig. 3) the model has an optical depth of 
unity at a radius near the transition between the central peak and the outer portion, 
and so the observations provide a sensitive test of the validity of the model. 
The aerial spacings used were 27A in the north-south direction and 307A in the 
east—west direction; and the ratio of the corresponding Fourier components of 
the strip brightness distribution is certainly less than 0-2, and probably not more 
than o-1. The values given by the model are o-11 and 0-24 for electron 
temperatures of 10000 °K and 20000 °K respectively. 
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The model with T,,= 10000 °K thus accords with the 1-gm result, but with 
T’,= 20000 °K the effective diameter is too small, and cannot be appreciably 
increased if the radio spectrum is to be preserved. Although this discussion 
relates to an exponential distribution of electron density, the result is not 
critically dependent on the choice of model, and it seems that the electron 
temperature of the nebula is probably less than 20000°K. This merely 
corroborates optical estimates, which are around 10000 °K; but it is possible 
that with further technical progress radio methods will enable narrower limits 
to be set on the value. 

(6) Comparison with optical data.—The overall shape of the M42 model, 
and the density of its outer portions, were chosen to suit optical data, while the 
emission measure of the central portion was determined from high-frequency 
radio observations. Further comparisons with optical data are possible. 

Westerhout and Oort (19) give the emission measure at the centre of the 
nebula as 8 x 10°, in agreement with the radio model. ‘The emission measures 
of Johnson (20) are averaged over areas of 8’ x 18’; similar averages taken from 
the model are comparable with his results, though somewhat smaller. The 
emission measure at the periphery of the model is consistent with the overall 
diameter recorded on long exposure photographs of M42. 

Electron densities in the Orion nebula have been derived by Osterbrock (21) 
from measurements of the relative intensities of the 43726 and A3729 lines of 
Ou. Within 5’ of the centre, the radial dependence of electron density agrees 
well with the parameter a, but the actual values measured by Osterbrock are 
greater by a factor of order five than those given by equation (6) with the adopted 
parameters. In the outer portion, more than 7’ from the centre, the ratio is 
roughly 20:1. However, the model takes no account of irregularities within 
the nebula, the presence of which may affect the comparison of the values of N 
(as given by the optical measurements) with the integral [ N?. dl. 

Another test of the model comes from the application of Strémgren’s (12) 
theory of the excitation of H11 regions by O and B stars. The stars of the 
Trapezium cluster, at the centre of M42, range from O7 to Br in spectral type; 
the value of the ionization parameter J (equation (5)) appropriate to their combined 
effect is found from Strémgren’s Table 5 to be about 100. The theory assumes 
that the ionization density throughout the sphere is uniform; the nebula is far 
from uniform, but it is plausible that the actual value of J should be comparable 
with that given by the simple theory. In fact, the model gives /=60. 

It cannot be claimed that such an idealized model represents the true 
composition of the nebula. It takes no account of the irregular structure*, nor 
of the presence of the solid particles which are responsible for the optical 
continuum radiation of M42. It may however be concluded that the radio 
emission of the nebula can be explained in terms of a simple thermal model 
which agrees reasonably well with optical data. 


5. Other features of the Orion region 

(a) The Horsehead group, IC 434, NGC 2023/4, etc.—This group is situated 
in the vicinity of the stars €Orionis and oOrionis, to the north of the Great 
Nebula, and takes its name from a well-known dark marking visible in IC 434. 


* Much of the apparent irregularity of M42 may however be due to foreground obscuration. 
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It is a prominent source at 3:5 m: the peak is situated at R.A.0537°8 + 0°3, 
Dec. —2° 30’ + 10’ (1950), near the centre of IC 434. The integrated flux density 
is (130+ 10) x 10°-*° w, m.~* (c/s)~!, greater than that of M42. 

The source is very extended: its half-power contour encloses an area of 
34 square degrees, much greater than the aerial beam, and greater than the 
visible extent of the nebula, which suggests that part of the object may be hidden 
from view behind the obscuring matter found in its neighbourhood. The 3:5 m 
isophotes also suggest that the Horsehead and Great Nebulae may be linked by 
a ridge of emission. 

The Horsehead Nebula has also been observed at 21 cm by Hagen, Lilley and 
McClain (15), using a 1° aerial beam; the flux density is 50 x 107° w. m.~* (c/s)~!. 
This appears to be a peak value; if the source has the same angular size at 21 cm 
and 3:5 m the integrated flux densities at the two wave-lengths are similar, as 
would be expected if the radiation were thermal. Westerhout (18) obtains 
95 x 10°* w. m.~? (c/s)~! at 22cm. 

Assuming for simplicity that the radio source is an optically thin H 1 region, 
500 parsecs from the solar system, and that it is roughly spherical in shape and 
of uniform electron density N, the 3-5 m results give N~ 10cm “and a mass of 
about 700 Mo. ‘The emission measure at the centre of this model is 1000, 
with +=0-06: the optical emission measure is given as 4000 (Westerhout 
and Oort (19)). The ionization parameter J (equation (5)) is about 40, which 
is consistent with excitation of this Hi region by ¢ Orionis (type Bo) and 
a Orionis (O9°5). 

Such an Hu region should appear in absorption at long wave-lengths, and 
should produce a dip equivalent to at least 40000 °K on the 15m records. ‘This 
dip is not observed, although the sensitivity should be just enough to detect it. 
If it can be confirmed that the dip is much less than the expected value, it would 
suggest that some non-thermal emission is received from this direction at long 
wave-lengths; though it is possible that its source is unconnected with the 
Horsehead Nebula. 

(6) The H 11 ring of Barnard.—This H11 region is roughly in the form of an 
elliptical ring, 17° x 12° overall, of which the major axis approximates to the 
galactic parallel of latitude — 18°: the angular thickness of the ring is nowhere 
greater than about 1°. It encloses M42 and the Horsehead complex. The brightest 
parts of Barnard’s Ring are in the following sector, towards the galactic equator : 
the preceding parts arg extremely faint, and are barely visible in photographs by 
Lower (22) in which the brighter parts are heavily over-exposed. ‘The maximum 
optical emission measure is about 2400 (19). 

It is probable that the visible ring is the periphery of a three-dimensional 
structure. ‘This interpretation is supported by observations of 21 cm line 
radiation by Menon (23), who detected at points inside the ring three concentra- 
tions of neutral hydrogen with different radial velocities. His results imply the 
existence of a central concentration of neutral hydrogen, surrounded by an 
expanding shell; assuming this shell to be spherical, he computes its thickness 
and inside radius to be 29 pe and 38 pc respectively, and the average density 
as 3:8 atomscm~*. The duration of expansion is estimated as a few million 
years. 

The following part of Barnard’s Ring is visible as a ridge on the 3-5 m 
contours, extending between declinations +4° and —7°. The resolution is 











XUM 








— Ee Tt. ct 





XUM 


No. 6, 1958 Radio emission from Orion 599 


inadequate to show much detail, but the observations yield values of { E. dé, 
the integral of the emission measure along a declination scan across the ridge. 
The corresponding optical quantity cannot be accurately evaluated in the 
absence of isophotes but a very rough estimate gives a result three times the 
3°5m value. At 15 m a trough of absorption coincides in position with the 
bright parts of the ring: the maximum depth is about 30000 °K, and the 
estimated value of { E.d@ is consistent with the 3°5 m result. 

For comparison with the 21cm results, the 3:5 m observations can be 
related to a model of a uniform transparent shell. With allowance for aerial 
smoothing, the thickness of this H 11 shell is found to be <15 pe: the maximum 
emission measure is about 1000, leading to an electron density N = 4 in the 
bright portions of the ring. ‘The preceding portions are unobservable at 3-5 m, 
and for these N<1°5. 

Application of Strémgren’s theory shows that the radiation of stars inside 
the H 11 ring is insufficient to produce the observed ionization (e.g. (23)). Instead, 
the radio observations suggest the presence of an H1 shell expanding into a 
relatively stationary medium. At the boundary, a thin shell of H 1 ionization 
is formed by collisional excitation: this is most intense in low galactic latitudes, 
where the density of the surrounding medium is expected to be greatest. 
Comparison of 21 cm and 3°5 m data shows that in these regions the percentage 
ionization may be high. Comparison with optical data suggests that there is 
no appreciable non-thermal component in the 3°5 m radiation from the edge of 
the shell. 

(c) An extended source in Orion—‘‘ Orion X’’.—The 3:5 m isophotes reveal 
an extensive area of faint emission covering much of the north-following part 
of Orion. A scan across the region at declination —2°17’ is shown in Fig. 5, 
and it is evident that the received radiation is stronger in the centre of the Orion 
region, at about R.A. 0545—between the peaks at 0538 and 0552, due respectively 
to the Horsehead Nebula and Barnard’s Ring—than outside the region, as at 
0500 and o605. ‘This effect is visible on all scans between declinations + 5° 
and —7. The contributions to the profiles of the other features can be 
estimated, as indicated in the figure: the separation is of course very rough, 
but if it is made for each scan an approximate contour diagram of the excess 
radiation can be built up. ‘The maximum excess temperature is 150 °K, and 
the angular extent is estimated as 7° in R.A. by 12° in declination. The source 
will be referred to as ‘‘Orion X’’: it is perhaps associated with Barnard’s Ring, 
but seems to extend outside the north-following boundary of the Ring. 

Several observers have reported a faint extended source in this position 
(Table IV). ‘Though the spectrum has not been measured accurately, there 
seems to be little variation of flux density with wave-length. 

Such a spectrum is consistent with thermal radiation from an optically thin 
H 11 region, and this interpretation appears reasonable in view of the presence 
of faint Ha emission over much of the Orion region. Application of the 
assumptions previously made (Section 5 (a)) to the radio source gives E ~ 300, 
N = 2, and a mass of 10°Mg. The ionization parameter J of equation (5) is 
about 100—equivalent to the radiation of sixty stars of type Bo. It does not 
seem feasible to attribute the ionization of this object to the stars actually found 
in the region, especially as the emission nebulae previously discussed account 
for the ionizing radiation of several of the brightest stars; and it may be shown 
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that this deficiency remains if the radio source is supposed to be distinct from 
the Orion complex, and to lie closer to the solar system. If the object is thermal, 
its excitation must be collisional, and not due to stellar radiation. 

Long wave-length observations, however, cast doubt on the classification of 
Orion X asa thermal source. It is seen in emission at 7-9 m, at which wave-length 
a faint Hit region would scarcely be prominent as the background brightness 
temperature in this part of the sky is close to 10000 °K. Moreover, if Orion X 
were an HII region it would have an optical depth of about 0-4 at 15 m and would 
thus produce appreciable absorption. The situation is complicated by the 
absorption due to Barnard’s Ring; however, comparison of 15 m radiation from 
Orion with that from neighbouring areas in lower galactic latitudes yields no 
evidence of absorption in Orion X. 






¥ 400 °K 








' { ' 
0600 0530 O500 
R.A. 1950 


Fic. 5.—Scan through the Orion region at declination —2° 17’ (A=3°5 m), showing features :— 
H=Horsehead Nebula, IC 434, etc. (Section 5 (a)). 
B=Barnard’s H 11 ring (Section 5 (b)). 
X = Extended object ‘‘ Orion X”’ (Section 5 (c)). 
S= Discrete sources (Mills and Slee 05-03, 05-06). 
The radio brightness increases eastwards from R.A. 0600 towards the galactic equator. 
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Observations of an extended radio source in Orion ‘“‘ Orion X”’ (Section 5 (c)) 


Observer Wave- R.A. Dec. Diameter Peak ‘Total Flux 
length 1950 1950 R.A.x Dec. temp. 10774 w.m.~? 
m aoe : ; °K (c/s)— 
Blythe (24) 79 of 40 ~+=+1 IOX5 3000 13 
Rishbeth (present observations) 3°5 0§ 40 -—I 7X12 150 10 
Bolton, Westfold, Stanley and 
Slee (25) 3°0 05 42 ° 10X5* — 15 
Ko and Kraus (26) I'2 05 42 ° 9x9 40 (~ 14)T 
Piddington and Trent (14) o5 0542 —I1 5x5 7 (~ 10) 


* Major axis parallel to galactic equator. 
+ Author’s estimate. 


These considerations do not conclusively disprove that Orion X is a thermal 
source, but they suggest that a non-thermal component may be present. It is 
therefore interesting that the extent of the radio source is similar to that of the 
dark nebula charted in the Skalnate Pleso Atlas (g) and indicated in Fig. 1. 
Such concentrations of interstellar dust in the vicinity of bright stars might 
afford suitable conditions for the production of non-thermal radiation. 

(d) Other objects in Orion.—The emission nebula M43 lies about 30’ to 
the north of M42, and its dimensions are 42’x 26’. At 3:5 m it is not 
resolved from M42, but a northward extension of the contours of the latter 
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is interpreted as the contribution of M43, which is assigned a flux density of 
(20+5)x 10-*% w. m.~*(c/s)-!, a quarter of that of M42. With simplifying 
assumptions as before, the average electron density of M43 is found to be 
25 cm~*, the mass of order 50Mo, and the central emission measure 3000. 
The last figure seems too low for an object as bright as M43, but the total 
ionization gives J ~ 20, quite consistent with excitation by the B1 star 42 Orions. 
The optical brightness of M43 shows a rapid radial decrease from the central 
peak, and so the uniform model used to derive the above figures is not very 
appropriate, and gives order-of-magnitude estimates only. 

Amongst the faint H regions in Orion is a group a few degrees to the west 
of M42. Some of these nebulae are comparable in extent to the 3-5 m aerial 
beamwidth. With assumed emission measures of order 500, their flux densities 
should be about 10 x 10-** w. m.~* (c/s)~1, great enough for detection. A source 
of flux density of 14 x 107° w.m.~*(c/s)~1, at R.A. 0513, Dec. —2°, may be 
identified with one Hi region; whilst on examination of the isophotes the 
nebula at 0512, —5° can be assigned the very small flux density of 
3 x 10-*6 w. m.~* (c/s)~!, corresponding to an increase of 50° K. 

A weak source at 0540, — 54° coincides with a mass of filamentary structures 
prominent on the Palomar Sky Survey red-sensitive plate but barely discernible 
on the corresponding blue plate: this object is two degrees east of M42. The 
remaining discrete sources in the region do not seem to be associated with optical 
features and most of them are doubtless merely aligned by chance with the 
Orion complex. 

A 15m source is thought to exist at R.A. 0550, Dec. —7° (see Fig. 2), but 
requires confirmation. It lies near a bright part of Barnard’s Ring. If real, 
this source is non-thermal, on account of its high brightness temperature. 

3°5 m records of the A Orionis emission region, which is outside the area 
covered by Fig. 1, have also been examined. A weak extended source is present, 
but it does not correspond well to the optical features. The 3°5m emission 
measure is about 800 at R.A. 0520, Dec. + 10°, and 300 at 0528, + 10° (1950). 
The optical emission measure near the latter position is 3100 (20). 

6. Conclusions.—Several emission nebulae in the Orion region have been 
observed at a wave-length of 3-5 m, and comparison with optical data indicates 
that the radio emission is thermal. In the case of M42, it is found that a simple 
thermal model can be made to agree with a variety of radio and optical observations. 
It is suggested, however, that the Horsehead group, IC 434 etc., may emit a 
non-thermal component of radiation at the longest wave-lengths. 

In several cases the radio emission measures are less than the optical values. 
A similar result was found by the author for H1 regions in the Vela-Puppis 
region (8). On the other hand, good agreement between optical and radio data 
has been found by previous observers, e.g. by Haddock (1) for the bright galactic 
nebulae in Sagittarius, and for the Rosette Nebula by Mills, Little and Sheridan 
(11). Many of the author’s observations refer to very faint objects and the 
differences can largely be attributed to selection of the optical data; the optical 
emission measures of a nebula apply to the brightest parts, the radio observations 
to averages over finite areas. 

A large faint extended source, referred to as Orion X, is detected in the region. 
Most of the radio data suggest that it may be an H1! region with collisional 
excitation; but long wave-length observations, though not conclusive, favour 
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a non-thermal origin. The position and extent of this object suggest some 
correspondence to the extensive dark clouds in Orion, and to the Hi! ring of 
Barnard. 
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IMPROVED MEASUREMENTS OF THE POSITIONS OF 
17 INTENSE RADIO STARS 


B. Elsmore 


(Communicated by M. Ryle) 
(Received 1958 July 28) 


Summary 


Some observations at a wave-length of 1-9 m, which have been aimed 
at improving the accuracy of the declinations of the sources, are based 
on the times of transit across two interferometers having inclined axes. 
The method of using the Cambridge Radio Telescope in this way to form 
a crossed axis interferometer is described, and the accurate positions of 
17 sources are given. The positions derived in the third Cambridge Survey 
are included for comparison. 





1. Introduction.—The Cambridge Radio Telescope (1) has recently been 
used at a wave-length of 1-9 m for a new survey of radio stars. A detailed 
account of this survey, which will be referred to as the third Cambridge Survey 
(3C survey) is to be published shortly (2) but as some preliminary results are 
available it was thought worth while to check the positions of a number of the 
intense stars by an independent method. ‘The present paper describes this 
method and gives the accurate positions of 17 sources. 

In the past, determinations of the times of transit have provided measurements 
of mght ascension more accurate than those of declination. The present 
observations, which are aimed at improving the accuracy of the declinations, 
are based on measurements of the times of transit across the principal planes 
of two interferometers liaving inclined axes (3). The method of using the 
Cambridge radio telescope in this way to form a crossed axis interferometer 
is described in the next section, and the results and a discussion are given 
in Section 3. 

2. The methods of observation.—The declination of each radio star was 
determined from two observations at a wave-length of 1-gm, taken on 
successive days, using first one diagonal pair of aerials of the Cambridge radio 
telescope and then the other. The pairs of aerials were connected in turn to 
a phase-switching receiver (4) to comprise two separate interferometers with 
their axes inclined at a small angle (10°5’), as shown in Fig. 1. The difference 
between the times of transit 7, and 7, in Fig. 2 across the two fictitious 
meridians is a measure of declination; the accuracy of the measurement is 
greatest for sources near the zenith Z and is least for sources near the horizon. 

The reception polar diagram of the individual aerials limits the observation 
to an interval near the central meridian transit, indicated by 7, in Figs. 2 and 3, 
which corresponds to axis O in Fig. 1. The transit 7, may be derived, however, 
by subtracting from the measured time ¢, in Fig. 3 (6) an integral number N 
of fringe periods; the number N is derived from a knowledge of the approximate 
position of the source. 

The periodicity of the interference fringes (shown dotted in Fig. 3 (6)) 
depends upon the mean wave-length, which is difficult to measure precisely 
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when the receiver pass-band is wide and possibly asymmetric, and errors in 
this measurement become serious when N is large. However, by inserting 
between one aerial and the receiver an additional cable, whose electrical length 
may be exactly determined, the fictitious meridian may be moved through a 
known angle to such a position that N is small for any given source; it then 
becomes a small circle parallel to the original meridian. In the example shown 
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Fic. 1.—Showing the Cambridge Radio Telescope connected to form two separate interferometers 
with axes 1 and 2. 


Celestial 
Equator 








Fic. 2.—Celestial sphere illustrating the fictitious meridians of a crossed axis interferometer. 
The times of transit of a star are T, and Ty. 


graphically in Fig. 3 (6), an additional cable of electrical length three nominal 
wave-lengths (3 x 188 cm), inserted between the SW aerial and the receiver, 
will move the collimation plane from 7, to a measured point near ¢,; similarly, 
in Fig. 3 (c), 7, will move to near t, when the additional cable is inserted in 
the path between the SE aerial and the receiver. In the interval between ¢, 
and ¢, any small difference between the mean wave-length and the nominal 
wave-length is now unimportant. 

In order to calibrate the system and to ensure that there is no difference 
in the phase angles introduced by two pre-amplifiers, a c.w. signal at the 
nominal wave-length of the receiver is passed down two equal cables and into 
the pre-amplifiers, the outputs of which are connected through equal cables 
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to the phase-switching receiver as in Fig. 4. An additional cable, whose length 
is determined by the declination of the source to be observed as described 
above, is inserted into one path near the receiver and adjustments are made 
by the addition of further small lengths of cable so that the two paths differ 
by a predetermined exact number of wave-lengths. The phase angle of the 
individual aerials is determined from observations of radio stars which have 
been identified with optical objects and whose positions are accurately known. 


(a) Pl mn 











To 
Fic. 3. 

(a) The reception polar diagram of the individual aerials. 

(b) The full line indicates the recorder output produced by the transit of a star across the 
NE|SW interferometer. The dotted line shows the maxima and minima and T, is the time of 
transit across the fictitious meridian. (Time increases towards the right.) 

(c) T, is the time of transit across the fictitious meridian of the NW|SE interferometer. 


Preamplifier Receiver Preamplifier 
WS 


>. ~~ ft. 
—_ 











Signal Generator 


Fic. 4.—The arrangement used for calibration. 


The times ¢,, ¢, may be determined accurately from the times of zero 
deflection; the output of the receiver is therefore applied to a photographic 
recording system, designed to operate only near the times of crossovers and 
to record on 35 mm film moving at a speed of one inch in six seconds. 
Fig. 5 shows a section of film produced during one of the transits of 2C.34. 
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The right ascension of each source is derived from the mean of the times 
of transit across the two skewed axis interferometers. 











Fic. 5.—A part of the film produced by the photographic recorder, 
during one of the transits of 2C.34. 


3. Results and discussion.—Errors in position may arise from four causes: 
(i) instrumental effects; (ii) limitation of reading accuracy (i.e. signal-to-noise 
effects); (iii) confusion effects due to weak adjacent sources; and (iv) refraction 
in the terrestrial ionosphere and atmosphere. 

(i) Instrumental effects.—The present determination of the declination of a 
source depends upon measurements of the times of transit, in contrast to the 
method used in the 2C and 3C surveys, which required measurements of the 
relative amplitudes of the records on successive days. Measurements have 
shown that the overall phase stability of the aerials, pre-amplifiers, cables and 
receiver remains constant to within 10, corresponding to 5cm of path 
ditference, over periods of several months. ‘The overall gain stability has been 
found to vary from day to day by up to 15 per cent, particularly in the presence 
of heavy rain or frost on the aerial feeder system. 

The errors in declination caused by such gain variations in the 2C and 3C 
surveys are likely to be greater than those caused by phase variations in the 
present method. For intense sources, for which the instrumental effects are 
more severe than those caused by the confusion effects of weak adjacent sources 
(see (iii) below), the new method, therefore, has considerable advantages. 

The apparent positions of the intense sources in the constellations of 
Taurus, Virgo, Cygnus and Cassiopeia (I.A.U. Nos. o5N2A, 12Nr1A, 
19N4A and 23N5A respectively) were compared with the corresponding 
positions of the optical objects and were used to derive the collimation error 
of the system in declination. ‘This error 3'-0 (+1':3) was then used to 
correct the observed positions of the other sources. 

(ii) Limitation of reading accuracy, and noise fluctuations—The probable 
errors quoted arise mainly from the noise fluctuations and were derived from 
the scatter of the times of at least eight crossovers, obtained during each 
individual transit observation, from which the time TJ across the fictitious 
meridian was derived. Many of the transit observations were repeated 
several times. 

(iii) Confusion effects—The errors due to the presence of weak sources 
in the reception pattern can be computed statistically from the spatial 
distribution of sources. Some indication of the extent of these confusion 
effects is also obtained from the scatter of the times of crossovers as in (ii) 
above, and these effects are included in the quoted errors of the less intense 
sources. 
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(iv) Ionospheric refraction.—The refraction of the radio waves in the 
terrestrial ionosphere has been investigated by Smith (3), who showed that 
the error in right ascension is a maximum at dawn and dusk, but the error at 
this wave-length is unlikely to exceed 0’-15 of arc. 

The error in declination may be larger but is unlikely to exceed 0'-5 of arc 
for zenith angles less than 60°. 

Since these errors are less than those discussed above no corrections for 
refraction have been applied. 


The positions have been corrected for precession, nutation, and aberration 
and refer to epoch 1950-0, and are listed with the 2C survey (§) numbers. 


TABLE I 
Approximate 
flux 
Number R.A. — Error Dec. Error AR.A. density Notes 
hms 8 : : s 10-*° MKS 
2C 34 0022383 3 63 53 3 16:2 110 IAU ooN6A 
00 22 37 2 63 52 2 (a) SN 1572 
94 0©10612°'7 I°0 13 02 4°3 3-2 58 
or 0613 =3°0 13 04 9 
379 94 15 00°3 1°5 37 51 2°5 4°9 60 
041505 2:0 37 50 4 
404 04 33 53°5 1°5 29 37 3°0 4°1 204 TAU 04N3A 
043355 4:0 29 34 3 
440 O05 0117°4 2°0 38 03 2°8 5°0 55 
05 OI 21 4°0 38 03 8 
724 o81000°9 I°5 48 23 22 7°0 95 IAU o8N4A 
08 10 03 * 3:0 48 22 3 
805 0917 52°1 2:0 45 53°5 2°0 6°4 42 IAU ogN4A 
091754 # 2:0 45 52 2 (d) 
855 09 58 56°3 3:0 29 04 2°8 4°1 30 (6) 
09 58 56 7:0 29 00 10:0 
1175 14 09 32°0 2°0 52 25 3°4 8-4 74 IAU 14Ns5A 
I4 09 32 3:0 52 26 50 
1185 141641°2 2:0 06 46°5 = 4°5 3°2 63 (c) 
141640 8:0 06 46 7:0 
1402 16 2656-2 1°8 39 39°5 _2°0 5'1 49 IAU 16N4A 
16 26 54 5:0 39 38 3:0 (d) NGC 6166 
1432 1648 42°91 1°5 05 04 6°3 311 300 IAU 16NoA 
16 48 43 950 0510 10:0 
1473. 1717596 1°8 —00 57 6°8 31 180 
17 17 59 3:0 —00 50 6-0 
1569 18 28 16°5 3:0 48 42 2°5 76 70 
18 28 12 3:0 48 43 3:0 
1686 201219°8 2:0 23 24 3°6 3°6 102 (b) 
201217 2:0 23 26 5°0 
1775 21 21 30°7 3°0 2449°5 4°2 3°8 62 (c) 
2F OT 91 4:0 24 48 o 
1870 22 43 33°9 3°0 39 19°09 28 51 50 
22 43 30:0 5:0 39 21:0 3:0 


Notes to Table I 
(a) Supernova of 1572 (6). 
(b) 2C position changed by a R.A. lobe. 
(c) 2C position changed by a Dec. lobe. 
(d) IAU designation incorrect. Should be 16N3A. 
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Mr D. O. Edge has kindly made available some preliminary positions 
obtained in the 3C survey by the alternative and largely independent method, 
and these are included in Table I for purposes of comparison; they are given | 
in italics beneath the results of the present observations. His measurements of 
flux density are also included in column 7 in units 10-6 watts m~? (c/s)-". 

Since the axis of the instrument is not exactly east-west, the area defining 
the probable error of a source has the form of a parallelogram, which may be 
derived for each source by reference to Fig. 6. The inclination of the 
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Fic. 6.—The parallelogram defining the probable errors of the radio sources. 





parallelogram may be derived from the value given in column 6 in the table, 
A R.A. which is the change in R.A. for a change of 10’-o of arc in declination. 

Possible identifications both with optical objects and with radio sources 
found in recent surveys are indicated in the final column. 
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providing the positions and flux densities from the 3C survey, prior to 
publication. He also wishes to thank Dr J. R. Shakeshaft for kindly reading 
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these observations could not have been undertaken. 
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RADIAL VELOCITIES OF FUNDAMENTAL SOUTHERN STARS 


William Buscombe and Pamela M. Morris 
(Received 1958 July 3) 


Summary 


Newly determined radial velocities are given for 12 standard-velocity 
stars and 63 other bright stars south of the equator. The plates, taken with 
the three-prism Cassegrain spectrograph of the 30-inch reflector, have a 
dispersion of 36 A/mm at Hy. The Stromlo velocities of the standard- 
velocity stars were found to agree satisfactorily with measurements at other 
observatories on the Lick system. Velocity residuals for individual spectral 
lines are discussed. 

The stars on the programme have been chosen from the FK3 and N3o 
position catalogues. At least three spectrograms of each star have been 
measured. A luminosity class is given on the revised Yerkes system for stars 
of early type. 29 new constant velocities are announced, as well as 18 whose 
variability has been detected for the first time. Observations have been 
resumed on 15 stars previously known to have variable velocity. 





Selection of programme stars.—In 1954 June, when repeated tests on I.A.U. 
standard-velocity stars showed that the performance of the three-prism Cassegrain 
spectrograph was reliable for velocity determinations, regular observing was 
commenced on stars from two principal lists. ‘The first included stars of early 
type in the region of the Scorpio—Centaurus cluster, and the other was planned 
for observations during the spring and summer, when the cluster cannot be 
observed. As the practical limiting magnitude is around 7-0 (photographic), 
stars in the standard catalogues FK3 (1) and N30 (2) were chosen, for which 
accurate proper motions are available, but reliable radial velocities lacking. 
From right ascensions 17 through o to 8 hours, south of declination + 20°, 
about 300 stars brighter than visual magnitude 6-5 were selected, for which either 
no velocity has been published, or else the velocity remains uncertain, i.e. quality 
c, d, or e in the General Catalogue of Stellar Radial Velocities (3). 

Standard velocity stars.—The observing experience of the initial test period 
showed that for an observatory as far south as Mt Stromlo (latitude — 35°) the 
standard-velocity stars recommended to Commission 30 of the I.A.U. (4) are not 
sufficient. About 20 additional stars south of declination —30° were provision- 
ally selected from those with velocities of quality a in the Radial Velocity 
Catalogue, and consistent agreement between Lick and Cape observations. 
The selection was made in such a way that at least one of these stars could be 
observed at any time of the night. Several of these have been observed on a 
significant number of nights with the stars of the N30 programme, and the 
measures are summarized in the first part of Table I. 

Observation and measurement.—Details of the observing routine at the 
telescope, and processing in the dark-room, have been outlined in an earlier 
communication (5) to the Society. To eliminate the curvature of the spectral 
lines on the plates, a curved slit has been in use since 1956 August. All the 
measurements reported in the present paper have been made with Zeiss Abbe 
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comparators. ‘The emission lines of the iron arc comparison source have been 
measured so close to the tips adjacent to the star spectrum that corrections for 
curvature of the lines on the older spectrograms were negligible. 15 to 20 
unblended lines from this source are measured on each plate as a precise scale of 
wavelengths. An interpolation curve is drawn showing the deviations from the 
standard settings computed by the Cornu formula for the prismatic dispersion. 
The derived velocities, corrected to the Sun, are listed in Tables I andII. The 
results are quoted only to the nearest integral kilometre per second, as the 
decimal fractions are considerably below the limiting accuracy for real signifi- 
cance. ‘The various participants in this programme hope to announce further 
results from time to time, whensufficient data are available for other groups of stars. 

Wavelengths for stellar lines.—Wavelengths for stellar absorption lines (denoted 
in Table III by the letter P), which Petrie (6) recommended after his studies 
with the single-prism dispersion of 30 A/mm at Hy, were used as a basis for 
the reductions. For each line the residuals from the mean plate velocities have 
been examined. Among A stars in particular, no standards have been selected 
for which the consistency of the average velocities for all lines on the plate could 
be verified. ‘The mean velocity residuals for stars of each range of the later 
spectral types are of the same order as those reported by Wright (7) and Evans 
(8). Measurements of a particular line were considered to deviate significantly 
if the algebraic mean residual AV, of its Doppler velocity from the mean 
velocity V for all lines on a plate, in a group of plates of the same spectral type, 
exceeds three times the standard error o,y of the mean residual. ‘The wave- 
lengths of lines which consistently give velocities in closer agreement with the 
simple mean for a plate are listed in the middle column of Table III. While 
a few of the lines recommended by Petrie give slight systematic deviations 
from the average plate velocities, some additional lines seem consistently reliable 
on the Stromlo spectrograms. 

TABLE | 
Stars with constant velocities 

Column 

I Number in Henry Draper Catalogue (Harvard Annals, 91-99). 

2 Constellation designation. For southern stars lacking a Greek Bayer letter, the 
Gould number is quoted from Uranometria Argentina (Cordoba Resultados, 1). 
Visual apparent magnitude, from the revised Harvard Photometry (H.A., 50). 


3 

4 Spectral type, from the Henry Draper Catalogue. 

5 Mean radial velocity relative to the Sun, from Mt Stromlo spectrograms, to 
nearest integral kilometre per second. 

6 Probable error of the mean, from deviations of individual plate velocities, in 


km/sec. 

7 Number of plates measured. 

8 Revised spectral type on MK system, determined from Mt Stromlo spectrograms 
by M. L. Woods (Mem. Commonwealth Obs., 3 (12), 1955), A. de Vaucouleurs 
(M.N., 117, 449, 1957) or P. M. Morris (unpublished). 

9 Mean radial velocity from R. E. Wilson’s General Catalogue of Radial Velocities 
(ref. 3). The lower-case letters denote quality by internal consistency of the 
published velocity measures from individual plates, as follows :— 

a Probable error of mean less than 0-9 km/sec from at least 5 plates 


b ” ” ” ” ” ” 2°0 ” ” ” ” 3 plates 
€ ” ” ” ” ” ” 5°0 ” ” ”» ” 2 plates 
d ” ” ” ” ” ” 10°0 ” ” ” ” I plate 


e Large probable errors for single plates of early-type stars. 
The succeeding digit indicates the number of plates measured previously at the 
Lick Observatory (designated by the capital letter L). 
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TABLE I (continued) 





Column 
10 Notes, including the suspicion of variable velocity when the probable error of the 
mean velocity, multiplied by \/n, exceeds the average probable error for the 
scatter of velocity measures from individual lines on the various plates. 








(1) (2) (3) (4) () ©) (7) (8) (9) (10) 
HD Star my Sp. V p.e. n_ Sp. Type GCRV Notes 
km/s + 
4128 BCet 2:24 Ko +13 o2 25 G8III +13 
20794 82Eri 4:30 Gs +89 o2 7 G5V +87 
36079 BLep 2:96 Go —14 05 g G2il —14 
81797 xHya 2:16 K2 -—- 3 o4 6 K3 Ill — 4 
108903 y'Cru 161 M2 +22 of 14 M3II +21 
109379 BCrv 2°84 Gs - 8 03 5 GslIl -— 8 
150798 xTrA 188 K2 -— 3 o3 14 Ke2iil —4 
168454 d5Sgr 2°84 Ko -—-—20 o2 19 Keiil —20 
196171 aInd 3:21 Ko —1 06 5 Kolll — 1 
203608 y Pav 4:30 F8 -—30 06 3 Fs V —30 
204867 B Aqr 3:07 Go +6 06 7 Golb + 7 
223647 y' Oct 5-10 Gg +12 0°5 4 G5 Ill +14 
4211 A*Scl 5°97, Ko +24 06 4 
6192 w Phe 6:00 Ko +13. o'5 3 
8498 109G Scl 5°82 K5 -—16 06 c 
9896 12GHyi 6:12 F2 —-19 0%3 3 F2V 
10042 14GHyi 6:06 G5 —1 oF 5 
13940 135G Phe 5°86 Ko +14 1°7 4 Var? 
15975 A' For 5°88 Ko +13. 0°5 4 
22789 7 For 6:08 Ao +39 #1°8 4 AoV 
24626 t Eri 5°12 Bs +25 1°0 5 B6olVv +18c7L 
29435 253G Eri 6:29 B8 +15 1°7 5 BgoIV-V 
| 30422 268G Eri 5-97 A2 +19 1°73. 5 A3ZIII-IV 
: 34172 €Mens°85 Ko —- 5 06 3 
| 34868 12GCol 5:75 Ao +18 1°71 5 AolV 
40953 « Mens5:56 B8 +10 2-4 5 Bos V 
| 51557 tVol 5°52 B8& +19 I°9 4 BO6IV +18¢5L 
183806 63G Tel 5°87. Ao — 8 17 4 Ap 
188097 76G Pav 5°82 A3z —10 o'9 5 Am 
| 194215 296G Sgr 5°97 Ko —14 13 5 Var? 
199623 8GInd 5°88 F5 -—21 0o'5 4 F6V 
| 200073 43G Mic 5:94 Ko +7 I°5 3 Suspected 
subdwarf 
| 203949 68G Mic 5°69 Ko —76 1°3 5 
| 207229 54G Ind 565 Ko —-6 03 2 —id1L 
208741 66G Oct 5:91 F2 +7 = =1°5 3 «+F3ZIII Var? 
| 216437 pind 6:14 Go — 3 06 4 
219077 23G Tuc 6:24 Go —26 14 3 
219765 7Oct 5°56 Ko +31 o9 5 
220096 11GScl 5°81 G5 +13 =#I'! 5 
221420 83G Oct 5°78 Ko +26 08 3 
222820 36G Tuc 5°66 K5 —13 #1°5 3 
223352 S5Scl 4°64 Ao +6 2°1 5 AoV +14d5L 
| 224022 27GPhe 6:01 F8 -—10 0o'5 4 F8&8V 
224361 42G Tuc 6:04 A2 ~—13 2°7 3 A2p 
| 224834 7Phe 5°66 Ko +8 0% 4 
42 
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Tas_e IT 


Stars with variable velocity 


Henry Draper number and constellation designation 
Visual magnitude and HD spectral type 

Julian date and decimal for time of mid-exposure. 
Radial velocity relative to the Sun. 

Internal probable error of deviations from the mean velocity for the plate to 


velocities for the individual lines. 


onl a 


(1) 
HD 
Star 
2884 
B’ Tuc 


4293 
70G Phe 


7312 
102G Scl 


10800 
3G Oct 


15233 
A Hor 


(2) 

my 

Sp 
4°52 
Bog 


6-00 
AS 


5°88 
Go 


5°47 
F2 


Revised spectral type 

Mean radial velocity in Catalogue 
Notes, including references to earlier publications from other observatories 
quoting individual plate velocities. 


(3) 

JD 
7 
5025072 
5315'277 
6036°276 
6094°157 


4966-211 
4968-229 
5325°262 
5716-175 
6035°270 
6120°093 
6144°016 
6150°970 
6154969 


4966°244 
4968°258 
5000°188 
5347°206 
5362-164 
5745 °132 
5760-073 
5776:022 
5812-954 
5825-938 
6036°295 
6067-245 
6086-181 
6120°117 
6125-140 
6151-000 


5357°203 
6151°034 
6182-974 


5001 °206 
5325 °334 
5711°249 


(4) 
V 
km/s 
+ 33 
+12 
+ 36 
+ 9 


+ 6 
+26 
+62 
+26 
+40 
-17 
—40 
—27 
— 36 


+52 
+39 
— 2 
+31 
+20 


— 24 
ag 
— 34 


(5) 
p.e. 


YNWNWNWWPRWUPDNHDNAD 
Me OWUSU AN HWW 6 CWbMNAI 


ON w 
COW a 


mw nN 
NW N 


Vol. 118 


i As in Table I. 


As in Table I, cols. 8-9. 


(6) (7) 
Sp GCRV 


B8 V +10c7L 


A7ZIII + 1c5C 


A7 Ill 


F2IIp +27ce5L 


(8) 


Notes 


Lick Publ., 16. 


M.N., 117. 
Under observation. 
for orbit 


Under observation 
for orbit 


Lick Publ., 16. 














— 











we 














—_—-— 


en ee 








XUM 


No. 6, 1958 
(1) (2) 
HD my, 
Star Sp 
20121 5°92 
71GEri Fz 
21722 5°96 
59G Hyi Fa 
21981 6:o1 
44G Hor Ao 
22233 5°60 
45G Hor Ko 
24706 5°77 


55G Hor Ko 


41843 5°72 
74G Col Ao 
46936 = 5°57 


50G CMa Bg 


54732 598 
28G Car G5 


59635 54! 
98G Pup B8 
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(3) 
JD 


a 


5414‘111 
6125°173 
6151°067 
6176-010 
6178-005 


5357°271 
6115°214 
6125208 
6151°116 
6183:017 


5357°283 
5414130 
5741°231 
5760°186 
6094°236 


5385°217 
5742°247 
6095 °238 
6126-218 
6148-120 
6177°063 


5348°317 
5361°254 
6128-205 


6154°133 © 


6156-082 
6212-969 
6216-010 


5025 °269 
5800-167 
6150°219 
6151°206 
6176°132 


5778°256 
5809°146 
6176°157 
6178°134 
6262°915 


5062°256 
6176°197 
6262°949 


5826°171 
6176°227 
6178-199 
6262-976 


TABLE II (continued) 


(4) 
V 
km/s 


| 
NW ON 0ODwW 


++ 
N 
ioe} 


i 


(5) 


p.e. 
a 


I°9 


(6) 
Sp 


F6 III 


Fs IV 


B8 V 


B3 III 


(7) 
GCRV 


+34¢4L 


+40c2L 


+42¢3L 


+26c4L 


(8) 


Notes 


Lick Bull., 15. 
Visual binary 


Lick Bull., 15. 


Lick Publ., 16. 


Lick Publ., 16. 
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(1) (2) 
HD my 
Star Sp 
67523 2°88 
p Pup Fs 
173117 5°76 
94G Sgr B8 
175510 5°03 
A Tel Bog 
188899 5°05 
61 Sgr Ao 
189198 5°80 
259G Ser As 
193571 = 5°64 
x’ Sgr Ao 
199443 5°95 
64G Cap A3 


William Buscombe and Pamela M. Morris 


(3) 
JD 
re 
5585880 
5585 °897 
5585928 
5585 °954 
5585976 
5586-001 
5586-022 
5778°170 
5778°185 
5778°205 
5778°222 
5778°244 
5778°267 


5612212 
5985°171 
5987°173 
5997°167 


5270°166 
5360-919 
5586°274 


5265°160 
5293°135 
5647°133 
6036-084 
6083 °967 


4967-046 
5591 °328 
5612°280 
5646-194 
5715°990 
5995 °202 
6035°118 
6038°097 
6067-610 
6083-987 
6085-946 
6092-006 


5000°937 
5298°126 
5734°015 
5741°929 


4974°029 
4999°961 
5028°917 
5290°199 
5293°179 
5987-262 


TABLE II (continued) 


(4) (5) 
V p.e. 
km/s + 
+40°5 0°9 
+ 39°2 0°9 
+4a°7 tt 
+46°3 0°7 
+48°7 1°2 
+42°0 I°0 
+ 44°3 1°0 
+47°3 1°2 
F4r°7 1°7 
+42°5 1°4 
+45°'0 I'l 
+47°8 1°6 
+52°5 1°9 
+34 2°0 
—- 6 69 
+48 4:0 
+15 3°6 
+10 4°6 
+18 7:8 
—-2 36 
+2 3°3 
+64 6°1 
—-2 2°'5 
+i2 3° 

o «6:2 
+37 3% 
+34 «4/1 
+24 3°2 
+13 58 
— 2s #3 
—23 17 
—-12 16 
+20 2°5 
+15 2°6 
—28 2°5 
+16 2:8 
—16 2°0 
+ tf a2 
—14 24 
—28 2:9 
+7 332 
—-15 2°5 
—20 34 
+24 3°4 
+25 3°9 
+16 3:0 


(6) 
Sp 


Bs V 


Bg III 


Az IV 


A7 IlIp 


A1V 


As IIIp 


(7) 
GCRV 


+ 47a 


— 6d5L 


— 4c8LY 


= 6c4C 


—17¢3L 


Vol. 118 


(3) 
Notes 


Lick Publ., 16. 
Cape Ann., 10(8). 
Ap. f., 124. 
Observatory, 77. 


Lick Publ., 16. 


Lick Publ., 16. 
Yerkes Publ., 7. 


M.N., 117. 
Under observation 
for orbit. 


Lick Bull., 15. 
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TABLE II (continued) 


(1) (2) (3) (4) (5) (6) (7) (8) 
HD my JD Vs pe. Sp GCRV Notes 
Star Sp 243.... km/s + 


203705 5°54 4974078 3+ 7 3°1 Ag Ill 
18 Aqr As 4999°984 + 4 2°0 
5024°947 —19 4:2 
5354012, +39 3°! 





204960 5°73 5266:275 +16 =I'1 
3G Gru Ko 5354025 — 1 1°2 
5716057 +6 13 

6035°149 +22 I1°2 


207052 5°43 5379°968 —24 2°'5 Az V + 1c6LY = Lick Publ., 16. 


ACap Ao 5742'°985 —23 2°8 Yerkes Publ., 7. 
6036149 +33 2°9 
6067°083 — 6 2:7 


211539 5°74 5057°937 +32 3°2 
vOct Ko 5348-092 +21 1°8 
5354°072 + 5 7 


—————— Oe ee 


211838 5°36 5025°970 —10 O6°'1 B8 V — gc7LY Lick Publ., 16. 
pAgqr B8 ~— 5360-053 Oo 4:2 Ap. F., 64. 
6067‘101 +30 5°8 
6086:035 + 8 68 





212132 5°82 4967:128 +43 #1°7 Fo IVp 
a Gru Fo 5269°283. +32 ° 
5759°967 —14 2°1 


214085 6-11 5711119 +10 2°58 A4 V Under observation 
51G Gru A2 5743°025 —1I 2°%4 for orbit. 
5771°958 —36 2:1 
5777°942 —18 3:2 
6038:149 +52 253 
6067'125 +31 2:8 
6086-061 —31 2:2 


GG a, SE 


215789 3°69 §=64971°133. «+11 1°9 A2V oc3L Lick Publ., 16. 
e Gru A2 5321°178 — 3. 2°91 
6036:203 + 7 2°0 
6084:071 — I 3°3 
6094°076 —23 273 
6115°009 — 6 1°74 
222095 4°86 4970°181 +26 3:5 A2V + 1ocsL Lick Publ., 16. 
11G Phe A2 5295'207 +51 3°6 Under observation 
5321°228 — I 14 for orbit 


5716:092 +31 2°5 
5749°050,  —32 3°4 





5987°318 +9 21 
6115°043 — 4 46 
6143°9388 + 1 3:8 
6154°049 — 7 I°5 
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William Buscombe and Pamela M. Morris 


(1) (2) (3) 


HD my, 
Star Sp 


222847 
106 Aqr B8 


225132 
2 Cet Bo 


Sp. 
Type 


Bo—Bg 


Ao-F2 


F4-Gs5 


G8&-Ks 


5°26 


4°62 


JD 
a 
4967°186 
5269°347 
5361°128 
6036°237 


5295 °347 
5361°152 
6035°255 
6115°065 


TABLE II (continued) 


(4) 
V 
km/s 
+48 
+ 32 
+23 
+17 


+ 33 
+ 33 
+ 33 
+26 


(5) 
p.e 


(6) 
Sp 


B8 V 


Ao :IV 


TABLE III 


(7) 
GCRV 


+14d6LY 


—5d13LY 


Wave-lengths of lines consistently measured 


Velocity 
Low 


3933 °664 P 
3964°727 P 
3968-465 

3970°074 P 
4130°876 P 


3933°684 P 
3968-404 P 
4118-678 
4202°042 
4325°765 


4077°710 P 
4132°230 

4282°640 P 
4325°777 P 


4005°205 P 
4034°494 
4101°743 
4132°230 
4282-640 P 


4009 °27 
4026'14 P 
4101°737 P 
4128°051 P 
4143°759 


4005 °205 

4030°615 P 
4045°796 P 
4063°555 P 
4071°687 P 
4077°710 

4101°750 P 


4005°205 P 
4030°673 P 
4045 °827 P 
4063-607 P 
4071°751 P 
4092°710 P 
4101°743 P 
4118-702 


4030°673 P 
4045 °827 P 
4063607 

4071°751 P 
4077°710 P 
4092°512 P 
4118-702 P 
4143°772 

4202°093 P 


Velocity 
Average 


4267°167 P 
4340°468 P 
4471°477 P 
4481°228 P 


4132060 
4143 °686 
4260°429 
4307°914 P 
4404°726 P 
4415°137 P 


4143°772 
4202°093 P 
4215°761 
4260-429 P 
4307°914 
4340°540 P 
4404°745 P 
4415°137 


4235°913 
4254°330 P 
4260-429 P 
4325°777 P 
4404°745 P 
4415137 P 
4427°319 P 
4454°929 

4461°793 


Vol. 118 


(8) 
Notes 


Lick Publ., 16. 
Ap. f., 64. 


Lick Publ., 16. 
Yerkes Publ., 7. 


Velocity 
High 


4069°766 P 
4072°162 
4075868 
4116°103 
4387-928 P 


4215°600 P 
4250°510 P 
4254°348 P 
4340°526 P 
4481°293 P 


4134°478 
4191°593 
4235 °913 


4134°478 
4191°593 
4215°761 
4307°914 
4340°540 P 
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‘TaB_Le IV 


Observers and measurers 








Name Standard Velocity N30 Programme 
Observed Measured Observed Measured 
W. Buscombe 41 65 142 183 
H. Gollnow 68 10 29 13 
J. Graham * I 13 I 22 
G. Hagemann 5 87 2 13 
W. Heintz t 15 I 35 ° 
P. Morris 6 17 41 125 
A. Przybylski 42 15 63 ° 
Totals 178 313 


* Sydney University: summer assistant, 1957-58. 
+ Munich Observatory: guest observer, 1954-55. 


Acknowledgments.—In presenting the results of a part of our extensive 
programme for determining radial velocities, we wish to thank the other observers 
and assistants for their help and co-operation. ‘The work could never have got 
under way without the painstaking attention by Dr Gollnow to many details 
of the final design and operation of the spectrograph, and his collaboration in 
choosing the observing lists of programme and standard stars. Messrs Banham 
and Thomas, with their associated craftsmen, have spared no effort in applying 
their skill to the construction and maintenance of the instrument. 


Mount Stromlo Observatory, 
Canberra, Australia: 
1958 May 15. 
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SPECTRAL TYPES AND RADIAL VELOCITIES IN THE 
GALACTIC CLUSTER NGC 3293 


M. W. Feast 


(Communicated by the Radcliffe Observer) 
(Received 1958 May 7) 


Summary 


112 spectra of 35 stars have been used for the determination of spectral 
types and 89 of these spectra for the determination of the radial velocities 
of 26 of the stars. The spectral types taken together with Dr T. E. Houck’s 
photometric measures show the H—-R diagram expected for a young cluster. 
The stars are mainly of early B type with luminosity classes from V to Ib. 
An Mo fab star is almost certainly a cluster member. The radial velocities 
show a dependence on magnitude (Fig. 3) which is large and unexplained. 
The Ca 11 interstellar lines are free of the effect, which may be possibly related 
to the Trumpler effect. 





1. Introduction.—The study of galactic clusters containing early B type stars 
is of interest for several reasons. It appears to be generally agreed that such 
clusters must be quite young and that a study of their Hertzsprung—Russell 
diagrams is of importance for theories of stellar evolution. In addition, the 
H-R diagrams are of interest in checking the luminosity scale of the MK system. 
Again, it would appear rather likely that a considerable number of these clusters 
are dynamically unstable and will expand eventually, passing through a stage in 
which they would be identified as stellar associations (cf. (1) etc.). A study of 
the motions of the cluster members is therefore desirable. Finally, a 
determination of the radial velocities of clusters of early B type stars is of 
importance for galactic rotation studies. 

NGC 3293 is a compact cluster in Carina (10"34-5-58° 00 (1950)) not far 
from, and perhaps associated with, the 7 Carinae complex (see for instance the 
photograph shown in (2)). Fig. 1 (Plate g) is reproduced from a photograph of 
the region of NGC 3293 taken by Dr A. D. Thackeray at the Newtonian focus of 
the Radcliffe 74-inch reflector. The cluster is also shown on the charts of the 
HD extension (3) (chart No. 127). A few stars Npr. the main cluster and 
which may be associated with it are included in the present discussion. Star 
No. 1, the only O type star in this region, is the centre of an emission nebula 
(Gum (4), No. 30) and the whole cluster is involved in faint nebulosity 
(c.f. (2), Fig. 6). Mrs E. v. P. Smith (2) has published polarization measure- 
ments of seven stars in the cluster. Miss D. Hoffleit (5) and Morgan, Whitford 
and Code (6) have published the spectral types of a few stars in the cluster. 
Objective prism surveys have revealed 10 OB stars in the cluster (7) and one 
Mo supergiant (8). The spectral types and radial velocities of the two brightest 
stars were included in the first Radcliffe catalogue (9) and of No. 1 in the second 
catalogue (10). Some years ago an important, unpublished, photometric study 
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of the cluster was made by Dr T. E. Houck at the Cape, Boyden and Radcliffe 
Observatories. NGC 3293, which is listed as Lacaille 4375 in the Cordoba 
work (x1), has been referred to as the II Carina association (12) though it is 
considerably smaller than the well-known O and B associations. 

2. Spectral types.—The present study of NGC 3293 is based on 112 spectra 
of 35 stars. The spectra were all obtained with the two-prism spectrograph at 
the Cassegrain focus of the 74-inch Radcliffe reflector. 86 spectra were obtained 
with the {/3°7 camera (49 A/mm at Hy), 3 with the f/6 camera (29 A/mm at Hy) 
and 23 with the f/2 camera (86 A/mm at Hy). In estimating the spectral types, 
all these plates were employed but the radial velocities dealt with later depend 
only on the spectra taken with the two higher dispersions. ‘The types were 
estimated in the usual way by comparison with standard stars of the MK system. 
‘Though types of the three stars, Nos. 1, 3 and 4, have been published previously 
(9) (10), the types given here (which differ very slightly from those given earlier) 
are independent re-determinations. ‘The 35 stars studied include the 31 stars 
for which Dr Houck has made photometric measurements. 

The results are given in Table 1. The HD, HDE or CPD number is given 
where possible and the stars are identified in Fig. 1. ‘The types of Hoffleit (5), 
Morgan, Whitford and Code (6) and the HD types are included where available 
for comparison. The fragmentary overlap with Hoffleit and with Morgan, 
Whitford and Code may be regarded as satisfactory. Columns 7 and 8 contain 
data kindly supplied in advance of publication by Dr Houck. The visual 
magnitudes, my, should be close to the Johnson V magnitude system. The 
colours (B-V) were converted from Dr Houck’s instrumental system using a 
relation supplied by him. Using the intrinsic colours of the O and B type stars 
given by Hiltner (13), the colour excesses and absorption corrections were 
determined (a ratio of total to selective absorption of 3-0 being adopted) and 
the magnitude corrected for absorption, m,, found. The absorption corrections 
range from o™-57 to 1™38 and an assumption of uniform reddening over the 
whole cluster would have resulted in a very poorly defined sequence in the HR 
diagram. Column g of Table I gives the values of m), and column to the 
distance moduli of the individual stars using the calibration of the MK system 
given by Johnson and Hiltner (14). 

Star No. 21 is that previously listed as an Mo supergiant by Blanco and 
Miinch (8). ‘The slit spectra now obtained confirm and refine this classification. 
The intrinsic colour of an Mo Iab star was taken from Johnson and Hiltner (14). 
The distance moduli of the B stars and of No. 21 agree sufficiently well to confirm 
the membership of the star in the cluster. The radial velocity results of the 
next section also support this result. Stars 1, 2 and 11 inthe Npr. “‘flare’’ do 
not have significantly different distance moduli from the cluster stars and were 
included with them in deducing the mean modulus. It cannot be claimed with 
any certainty that these stars do belong to the cluster though a “‘ flare ’’ of this 
kind would not be unique (cf. NGC 2169 (15)). However the mean distance 
modulus of the cluster is not appreciably affected by including or omitting these 
stars. A mean modulus of 12-08 is found from 31 stars. A standard error, due to 
scatter of the individual moduli, of +o™-1 is found. This, of course, does not 
include the uncertainty due to systematic errors in the calibration of the MK 
system, etc. This distance modulus corresponds to a distance of 2-6kpce with 


an ‘‘ internal ’’ standard error of about + 0-1 kpc. 
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Star 


Mab WN 


“ 
oo ON SO 


31 
32 
33 
34 
35 





HD HDE 
or 

CPD No. 
91824 
91850 
91943 
91969 
92024 


92007 
91983 
92044 
303065 
303067 


303068 
303075 
—57° 3522 
—57 3507 
—57 3504 


— 57° 3500 
—57° 3514 
3524 
3524 
3523 


7 3502 
7 3506 
7 3506 
7° 3517 
7° 3521 


—57 3526 
> 3527 
—57 3531 
° 3518 


3528 
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TABLE I 


Spectral types and photometric data 


Hoffleit M.W.C. 

Radcliffe type type 
MK type (5) (6) 

O7 O7 
Br Ill Bi IV: 
Bo:s Ib Bo:s I 
Bo I(b?) Bola 
Br Ill 


Bos Ib 
Bo Ib 


Bo Il 
Bx Ill 
Bos II 


Bo-s III 
Br Ill 
Bo:s III 
B2 V 
Br V 


B1 V 
Bi V 
B2 V 
Br V Bo-5 V 
Bi V- 


Br Ill Bo-s5 V 
B2-3 V: 

Bi V 

Bx Ill 

Br Ill 


Mo Iab(b ?) 
Br II 
B: III 
Br IIT 
Br III 


Bx II 


Br Ill 
Boss III 
Bri V 
B2 V 
Bs V:: 


HD 
type 
Oe5 
B2 
Bo 
Bo 


My 
815 
g°12 
6°72 
6°52 
9°06 


8-22 
8°59 
8-27 
9°99 
9°56 


10°02 
10°03 
10°18 

9°28 
10°OI 


8°74 


9°30 
g‘1Io0 
8:03 


7°31 
7°62 
9°21 
9°22 
8°15 


9°28 
8°95 
10°25 
10°67 


11°55 
12°97 
10°65 


B-V 
—o0-05 
+0°14 
+0-06 
+o°o!I 
— 0°04 


+0:07 
+0°04 
+O°15 
— 0°02 
+ 0°04 


—6°OL 
+ 0°05 
—0°03 
— 0°04 
+0°'o! 


—0'02 


— 0°07 
+018 
+ 0°02 


+ 2:06 
+0:'07 
+0:°03 
+ 0°03 
+0:°06 


—0'04 
+0°05 
+0°04 
+0:02 


0°00 
+ 0°23 
+0°04 


Vol. 118 
My 2%9-M 
7°34 =12°4 
7°92 123 
5°79 «11'S 
5°59 11°8 
8-40 12°8 
7°59 12°5 
7°69 «12°! 
6-98 11°7 
9°33 «I1'3 
8-66 11°8 
9°27 12"4 
g'IOo. =12°2 
9°55 11°6 
8-62 I1°7 
9°20 12°3 
8:02 12°4 
93 11S 
9-92 2h 
719 «11°6 
6°17 I1'5 
6°63 «11°5 
8°34 12°7 
8°35 12°8 
719 «11°6 
8:62 13°0 
7°96 12°7 
9°35 12°4 
9°89 «11°9 
10°83 _-12°8 
I2°OI I1°5 
9°81 11°8 


Note.—The following pairs of stars have the same CPD numbers : 
6 and 26, 18 and 19, 22 and 23. 


The HR diagram for NGC 3293 is shown in Fig. 2. The luminosity class 
assigned to each star is indicated in this figure and it will be seen that the 
correlation of luminosity class and my (or M) is satisfactory. The scatter in m, 
at a given type for each luminosity class is no greater than that found in other 
works of this kind (e.g. h and x Per (16)). It is not clear whether this scatter 
could be reduced by refining the process of spectral typing or whether most of it 
is inherent in the method of spectroscopic parallaxes. 


There appears to be 
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no tendency within the accuracy of the results, for different subgroups (of 
type and luminosity class) to give systematically different distance moduli. That 
is, there appears to be no significant deviations from the Johnson and Hiltner (14) 
calibration of the MK system. It should be pointed out that we have not taken 
into account the fact that a good proportion of the stars are probably binaries 
(see Section 3). This will, of course, affect the magnitudes to some extent. 
However, the stars known to be variable in velocity do not appear to differ 
systematically in the HR diagram from the other stars. 


SPECTRAL TYPE 
B3 BS 








O7 BO BO.S_ BI B. 
$.0 T T T T T T 8S B7 8 “wo 7~7:0 
@Ie 
ele 
6-OF 
Tab @ 4-6-0 
7-0 4-s.o 
™Mo M 
8.0 4-4.0 
9.0 4-3.0 
10.0 42.0 
1.0 1-0 
12.0 0.0 
13.0 4 1 1 1 1 1 i 1 1 4 











i i L. 
=.24 -.20 =216 ~.12 -.08 + 1.68 


INTRINSIC COLOUR 


1 
—-32 --28 


Fic. 2.—Hertzsprung—Russell diagram for NGC 3293. The absolute magnitude (M) scale 
depends on a distance modulus of 12:1. The luminosity classes of the stars are shown. The line 
is the ‘ zero age’ main seq e of Foh and Hiltner (14). 





The sequence defined by the B stars in the HR diagram is like that found 
in other groups of similar stars, e.g. the Orion region (17) and h and x Per (16). 
Its interpretation as evolution away from the main sequence for the more massive 
stars is well-known. The one M type supergiant found is similar in type to the 
M stars found in h and x Per and similarly fits into the present theories of stellar 
evolution. 

The O7 star (No. 1) in the “‘ flare ’’ does not follow the general trend but 
probably lies on or near the main sequence. This fact, however, cannot be used 
as evidence against its membership in the cluster. In fact O type stars which 
deviate from the general sequence in this way are found in several clusters and 
associations, e.g. Orion Nebula (17), h and x Per (16) and I Gem (18). 

It is unfortunate that magnitudes and colours are not available for four faint 
stars in Table I, as these would have strengthened the lower part of the HR 
diagram. Extension of the photometry to fainter stars is clearly desirable in 
view of the possibility of finding deviations from the main sequence, e.g. (19). 
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However, accurate spectral types for many fainter stars might present some 
difficulties, as the faintest stars in the present programme reach about 13™(my). 
If we adopt a distance of 2-6kpc and a radius of about 4-5 minutes of arc, then 
we obtain a linear radius of 3°5 pc, a fairly average value for an open cluster. 
For convenience, a number of stars in the field of Fig. 1, which have not been 
included in the present work but for which HD or HDE types are available, 
are collected together in Table II. Few, if any, are likely to be cluster members. 


TaBLe IT 
Field stars (probably non-members) 

No. HD or HDE Type 
36 303066 A2 
37 303069 B8 
38 303070 A 
39 303071 A 
40 303072 G5 
41 303073 B8 
42 303074 Ao 
43 303076 Ks 
44 303061 B 
45 303062 B8 
46 303060 Ao 
47 302981 Ao 
48 92121 A3 
49 303077 Fo 
50 303064 B3 
51 303063 Ao 


3. The radial velocities of the stars.—The radial velocity data are collected 
in ‘Table III. In the case of star No. 3 the first six observations, and in the case 
of star No. 4 the first five observations are taken from (g) whilst in the case of 
No. 1 the first four observations are taken from (10). The observing and 
reduction procedures employed were those previously described (9, 10). It is 
noted that, since all calculations were originally carried out to the nearest 
o-1 km/sec and only rounded off in compiling Table III, the quoted means do 
not always agree exactly with the means of the tabulated velocities. Of the 
26 O and B stars in the programme, 9g, or one star in three, give evidence of 
being binaries, either from showing variable velocities or doubling of lines. 
It is possible that one or two of these stars might, on further examination, turn 
out to have sensibly constant velocities. However, since we detect binaries only 
when their orbits are suitably inclined to the plane of the sky, we may in any 
case conclude that NGC 3293 is rich in binaries. 

The mean velocity of the cluster is — 14-1 +1-8(s.e.)km/sec. This was 
determined from those B stars for which two or more spectra were available and 
which were not classified as variable velocity stars. The three stars in the flare 
(Nos. 1, 2 and 11) were omitted from the mean, which therefore depends on 
15 B stars. It will be seen that this mean agrees very well with the Mo super- 
giant (—12+1(s.e.) km/sec) supporting the latter’s cluster membership. The 
velocity of the Mo star was determined using the wave-lengths recommended 
for the late type stars by the Victoria workers (20). Though these wave-lengths 
were not tested for stars later than K8, their use for an Mo star will not have 
introduced any appreciable errors. 
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B TaBLe III 
‘ Radial velocities in NGC 3293 
4 Star Call Star Call 
km/sec km/sec km/sec km/sec 
No.1 O7 No. 7. Br III 
1952 IV. 3°801 = +27 (8) + 7 (1) 1956 IV. 11-798 —36 (8) + 4(1) 
IV 5°750 -—22 (4) + 7) 1957 IV. 13°828 +48 (9) — 1 (2) 
IV 6792 +7 (4) VI 9802 —18 (9) —11:(1) 
Voo715 39-22 (7) «2+ 5 (1) —_ =" 
1957 IV 20°867 —36 (8) + 1 (1) ' ’ : 
IV 23794 —4s (8) + 8(2) Variable velocity. Range 84. 
VI 9783. —14 (6) + 7 (2) 
— 15: + 6 No. 8. Bors III 
Variable velocity. Range 72. 1956 IV. 11°820. —15(15) + 4(1) 
1957 IV 137753 —12(15) — 3 (1) 
No. 2. Br III VI 9°768 = (9) -17 (1) 
1957 IV 20°897 —16 (7) — 2(1) cae 5 
IV 23°817 — 4 (5) o (2) 
VI 12-742 — 7 (5) «+11 (2) 
—*9 +3 No. 9. B2V 
: 1957 IV 20°979 «+ 4 (8) +16 (2) 
No. 3. Bos Ib 1958 III 2°872 = +‘ 5: (6) o (1) 
1952 IV 9 ‘829 o (8) +11 (1) III g:097 ~—24 (9) +10 (1) 
IV 16-838 —21 (15) 0 (1) —5 +9 
V 2827 -— 6(14) + 4(1) : : 
1953 I 19069 — 9 (8) + 1(2) Pag variable rcomig: a 
I 24'029 — 7 (18) — 3(1) pectra show some signs of line doubling. 
IV 9850 — 2 (14) + 3 (1) 
1958 IIT = 37125 93=-— 7 (12) «2+ 7 (1) No. 10. B1 V 
Ht St 6 - 300) 06+ 1G) 1957 II] 22-791 + 6 (7) + 1(2) 
“—— = IV 23804 —20 (8) — 5 (1) 
No. 4. Bo I(b ?) 1958 III 9002 —18 (9) — 2(1) 
1953 I 22°071 — 1 (11) + 7 (1) ——e ~ 
I 27056 + 4(17) — 4(1) 
I 31°051 o(15) — 1 (1) Ne: 11.. Bx V 
III 28-843 o(10) + 1 (2) 1957 IV. 23°772 — 7(10) — 2(1) 
Vi 24°750 — 4(10) + 6(2) 
1958 IIT = 3°139 )S — 5 (11) + 9 (2) 
II] 8831 — 8 (8) +121() No. 12. BrV 
— + 3 1957 IV 22-940 —56: (3) — 5 (2) 
Lines broad and indications of doubling. 
No. 5. Bz Ill 
1956 IV. 18-788 +18 (8) — 3 (1) No. 13. B2V 
1957 IV 23856 —33 (8) — 2 (2) 
, “ 1957 IV 22975 — 6(10) — 2 (1) 
VI 12-781 34 (10) +14 (2) 1958 III —- o (8) + 5(1) 
— 16: + 3 =" i ie 
Variable velocity. Range 52. 
No. 6. Boss III No. 14. B1 V 
| 1956 IV 11°774 25 (7) «+ 3 (2) 1957 II] 22°814 —16 (6) +12 (1) 
1957 IV. 13°843 -— 9 (6) — 1(2) V 28-771 —28 (9) + 5) 
| VI 9730 «6-28 (5) 3+ 10 (1) VI 9689 — 4 (8) +12(1) 
—20 +4 VI 12822 — 9 (7) +20:(1) 
Lines rather poor. —14 +11 
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TaBLe III (continued) 
Star Call Star Call 
km/sec km/sec km/sec km/sec 
No. 15. BrV: No. 23. Br III 

1957 IV 23°:701 —10:(5) + 1 (2) 1957 IV 20°837. — 1(11) — 9 (1) 

1958 III 3:054 : — g(1) V 28694 +19(12) + 2 (1) 

— 4 VI 15°775 -—10 (9) — 4(1) 

Lines scatter widely and there are + 3: — 4 


indications of a binary nature. No 
satisfactory mean can be determined for 
the second spectrum. 


Possibly variable velocity. Range 29. 


No. 24. Br III 
No. 16. Br III 1957 IV 7880 -— 5 (6) — 1 (1) 
1956 IV 18814 —26(13) + 1 (1) m4 -— a se +a 
1957 IV. 7-909 17 (17) — 1 (1) 5 “ 
VI 9829) «©6027 (12) —12 0 
—23 0 
No. 25. Bz III 
me 6. OV 1956 IV 18-754 —22 (6) — 4(1) 
1957 IV 7815 3— 4(11) + 1 (2) 1957 IV 20°792 —14 (7) — 2 (2) 
V 28808 —23(11) + 4(1) VI 9°744 —24 (5) — 1 (1) 
VI 15°731 —22 (7) — 6(1) _ — 
— 16 a Lines rather diffuse. 
No. 19. Br III : 
1957 IV 7850 —18 (5) + 2 (2) No. 26, Br III 
V 28-849 —-19 (4) + 2(2) 1956 IV 18-838 + 2 (6) + 3 (2) 
VI 15688 —25 (6) + 6(1) 1957 IV 20808  —32 (5) + 4 (2) 
~_ +3 VI 12°674 —18 (7) + 2(1) 
Lines poor. — 16: + 3 
Lines rather poor. Probably variable 
No. 20. Br III velocity. Range 34. Possible indications 
1956 IV 18-762 —20 (9) — g(1) of line doubling. 
1957 IV 20-701 —28 (9) + 1 (1) 
VI 9°716 =-—20(11) — 6(1) No. 27. Bors III 
— 23 = sa 1956 IV 18-862 — 24 (10) —19 (1) 
: 1957 IV 13'737, — 3(11) — 1) 
No. 21. Mo Iab(b ?) , VI 12:706 —18(12) + 5 (1) 
1956 IV 11-840 —14 (6) ~15 —5 
1957 Il 7905 —14 (7) 
IV 20°938 —11 (13) . 
VI 9862 ae (5) No. 28. Br V 
~12 1957 IV 23°735 —13:(3) + 4 (2) 
1958 III 8-860 : + 1 (1) 
No. 22. Br II + 2 
1956 IV 18-769 —11(15) — 5 (2) Individual measures scatter too widely 
1957 IV 23°838 — 8 (17) 0 (2) for a meaningful mean in the second 
VI 9°756 = 14 (11) spectrum. There are indications of line 
—II —2 doubling. Probably a binary. 


Explanation of Table III 
The observations are listed chronologically under each star. The time of each 
observation (U.T.) is given followed by the derived stellar and interstellar velocities 
(reduced to the Sun). The figures in brackets are the numbers of lines contributing 
to each velocity. The adopted means are given in heavy type followed by brief comments 
on variability etc. 
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An analysis of the velocities of the 15 B stars defined as above, shows the 
mean square deviation from the mean to be o,?=46-6(km/sec)? whilst the 
inter-agreement of the individual determinations of the velocities of these 
15 stars yields a mean square observational error of o¢?= 15-4 (km/sec). The 
difference between these two quantities is large and if interpreted as a velocity 
dispersion in the cluster (or an expansion of the cluster) yields a true mean 
square deviation o? = 31-2 (km/sec); o=5-6km/sec. However, it is necessary to 
investigate whether effects in the stars themselves can give rise to line shifts 
which might be misinterpreted as Doppler shifts due to space motion. That 
there may well be some effect of this kind is seen if the measured radial velocities 
are plotted against the values of m, (Fig. 3a). It is clear that the bright and 
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Fic. 3.—Relation of radial velocity to my for (a) stars with constant velocities, 
(6) Ca II interstellar lines. 


faint stars give more positive velocities than those of intermediate brightness. It 
will be clear that the scatter of the points about the line which has been sketched 
in is no greater than that to be expected from the observational error estimated 
above (ce=3°9km/sec). For comparison, in Fig. 3b are plotted the velocities 
determined from the interstellar Cam lines. The absence of any magnitude 
dependence in these velocities is evidence against the effect of Fig. 3a being due 
to observational causes such as guiding errors etc. If the stellar velocities are 
plotted against my rather than mp, the effect is still present though there is some 
slight indication that there is more scatter in that case. 

One of the most conspicuous features of Fig. 3 is that the velocities of two 
supergiants (No. 3 and No. 4) are considerably positive with respect to the 
giants. Now these two velocities depend mainly on the measures previously 
published (g). Whilst the general observational procedure and methods of 
reduction were the same in that case as in the present work, there were two 
changes: (1) in the present programme Kodak 103a0 plates were used for all 
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the stars as against the finer grained Ilford SBS or I.5.R. plates used in (9); 
and (2) in the case of most of the stars on the present programme a spectrum 
width of 0-42 mm was used as compared with a width of 0-20 mm for the earlier 
programme (cf. (10)). ‘Though it would appear extremely unlikely that either 
of these changes could have introduced any systematic velocity errors of the 
magnitude shown in Fig. 3, it was felt necessary to check the point. The last 
two results in the cases of stars No. 3 and No. 4 in Table III were taken under 
conditions identical with the rest of the programme (103aO plates, 0-42 mm 
spectrum width). It will be clear from Table III that if these velocities had 
alone been used essentially the same results would have been obtained. In one 
case the velocity of the star would be increased and in the other decreased. 

It is natural to enquire whether erroneous wave-lengths of some of the 
spectrum lines could produce a spurious effect of this kind. It will at once be 
obvious that since all the observed B stars in the cluster have approximately the 
same spectral types (an extreme range of B2 to Bo with most stars Bo-5 or Br) 
any effect on wave-lengths is almost certainly due to luminosity effects. Since 
it seems that luminosity effects in the B stars vary monotonically from the dwarfs 
to the supergiants it is obviously difficult to explain Fig. 3a in this way. The 
adopted wave-lengths have been discussed previously (10); in particular, 
certain corrections to the I.A.U. wave-lengths of He 4026, 4471 were applied 
following the work of Petrie (21) and Radcliffe investigations. The corrections 
were 

(1) He4026 Luminosity classes I-II : no correction 

Luminosity classes III-V: apply correction +4 km/sec 


(2) He4471 Luminosity classes I-IV : no correction 

Luminosity class V : apply correction + 10km/sec. 
The + 10km/sec correction to 4471 for the dwarfs is the most important of the 
two corrections and if it had not been applied to these stars they would be plotted 
with smaller (more negative) velocities in Fig. 3. However, so far as the effect 
we are dealing with is concerned, such changes are negligible. For the four 
dwarfs plotted (Nos. 10, 13, 14 and 18) the effect of this correction on the means 
is 1*2, II, 1-4 and 1-0km/sec respectively. This example, furthermore, makes 
clear that in order to explain Fig. 3 as due to erroneous wave-lengths very 
large shifts of a few lines, or shifts of the order of 1okm/sec for a considerable 
number of the lines, would be necessary. Neither alternative appears at all 
likely. ' 

To investigate further whether the effect could be traced to erroneous wave- 
lengths, or alternatively to a deviation of lines due to certain elements only, the 
velocities derived from individual lines were investigated. The results are 
summarized in Fig. 4. The stars were arranged in four groups: 


(1) Stars No. 3 and No. 4: mean my = 5:6 
(2) Stars No. 8 and No. 22: mean my = 6°8 
(3) Stars Nos. 20, 25,6, 19 and 16: mean my=7°5 
(4) Stars Nos. 27, 24, 14, 10 and 18: mean my)=8°5. 


Star No. 13 which is shown plotted in Fig. 3a (m)=9°55) is omitted in Fig. 4 
as it is of rather too different a magnitude and velocity to be included in the last 
group and single stars are of too low weight to be plotted separately. 
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In Fig. 4 the dependence of velocity on my is shown for the H lines, He 1 lines, 
a group of lines of Si, C and O and also the Ca 1 interstellar lines for these same 
14 stars. The numbers in the diagram are proportional to the weights of each 
point and the horizontal lines indicate the magnitude of the standard errors. 
Both H and He! reproduce the effect shown by the mean velocities. In the 
case of the Si, O and C group, which is of lower weight than H and He, the three 
brighter groups of stars show the effect though the fainter group fails to show the 
increase in velocity shown by the mean curve. However, the standard errors 
of the quantities involved indicate that there is no significant deviation from the 
mean curve. As in the case of Fig. 3 5, where more measurements are plotted 
the interstellar Call results of Fig. 4 do not show any significant dependence 
On Mp. 




































t 1 t 
6.0F 
ms 82 
7-OF 
| '22 MEAN 
6-OF 
134 
a 
1 1 1 1 1 1 1 a 4 ‘. 1 | _————— 
-25 -20 -I5S -10 -5 O -25 -20 -I§ -10 -5 Oo -10 -S O +3 
kmfsec 
s T T T T “T T i rT a 
6.0 = 
7 
Mo 
7-0 + 
MEAN Hel MEAN) J 
ol, ci 
8.0 si,si® 4. 
i2 
! ! 1 i | ! l ! ! — 
-30 -25 -20 -I5 -IO <5 O -25 -20 -i5 -0 -5 O 
km/sec 


Fic. 4.—Relation between radial velocity and my for stars arranged in groups and for lines 
of different elements. The numbers are proportional to the weights and the lengths of the 
horizontal lines to the standard errors. 


The dependence of the velocities on magnitude clearly rules out a discussion 
of the internal motions of the cluster. No convincing explanation of the 
phenomenon can be suggested and it is most desirable to confirm the effect in 
other clusters. The published data on velocities in open clusters containing 
stars of similar spectral type to NGC 3293 do not appear sufficient to carry out 
further tests. It is, however, of some interest to note that in the Pleiades a 
similar effect may perhaps be present (cf. Smith and Struve (22) Fig. 1), though, 
of course, the spectral types and luminosities of the stars are quite different in 
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NGC 3293. Further work would probably be necessary to decide conclusively 
whether the effect is definitely present in the Pleiades. 

The positive displacement of the supergiants with respect to the giants in 
NGC 3293 is reminiscent of the ‘‘ Trumpler ”’ effect (23): the positive velocities 
exhibited by some O type stars in galactic clusters with respect to the B type 
members. ‘The interpretation of this effect as a gravitational red-shift has been 
much discussed on account of the exceptionally large masses that this implies 
for the O type stars involved. It is desirable to extend the work on the clusters 
showing the Trumpler effect to fainter magnitudes to determine whether the 
fainter cluster stars also show positive residuals. It will be recalled that the 
trapezium stars of Orion show a positive velocity with respect to the nebulosity 
(24). Mass motion in the stellar atmospheres of course remains a possibility 
but it is not at all clear why it should result in a magnitude dependence of the 
form of Fig. 3a. 

It is of some interest to calculate the velocities which NGC 3293 and the 
clusters showing the Trumpler effect are expected to have due to galactic 
rotation. Using the standard value of the solar motion and the Camm linear 
approximation without the K term and adopting A=17-5km/sec/kpc, 
R,= 8-9 kpc, 4) = 327°°7 (25), we obtain the results given in Table IV. For 


TasBLe IV 


Observed velocities and velocities calculated from galactic rotation constants for NGC 3293, 
the Orion region and clusters showing Trumpler effect 


Vo Vz Vo 
Cluster Distance (calculated) (23) (23) 
NGC (pc) km/sec km/sec km/sec Vo—Vpz Vo—Vo 
2244 1660 (26) +39 +35 (6) +44 (3) + 4 ~% 
2264 870 (26) +25 +19 (11) +33 (1) + 6 —8 
2353 760 (23) +30 +20 (5) +35 (1) +10 —5 
2362 1450 (26) +43 +34 (7) +43 (1) +9 ° 
6871 1340 (23) — 8 —21 (3) —11 (2) +13 3 
7380 1840 (23) — 36 —43 (4) —40 (1) +7 +4 
Mean + 8 —2 
Orion 400 (26) +22 ‘Trapezium stars (24) = +32 km/sec 
Nebula (24) = +18 km/sec 
3293 stars 2600 + 1 Measured stellar velocities : —23 to —2 km/sec 
gas 1300? + 4 Measured Catt velocities : +0°5 km/sec 


the Trumpler clusters (23) Vy and V, are the mean velocities for the B stars 
and the O stars respectively ; the difference (V,—V,) is therefore the value of 
the Trumpler effect. The number of stars used in each mean is given in 
brackets. In the case of NGC 3293 the predicted velocity is more positive than 
that of any of the stars measured, though the difference from the supergiants 
and faint dwarfs is relatively small. The results of Table IV show that so far as 
galactic rotation alone is concerned the O stars of the Trumpler clusters have 
as much claim to represent the true cluster velocity as the B stars. A similar 
situation exists in the Orion Nebula where again the predicted velocity falls 
between that of the trapezium stars and that of the gaseous nebula. 
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In concluding this section it should be remarked that the velocities of the 
three stars in the “‘ flare ’’ of the cluster do not materially add to the discussion 
already given of their possible cluster membership. 

4. The Ca1z interstellar velocities—Listed in Table III are the velocities of 
the interstellar Cait lines. Generally only K was measured but in some cases 
H was resolved from He and could also be measured. The results yield a mean 
value of +0°5 +0°8(s.e.) km/sec. For a comparison with the stellar results we 
find o ?=18-2 (km/sec)? oe?=9g-5 (km/sec)* and therefore o?=8-7 (km/sec)? 
o=2-9km/sec. Whilst there is no dependence of the velocity on the magnitude 
of the star in which it is found there is some indication of a dependence on the 
star’s position. This is shown in Fig. 5a where the velocities are plotted in the 
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Fic. 5.—Radial velocities for (a) interstellar Ca 11, (b) stars plotted on the plane 
of the sky (cf. Fig. 1). 


positions of the stars in the sky. On the other hand, there does not appear to be 
any noticeable relation between the stellar velocity and the star’s position 
(Fig. 56). The line in Fig. 5a which roughly divides the more positive from 
the more negative velocities (omitting the three stars in the flare) is approximately 
parallel to the galactic equator. It might appear likely that any Ca 11 absorption 
from gas associated with NGC 3293 contributes only a fraction of the total line 
absorption observed since the velocities given by the stars and gas are quite 
different and, further, the difference in the velocities is of the wrong sign to be 
interpreted as an expanding cloud of gas about the cluster, as has been found in 
the case of some groups of B type stars (27). On the other hand, the mean Ca lI 
velocity agrees better with the velocity predicted from galactic rotation for 
2°6kpe than for 1-3kpe (Table IV). However, the difference is relatively small 
and does not allow us to fix the mean distance of the Call absorption with any 
accuracy. A gradual change of Call velocity with position in the sky is not 
surprising, and the effects of differential galactic rotation of clouds at different 
distances, together with the known velocity dispersion amongst individual clouds, 
will be sufficient to account for the value of o found above. 
Table V summarizes some of the results found in this work. 
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TABLE V 
Some data on NGC 3293 
Distance modulus (m)— M) = 12°08 
Distance = 2°6kpc 
Approximate diameter = 3°5 pc 


Mean velocity 
Mean Car! interstellar velocity 


—14:1+ 18 (s.e.) km/sec 
+0°5+0°8 (s.e.) km/sec. 
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ON APPARENT LUMINOSITY IN GENERAL RELATIVITY 
V. Foseph 
(Communicated by F. A. E. Pirani) 
(Received 1958 July 4) 
Summary 
A formula for apparent luminosity, valid in an arbitrary space-time, is 


derived from the conservation equations of the electromagnetic stress-energy 
tensor for a null field, by the use of optical coordinates and tetrads. 





1. Introduction.—One of the major contributions of relativity theory to 
modern physical thought is the notion of field, the most important example of 
which is the radiation field. The electromagnetic radiation field is described 
by a skew-symmetric tensor F’,, satisfying Maxwell’s equations 


oF, W, , Wa _ a 
=< + ee + eats FY" =0, (1.1) 
while the electromechanical properties of the field are expressed by the 


electromagnetic stress-energy tensor E+’ 











Ev — —g"®F Fee +g" FOF. (1.2) 
which satisfies the conservation laws 
E« =o. (1.3) 


‘ 


From (1.1) we can derive (1) a ‘‘ wave’’ equation in general relativity, for the 
propagation of the field quantities F,, : 

oF. op = ARyrael  — RF + RF os (1.4) 
the symbols having their usual meanings. 

In flat space-time (1.4) reduces to the ordinary wave equation, showing 
that the field quantities are propagated with fundamental velocity. In a non-flat 
space-time the solution of (1.4), subject to specified boundary conditions, presents 
a formidable problem (2), except in some simple cases (3). In view of this 
difficulty there have been various attempts to obtain particular solutions of the 
electromechanical equations (1.3) rather than the field equations (1.1). It has 
been shown (4) that a test particle (corpuscle) arising, through certain limiting 
processes, from the stress-energy tensor E“”, moves along a null geodesic. This 
result is in accordance with the principle of relativistic geometrical optics, usually 
expressed in the variational form: 


dine, 8 | de=o. (1.5) 


wUvae 


From quantum theory we are familiar with the idea of a photon as a particle 
whose relative energy is hv and whose 4-momentum is a null vector. In general 
relativity the formula for the Doppler shift of a radiating star can be derived 
either (a) from geometrical optics, using the quantum hypothesis for corpuscles 
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moving along null geodesics or (6) from kinematical wave optics based on the 
properties of null 3-surfaces in space-time. Kermack, McCrea and Whittaker (5) 
have shown that the two viewpoints (a) and (b) lead to the same formula. They 
have thereby established the internal consistency of the particle and wave 
descriptions of the electromagnetic radiation field in general relativity theory 
in so far as the formula for Doppler shift is concerned. 

In the present paper it is shown that the particle and wave descriptions of the 
radiation field lead to identical results also for the apparent luminosity of a 
radiating star. A general expression for apparent luminosity was first obtained 
by Walker (6) using geometrical optics, although particular cases of this expression 
had been given earlier by Milne (7) and Telman (8). Here Walker’s general 
formula is derived from the conservation laws of the electromagnetic stress-energy 
tensor for a null field, without employing the language of geometrical optics. 
This approach to the problem of luminosity seems the more natural one in a 
field theory. 

The result established here is exact and invariant. It has been found 
convenient to use the physically significant system of coordinates known as 
optical coordinates (g) but the results are independent of this choice. ‘The 
tetrad formulation is used to describe in an invariant way the local experience 
of the observers concerned. 

2. Optical coordinates and tetrads.—We begin with an arbitrary space-time 
with metric 

ds? = g,, dx“ dx” (2.1) 
of signature (1, —1, —1, —1). 





of” 


Fic. 1. 


Let C be the world line of a star in space-time (2.1). Consider a set of four 
unit orthogonal vectors ,“ at O, Aq* being time-like and tangent to C, and 
Xa)" being space-like. (Latin letters a, b,... range and sum over 1, 2,3.) A tetrad 
field along C can be generated by Fermi-propagation of the vectors i)“, thus 


(8,0—2,) ro” <0 (2.2) 
or 


where 4 is the unit tangent vector to C and 7 is the arc-distance as measured 


from O. 
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Physically, the tetrad 4,“ at any event on C is interpreted as the local 
Minkowskian frame of an observer whose 4-velocity is A“ and who uses local 
cartesian axes ,,)" (10). 

The coordinates x“ of an event P on the forward null geodesics through any 
point S on C can be written in the form: 

xt = fu(P, 2, 3, 7,2) (2.3) 
where /', /?, [3 are the direction cosines of the null geodesic SP at S with respect 
to the cartesian axes ,,)“, that is: 

Je = Pua" 

PAo" 

p* being any null vector tangent to SP at S, 7 is the proper time as measured by 

a clock moving along the world line C, and A is the projected distance of the 
null geodesic SP with respect to the 4-velocity A" of the star at S, that is 


a 
A= [doy de (2.5) 


the integrand at any M point on SP being obtained by parallel propagation of 
Aw" along the curve SP to M; if Aw" is the vector at M on SP, then its value at 
a neighbouring point N on SP is given by 


Diol = oI — PagAco ax? 
where dx is the infinitesimal displacement vector MN and the I’",,“ are Christoffel 
symbols. 

Since /4/* = 1 by (2.4), the equations (2.3) define a transformation of coordinates 
from the system x* to ‘‘ optical coordinates’’ (7) A, 7, «, B («, B being any two 
independent functions of /*) in a domain D of space-time where the forward 
null geodesics through points on C do not intersect one another. Let the 
metric of space-time in the new optical coordinate system be 


ds* =h,,, dy" dy’, (2.6) 


where y° =A, y!=7, y?=«a and y®= 8B. 
It is shown in Appendix I that 


(2.4) 


hy, =O (u#1) and hoy =1 in D. (2.7) 
That is, the metric (2.6) can be written in the form 
ds? = 2 dy® dy! + h,; dy' dy! (i,j = 1, 2, 3). (2.8) 


It should be realized that the optical coordinates A, 7, «, 8 are observable quantities 
with direct physical significance in the sense that an operational procedure can 
be specified for measuring them. They are constructed in an invariant manner 
depending only on the choice of world line C, and not on any pre-existing 
coordinate system in space-time. 7, «, 8 can be realized by the single observer 
moving along the world line C, while A can be evaluated with the cooperation of 
a set of observers stationed along the null line SP at relative rest with respect to 
the observer with world line C. 
The concept of relative velocity of two separated observers may be generalized 
conveniently to general relativity in the following way. 
Let X,,)", 1" be the coordinate tetrads of two observers S, P whose world lines 
are C and C’ (Fig. 2); the backward cone of null geodesics at an event P on C’ 
intersects C in S. Consider a time-like 4-vector A at P which is obtained 
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by parallel propagation of Ao from S along the null geodesic SP. Then the 
scalar v given by 


v=1- or" Nout < (2.9) 


evaluated at P, is the magnitude of the velocity of the observer S relative to the 
observer P. The direction cosines |, of the velocity are 


L, a No" Nano" Mo} (2. I 0) 





Fic. 2. 


The concept of ‘‘ relative rest’’, due to Kermack, McCrea and Whittaker (5) 
is a special case of the concept of relative velocity (2.9) of two separated observers. 
S is said to be at rest relative to P if his velocity relative to P is zero. 


3. Propagation of null field.—In relativity theory there are two physically 
distinguishable electromagnetic radiation fields: (i) null-field and (ii) non-null 
field (11, 12). We shall consider the propagation of null fields in general relativity. 
The stress-energy tensor E*’ for a null field can be written in the form 

Ev = Untn’, (3-1) 
where n* is a null vector and U a scalar function. It is readily seen that the sum 
of two null fields is in general not a null field. 

L. Mariot (13) has shown that nm“ generates a null geodesic in space-time. 
In optical coordinates (2.6) m“ may be taken as 

n* = [1,0,0,0], (3.2) 
since SP (Fig. 1) is a null geodesic. 

The conservation equations (1.3) for FE” read: 


(— Wy (meee + {fbr =o, (3-3) 


For a null field with stress-energy tensor E“” given by (3.1), equations 
(3-3) reduce to the single equation 


£ ((-hy!8U} =0, 


whence 


=_ A . 
U= ym (3-4) 
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where A is an arbitrary function of 7, « and B only. From (2.8) 


det (h,,,)= — det (h,,,) where p, q=2, 3 


that is 
h= lA, vl asl lh yal: 
Therefore (—A)#?2=|h,.|*2. 
Hence ? A 
sien hp" (3-5) 
| eal 


The intensity of radiation as measured by an observer P (Fig. 2) moving 
along the world line C’, is given by the scalar J, 


[= E*'N@,.N), = Un'n'no,.Nov» I= U[non]?: (3-6) 

We can express 79 in terms of two observable quantities with respect to the 

observer P: (i) the velocity v of S relative to P and (ii) the angle 6 between 
the direction of the ray PS and the direction of the relative velocity v. 

The direction ratios of the null line PS with respect to the tetrad y,)" are: 


1,.1Na)" (3 7 ) 
and the direction ratios of the relative velocity of S with respect to 7)“ are 
— Xo" Nadu (3 8 ) 
Thus the angle @ between these two directions is given by 
COS A= — 0, Na"Ao)’ Nap LMM), Ney NOx Mo Nore)” (3-9) 


This expression may easily be rewritten in terms of v and yop» only 
(Appendix II). On solving the resulting equation (A2.7) for yo» one obtains 





Noo = (1 — v?)(1 + vcos 8). (3.10) 
The expression (3.6) for the intensity of radiation becomes 
~y2 
I £ a (3.11) 


[hy glt’? (1 +008 8)?’ 

Thus we have established a formula for the apparent luminosity (intensity 
of radiation) of a radiating body S in terms of observable quantities A, 7, «, B, v 
and 6. The parameter A depends on the radiating body only. 

This result has been found without the introduction of geometrical optics 
or entities like photons and the relation E=hv. 

We shall now show that this result agrees with that obtained by Walker (6) 


from geometrical optics. 

The cross-sectional area V of a thin pencil of null geodesics issuing from the 
radiating body at S (Fig. 2) and meeting the instantaneous space of P is, for 
the metric (2.8), given by 

V =|h,,q|''? da dB. (3-12) 
Let p“ be the momentum of 4-vector of a photon whose world line is SP; then 
the Doppler sh.ft D is 
— frequency of photon relative to 
frequency of photon relative to 7)’ 

_ energy of photon relative to 

~ energy of photon relative to y"’ 
by the quantum hypothesis. 
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Therefore 


p= Pron | 
P*N0 
But 
p" =hv[1, 0, 0, o], 
and 


Ao, = [1,1,0,0] at S, 
where v is the frequency of the photon relative to Aw,. Hence 
D=[noo)-. (3-13) 
From (3.5) and (3.6), we obtain 
I= B/V D~ 
by (3.12) and (3.13). Thisis, in fact, result (1) in Walker’s paper (6), B depending 
on the source only. 

The derivation of J as presented here confirms the validity of the corpuscular 
description of the radiation field in the case of apparent luminosity. 

4. Application to special theory.—We shall illustrate the use of the formula 
(3-11) in the special relativity case when the world line of the radiating star is 
a geodesic in Minkowskian space-time. With the world line C of the star as 
base line, the space-time metric can be written in optical coordinates as 


ds? = dr? + 2dr dX—d*(da2 + sin? « dp?). (4.1) 


The intensity of radiation as measured by an observer P (Fig. 2) is given 
by (3.14) as 
A 1—v* E I—v 
> > 5 tre, (4.2) 
sin? a (1+vcos6)? A*(1+vcos6)? 
where « depends on the source only, and J is the projected distance of the null 
geodesic SP with respect to the 4-velocity of the radiating star. In order to 
relate J to the observer’s notion of distance, it is convenient to give the result 


in terms of the projected distance » of SP with respect to the 4-velocity of the 
observer P. Thus 


P dy# 
= | WW) dy *= I U1) dX dv= Muy mI dd, 
Ss 


= Anjan =A(1 — 08)}!2(1 +0 608 8) (4.3) 
Therefore (4.2) reads 





i.e. 


\ [=en-*(1 —v?)?(1 + vc0s 8). 
Putting v=tanhu we obtain 
I=ep~(coshu + cos @sinhu)~*. (4.4) 
This result agrees with that obtained by Robertson (14). However, 
Robertson’s treatment of the subject involved the use of a particular coordinate 
system (rectangular cartesian) in extended regions of space-time, the formulation 
apparently being non-invariant. 
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APPENDIX I 


Form of metric in optical coordinates.—We now derive equation (2.7) 


ho, =0 (u#1), ho =1 
for the metric in optical coordinates (2.6) 
ds* =h,,, dy" dy”. 

The proof is similar to that given in Temple’s paper (9). 

Since the equation of C is given by A=0 we have, for a displacement along C, 

h,dy'dyi=dr? if A=o [i,j = 1, 2, 3]. 
That is 
hy =1 , 
\ if A=o (At.1) 


=0, hy, =0 


hog 
where p, q assume the values 2, 3. 
For a displacement along any null geodesic at S, r=const., 
h,,, dy" dy” = hg dd. 
That is 
hyy=0 if A=o. (A 1.2) 
The equations of the null geodesics from S, namely y/=const., y°=A, must 
satisfy the differential equations 


dye r ay" dy? _ 








ae tT ex ax =? (Ar-3) 
and 
h,, dy dy’ =o. 
The parameter \ has been chosen.in virtue of (2.5). It follows that 
Too" =0 and hog =0 in D. (A 1.4) 
That is, 
hyo — Moo _ _ 
2—* oy" = 2D un, y= 
Therefore by (A 1.4) we have 
Ohio 
a =O im D. 
Hence hy is independent of A } (A 1.5) 
and hyyp=o0 by (At.2) 


By definition of A (2.5) 
P P dy” P 
A= [ dondyt= [doth a—k [dd 
or s ° a 8 
since k=Nq""h,,(dy’/dA) is constant along SP in virtue of 
Oo? 9 and - (5) = 0, 








dA SA\ AA 
Thus 
Aw*h,o= 1 along SP. 
At S, 
No =8," by (A1.1), 
and so; 


hyy=1 at S. 
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Hence by (A1.5) 


hyy=1 in D. (A 1.6) 
Relations (2.7) follow from (A1.4), (A1.5) and (A1.6). 


APPENDIX II 
We shall derive the relation 


1 = Moot - 08)-!9 


08 l(t = 8) TPE 





Introducing the notation 
n™ = eM n19),5 (A 2.1 ) 
where 
e*8 = diag {1, —1, —1, — 1}, 
we can write (3.9) in the convenient form 


COS O= 1, Nq)"Awy?n,[{ — Menyn7N, }{ — AON" Ao)’, O}] 2. (A 2.2) 


Now 
1Na"Aorn = 1a)" Ao)” =r 1,1? Aoy? = 1,707, Ao) « 

But 

Mo" = 54, (A 2.3 ) 
and so 

_ f 

1,.Ma)" Aon,” = nA? o {nnn Aw? . 

Also 


—— Siet = 
M(H = NN HN ON = N,N), ON" — NHN, On =n,n? — [n, no}. (A2.3) 
Therefore 


b) — 
nmoynn, = — [nos 
since n“ is null. Similarly 


Noreen’, © = Ao”Aoy, — [An]? = 1 — [orn }?- 
Thus 


cos 4 = [n, Aw? — (2, n0)")(MoAw’ )I[Awn)? — 1]-?(a,.n0")- = (A2.4) 


From (2.9) we have 
100) Nop = (1-2?) “*._ (A2.5) 
Since A? is propagated parallel with SP and 5n°/SA=o0, n,Ao? is constant along 
SP. Therefore 
[2,do lp = [Ao ls = 1; (A2.6) 
since hy, =1 at S by (A1.1). 
Substituting (A2.5) and (A2.6) in (A2.4), we obtain 
cos 8= [1 (non, )(1— 02)! ][(1 = 08)-1— 1] (09) 
that is, 
I — Noo(1 — v?)-*? 


cos 6 = 
Tol (= 8) 1? 





(A2.7) 


by (3.2). 
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PE3IOME JJOKJIAJIOB, B ITEPEBOJIE HA PYCCKHM A3SbIK 
PACCEAHHE PAXMOBOJIH B COJIHEUHOM KOPOHE 


A. Xowum 


M3 naOmogennii paguousnyuenua KpaSopyyqHoii TyMaHHOCTH B HIOHe, Kora OHa HaXOAMTCA 
Ha HeOOJUIOM yrOBOoM paccTosHHM oT Cosma, CilefyeT, 4TO u3JIyueHHe Mo_BepraeTca 
3HAYHTEJIBHOMY pacCeHHHIO MPH MpoOxoxKeHHM uepes coHeUHYyIO KopoHy. B yaHHo cTaTEe 
OMMCbIBaIOTCA e?KeTOTHbIe M3MepeHHA 3a Mepwom c 1952 m0 1958rr. Bim momyueHpr 
cyleqtyrouywe HOBbIe pe3yJIbTaTHI: 

(1) u3MepeHHA Ha TJIMHe BOJIH B 7,9 M YKa3bIBaIOT Ha HasIMUuHe petbpakKUMOHHbIx adpdekTos, 
CBepx pacCeAtHHA, KOTOPble MOTyT aTb CBeeHHA O CpeHeii IIeEKTPOHHOM MIOTHOCTH B KOpoHe; 

(2) pacmpoctpaHeHve u3MepeHHii Ha Goslee KOpOTKy!O WJIMHY BOJIH B I,9 M NosBosmsO0 Gonee 
TOUHO ONpeyesMTb Gombuioe paccesHve Ha paCcCTOAHHM HaMOobUIero NpHOMKeHHA ; 

(3) HaOmofanoch OTUeTIMBOe M3MeHeHHe B Te4YeHHe WMKJIa COJIHCUHbIX MATeH B palioHe 
KOPOHBI Ha paccToAHHH 8RQO; u 

(4) paccesMe aHH30TpOMHO. OTOT pe3yibTaT NOATBep»KyaeT MperkHee 3akIOUeHHe, U4TO 
pacceAsHwe MpOMSBOAMTCA HUTeCBHHbIMH HelpaBHJIbHOCTAMH, BbIpaBHCHHbIMH B MarHHTHOM 


moje H yKa3biBaeT Ha TO, 4YTO HallpaBJeHHe NMOJA MpHONMSHTeEIbHO paMasIbHO Ha pacCTOAHHH B 
1520Ro. 


M3MEPEHHE IIPO@MJIA XPEBTA CKJIAZKM HA JIYHHOM MOBEPXHOCTU 
Tuavbepm Puavdep 


XpeGeT ckaqKH Ha JIYHHOM MOBepXHOCTH M3MePAJICA MeTOJIOM, COCTOALIMM M3 H3MepeHHit 
MeHAIOUIelicA JIMHbI TeHH, Opocaemoii Ha Hero coceyqHei ropoi. BepoaTubre ommoKn 
MeCTOHAXOPKeHHA TOUeK +330MB MJIOCKOCTH MOBepxXHOCTH HM TOJIbKO +14 MB HallpaBJICHHH 
NepneHMKyJIApHOM moBepxHocTu. Ilo mpodumto xpeOta cKIaqKH BHHO, 4YTO OH He MO?KeT 
CU4HTATLCA rpeOHEM BOJIHbI 3aTBepAeBIUIeH WaBbl. 


YCKOPEHHbIM METOJ,L OBHAPYXKEHHA HEMTYHA 
P. A. JIummaomon 


OnvcaH MaTemMaTMuecKHii MpHHUMM MeTOga oOOHapyxxkeHua HentyHa, Hykyalouleroca B 
3HaUMTeJIbBHO MEHBILIEM UMCICHHOM BbIUMCIeHHH. Bpema remMOleHTpHuecKorOo coeqMHeHHA 
MO?KET ObITh HaiieHO HMCKIIOUNTeIbHO pa3sOopom pacxooKaeHHi gonrotTh: Ypana. Oyun oTu 
CBeJIeHHA alOT BOSMOXKHOCTh LpewBuyeTh DoNrotTy, coryacHo sakoHy Bboyja, c TOUHOCTbIO B 
mpeyemax 15°. Haiigqa paccrosHue Ha KpyroBoii opOuTe, HanOomee cormacyroulelica c HaOmOye- 
HMAMH, MO?*KHO MIpes{CKasaTb 9Ty AONTOTy C TOUHOCTbIO CpaBHHMOM C TOUHOCTHIO JOCTHTHYyTON 
JIe Beppse. Ykasbipaetca Kak, OnaroyapA noy_xo_AUemMy KOMOMHMpOBaHHIO = yCJIOBHbIX 
ypaBHeHHii, UHCIIO HeEH3BeECTHbIX MO)KET OBITh DOBeeHO JO Tpex (M0 CpaBHeHHIO C BOCbMbIO B 
MlepBOHaUaJIBHbIX MeTOMax) DJIA mOGoro pacueTHOro cpefHero paccTtonHua. ToT xe mporecc 
ycTpaHAeT HeEKOTOpbie HeyfOOHbIe OCOOeHHOCTH, KOTOpbie HHaue BOSHHKIM ObI AA OpOuTHI 
O1M3b pesoHaHca 2: I. ‘ 


SBESHBIE PVYIIIIbl. II. QBMKEHMA B IIPOCTPAHCTBE 
HECKOJIbBKUX 3BE3][-KAPJIMNKOB PIABHOM MOCJIEZOBATEJIBHOCTU 


Oaunb Jzeen 


IIpH MoMOUH 3Be3QHbIX BeJIMUHH HW UBeTOB, MOUYYeHHbIX oTosIeKTpHuecKuM cmocobom, 
olpesesieHbi cboTOMeTpHuecKHe MapaswiaKCbl DJIA 359 3Be3], C MSBeECTHbIMH JIYYeEBbIMH CKOPOCTAMH 
M cCOOCTBeCHHbIMH JBHKCHHAMH, ABJIAIOUIMXCA KapJIMKaMH TyiaBHO NMOCIeOBaTeNbHOCTH 0 
cieKTpockommueckoH Klaccudpukalyuu oOcepBatopuii Mayxt Busscona unum EpKeca. Om 
MapasiaKCbI MPHMeHAIOTCA JIA BbIMMCIeHHA JBWKeHHM 3Be3, B MpocTpaHcTBe MO BeKTOpam 
U, V u W. 3Be3qbI mompasyenAIoTCA cormlacHoO BeJIMUMHe JBMWKeHHA nO W (nepneHqukKysAp- 
Horo K miocKkocTH TamakTuKH), a TakK>Ke 10 HX JBH)KeHHIO, MOCTpoeHHOMY B mIocKocTtu U,V B 
PHcyHKax 2 yo 5. W3 oTnx AMarpamMM BBbIBOAHTCA, 4TO pacipeyseseHHe ABMWKeHHA B MpocT- 
PaHCTBe He CJlyualiHO, NO KpaiiHeli Mepe DJIA MayIbIx BenmunH W 
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YPOBHH SHEPrMN WU BEPOATHOCTH MEPEXOJIA IUIA FE iv 
P. I’. Tapcmane 


PacuuTaHbl KBaHTOBOMexaHWueckHe MapaMeTpbI KOHdurypauun 3d° cnektpa Fe Iv u BbIUH- 
CJI€¢HbI yPOBHM 9HepruM 9sTO KOH(burypayuu. BpruncieHbI BepOATHOCTH Mepexoya AIA M3IIy- 
YeHHA MarHHTHbIX JMMOei WH 9IeKTPH4eCKHX KBapynosei WA Mepexoq0B MexkKy ypoBHAMH 
KoHdburypauuu 3d°. Paspickupaiotca suHuM [Fe Iv] Ha cnextpe RR Telu mpegnararotca 
HEKOTOPble OTO?KeECTBIICHHA. 


@POTOSJIEKTPHUECKAA KPHBAA BJIECKA RT COSBE3IZMA CKYJIBIITOPA 
I. I. Cuasce u E. M. JTundcet 


TlonyueHbl, “cnoub3yA 150 cM’ OBbIii pedbuexTop Pokdennepa Boiigenckol oOcepBaTopunu, 
cboToseKkTpHuecKHe KpHBble OylecKa 3aTMeHHOM MepemeHHOt RT Ckxyssntopa, rommyOoro 
M >KesITOrO yBeTa. O6a MaKCHMyMa He paBHbI HM He ABJIAIOTCA MIOCKHM rpeOHeM Kakoro-s1H00 
Buga. B TeueHve nmoceqHMx 60 NeT NepHoy sTOM 3Be3MbI CilerKa yMeHbIUIAaJICA HM, IPHHHMaA ITO 
YMeHbUIeHHe 3a JIMHeliHOe, HAXONMM Cllefyroulee BbIparxKeHHe KPUBO OsecKa, HaHOosIee Corlacy1- 
ujeeca c HaOIOeHHAMH: 


Min. = JD 24237369°51145 +(04'51156702— 10568 x 107 "E)E. 


PATMOU3SJIYYEHHE H3 OBJIACTH OPHOHA 
X. Puw6em 


IlpousBeyqeHa pasBeaqkKa oOacTH OpvoHa Ha DBYyX JJIMHaX BOJIH Mp NMOMOLWIM pasMoTesecKo- 
nop B CwaHeii. PaguonusiyueHue U3 APKHX TyMaHHOCTeli B 93TOM OOnacTH NOBHHMOMy 
TepMHYeCKOrO MpOvcxoo+KTeHHA; HO cyaObiii pacTAHHyTbIii HCTOUHUK B TOM »Ke OOACTH MODKeT 
M3JIyuaTb He TEPMHYeECKHM OOpa3sOM H MO?KeET ObITb CBA3AH C UIMPOKO TeMH 2M TyMaHHOCTbWO. 
pase vege TIpH NOMOMM Mpocroi MOfesIM, paqHo u onTHuecKkHe HaOJIOAeHHA TyMaHHOCTH 

pHoua M42. 


TOUHBIE TIOJIOKEHHMA 17 CHJIBHbIX PATMOSBE3]L 
B. Davcmop 


JlatoTcaA HeCKOJIbKO HaOJIOZeHHii Ha [JIMHe BOJIH B 1,9 M, MMCIOLIMX L{eJIbIO yJIyULueHHe TOUHOCTH 
CKJIOH€HHii MCTOUHHKOB, OCHOBaHbIX Ha MepHoslax NpoxoKeHHuA Nepey, AByMA HHTepdepomeTpamu 
C HaKJOHHbIMM OCcHMH. OmnMcbiBaeTCA Takoe MpHMeHeHHe KOMOpHJPKCKOrO pajy{MoTesecKona, 
oOpasyrolsero HHTepdbepomeTp cO CKpelMBaIOWIMMHCA OCAMM HM J{alOTCA TOUHbIe MOMNOKCHHA 
17-TH HMCTOUHHKOB. JIA cpaBHeHHA MIpHBeeHbI MosOxKeHHA, OMpeseeHHbIe TpeTbeii KIMO- 
PHIPKCKOM ChEMKOM. 


JIVYUEBbIE CKOPOCTH ®YHTAMEHTAJIBHBIX IOXKKHbIX 3BE3]I 
Buascam Bacxom6é u Ilameaa M. Moppuc 


JlaioTcA BHOBb ONpeseJICHHbie JIYYeBbIe CKOPOCTH 12-TH 3Be3/, C HOPMA@JIbHbIMM CKOPOCTAMH, 
a Tak7Ke 63-X APKMX 3Be3jl, HAXOJAUIMXCA 10>KHee 9KBaTOpa. Jucnepcua CHMMKOB, MOJIYYeCHHbIX 
IIpH MOMOLIM TpexlIpusMaTHYecKoro cneKTporpada Kacrpsitua, 36 A/mm upu Hy. Ckopoctu, 
HalijjeHHble OOcepBaTopHeit CrpomJio, JI 3Be3/, HOPMAJIbBHOM CKOPOCTH OKa3aJIMCb Y/IOBJICTBO- 
PAIOIIMMH M3MepeHHAM B Ypyrux oOcepBaTopuHaAx cuctembI JIuKa. Pas0upaloTcA OcTaTOUHbIe 
CKOPOCTH]VJIA HHOIMBUyasIbHbIX CMeKTpasIbHbIx JIMHMI1. 

Jina uccneyqoBaHua ObiIM BbIOpaHbI 3Be34bI M3 MO3MYMOHHbIxX KatTasioroB FK3 u N30. 
Msmepanucb mo Kpaitueii Mepe TpH CIIeKTporpaMMbl [WIA Kako 3Be3yqbI1. T[Ipupefen ksacc 
CBeTHMOCTH 10 MepecMoTpeHHOit cucteme Epkec JIA 3Be3, paHHero Tuma. J]atoTca 29 HOBBIX 
MIOCTOAHHBIX CKOpoctei, a TakxKe 18 CKOpocTeii, MepeMeHHOCTb KOTOPbIx OOHapyrKeHa BIIepBbIe. 
BosoOHoBJIeHbI HaOOeHHA Hay, 15 3Be3aMH, CKOPOCTH KOTOPbIX 3aBeJJOMO IIepeMeCHHBI. 
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CMNEKTPAJIBHbIE THIIbI U JIVUEBLBIE CKOPOCTHM B FAJIAKTHYUECKOM 
CKOIIJIEHHU NGC 3293 


M. B. ®ucm 


112 CN€KTPOB 35-TH 3Be3, NOCNy+KHIIM JIA OMpeyeueHHA MX CICKTPasIbHbIX THMOB H 89 Tex 
7Ke CIEKTPOB—JIA OMpeesieHuA JIyueBbIX CKOpocteii 26-TH M3 9THX 3Be3Q. OTH CleKTpasIbHble 
THMIbI, COBMe€CTHO C (OTOMeETPH¥UECKHMH M3MepeHHAMN Y-pa T. 9. Xoyka, MpHBOAAT K QuMar- 
pamme Tepuumpynra-Peccena, uero cieqoBasio OKMAaTb WIA MONOM0TO CKOMNeHHA. SBesyEI 
IIpHHayiexKaT, Goble 4acTbIO, K paHHemy TuMy B c KuaccamMH cBeTHMocTH OT V go Ib. SBespy 
Mo, Iab mMo>KHO C yBepeHHOCTbIO BKJIIOUNTE B CKomIeHve. OOHapy2KeHa Oobulad MW HEOObACHHMaA 
3aBHCHMOCTS JIYYeBbIX CKOpOcTeli OT 3Be3QHOM BeJIMUHHBI (pic. 3). OTO ABJIeHHe He BIIMAeT Ha 
Me@*Ky3Be3qHbIe MHHM Ca IT wu BO3MO>KHO CBA3aHO C acpcspeKTOM Tpammuiepa. 


O BHXMMOM CBETHMOCTH B OBLIEM TEOPHH OTHOCHTEJIBHOCTH 
B. owoseh 
M3 ypaBHeHHii COxpaHeHHA TeH30pa 3JIEKTPOMarHHTHOM S9HeEprHv-HallpAyKeHHA DIA HyIeBOrO 


NOJIA, IPH MOMOWIM ONTHYeCCKHX KOOPAMHAT HM TeTpayOB, BbIBOAHTCA *opMyJla WA BUAMMOli 
CBeTHMOCTH, IIpHMeHHMaA B MIPOH3BOJILHOM MpocTpaHcTBe-BpeMeHH. 
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